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* 


SUMMARY 

This report describes a Phase A performance-analysis computer program, 
HILTOP, that has been developed explicitly to generate optimum electric propul- 
sion trajectory data for missions of interest in the exploration of the solar system. 
HILTOP is a double-precision, FORTRAN IV, EBM 360 production program which 
is primarily designed to evaluate the performance capabilities of electric propul- 
sion systems and which may, in the hands of a skilled analyst, perform efficiently 
in the simulation of a wide variety of interplanetary missions. HILTOP uses nu- 
merical integration of the two-body, three-dimensional equations o£ motion and 
the Euler- Lagrange equations. It contains transversality conditions which permit 
the rapid generation of converged maximum-payload trajectory data, and allows 
the optimization of numerous other performance indices for which no transversa- 
lity conditions are included. In addition to optimizing the thrust direction and on- 
off switch times, other significant performance parameters that can be optimized 
are jet exhaust speed, power level, hyperbolic excess speeds, launch asymptote 
geocentric declination, flight time and launch date. The ability to simulate con- 
strained optimum solutions, including trajectories having specified propulsion time 
and constant thrust cone angle, are also optionally available. The program is de- 
signed to handle multiple -target missions with various types of encounters, such as 
rendezvous, stopover, orbital capture, and flyby. Performance requirements for 
a variety of launch vehicles may be determined. The documentation includes problem 
formulation, program usage specifications, sample problems, and detailed subroutine 
descriptions. 


iii 


Page 


intentionally Left Blank 



TABLE OF CONTENTS 


Summary iii 

Nomenclature ix 

I. Introduction 1 

II. Formulation . 7 

A. Spacecraft and Trajectory Models 7 

1. Spacecraft Mass Components 7 

2. Electric Propulsion System 10 

3. Target Ephemerides 13 

4. Differential Equations 13 

5. Optimality Conditions 16 

6. Integration 17 

7. Units 21 

B. Boundary Value Problem 23 

1. Boundary Constraints 23 

2. Initial Conditions for State and Adjoint Equations 24 

3. Target Conditions and Spacecraft Constraints 25 

4. Transversality Conditions 33 

5. Partial Derivatives 48 

C. Extension of Solution for Power Degradation 49 

D. Auxiliary Computations 57 

1. Standard Block Print Variables. 57 

2. Extremum Point Summary Print 62 

3. Swingby Continuation Analysis 66 

4. Spiral Capture 70 

E. Special Program Features 73 

1. Perturbation Step Size Selector 73 

2. Avoiding Corners in the Propulsion-Time Function. ... 74 

3. Lagrange Multiplier Scaling 75 

4. Generation of an ASTEA Tape 75 

5. Scratch Pack Output 75 

6. Normal Run Termination 76 

7. Ballistic Trajectory Option 76 

8. Rotation of Primer Vector with Launch Date 77 

9. Limitation of Power to Power Processors 78 

10. Housekeeping Power 78 

11. Imposed Coast Phases 79 


v 


in. Program Description 81 

A . Program Structure 82 

1. Main Program 82 

2. Subroutine Glossary 83 

3. Subroutine Calling Sequence 88 

4. Common Array Information 93 

B. HILTOP Input 135 

1. NAMELIST 135 

2. Definitions of Input Parameters ..................... 136 

3. Default Values of Input Parameters 156 

C. Program Output 158 

D. Job Control for HILTOP 167 

1. Program Execution 167 

2. HILTOP Program Machine Requirements 167 

IV. Sample Problems and Results 169 

A. Mercury Orb iter 171 

B. Ceres Rendezvous 183 

C. Deep Space Probe ". . 201 

D. Jupiter Flyby with Ballistic Swingby Continuation 211 

E. Encke Rendezvous with Double Asteroid Flyby 231 

F. Extra Ecliptic Mission 244 

G. Comet Rendezvous Mission 257 

H. Multiple Ballistic Swingby Mission 269 

V. References 287 

VI. Subroutine Descriptions 289 

AEINWT 291 

ALBEDO * 301 

ANSTEP . 305 

BEGIN 315 

CARKEP 329 

CDERIV 335 

CHECK 351 

CONVER 383 

CONVRT 389 

CORNER 393 

DATE1 403 

Vi 


Subroutine Descriptions (eont) 

DECLIN 

DERIV 

EFM 

EFMPRT 

ETA 

EXTAB 

FINISH . . 

' FUNCT 

GETI 

GETQ 

GETRV 

GUESS 

GUNTHR 

IMPRNT 

IMPULS 

INCOND 

INPUT 

INTERP 

LOAD 

MAIN 

MINMX3 

MORE 

OMASS 

PARINC 

PDATE 

PMPINT 

PRINT 

PRINTR 

PRIOR 

PUNCH 

QPRINT 

QSTART 

RADAR 

REMTIM 

RETINJ 

RKSTEP 

SCOMP 

SETUP 

SMQINT 

SOLAR 

SPRINT 

STEP 

STORE 

SUMMRY 


407 

411 

423 

441 

451 

455 

467 

475 

487 

495 

539 

545 

555 

561 

565 

571 

577 

593 

607 

617 

623 

653 

669 

683 

693 

699 

703 

723 

.745 

753 

769 

805 

837 

847 

849 

861 

871 

877 

901 

909 

941 

963 

967 

977 


VI 


Subroutine Descriptions (cont) 


SWING 983 

SWSTO .. 1009 

SWTRAJ 1017 

TAP 1023 

TAPSET 1049 

THANGD 1057 

TIKTOK . . 1067 

TRAJ 1075 

TRAJI 1105 

TRAVEL 1119 

TWINKL .... 1123 

VMAG 1127 

VPRINT 1131 

VSCAL 1137 


vill 



NOMENCLATURE 


Generally, upper-case symbols denote vectors and lower-case symbols 
denote scalars. Lower-case symbols with bars denote unit vectors. The abbre- 
viations EPS for electric propulsion system and BVP for boundary Rvalue problem 
are used. 



EPS instantaneous thrust acceleration; semi-major axis 
Semi-major axis of primary-target capture orbit 
Solar power law coefficients 

Arbitrary unit vectors used in (132) and (139) 



A coefficient in the efficiency law 

► Launch vehicle coefficients 

Vector constant of optimal rocket problem, expression (63) 
Radians-to-degrees conversion factor 

EPS jet exhaust speed (constant); abbreviation for cosine function 

Retro stage jet exhaust speed 

Auxiliary quantity given by expression (74) 

1 Coefficients in quadratic expression for Av. , expression (78). 


d A coefficient in the efficiency law; an auziliary quantity in the coast-phase 
solution; solar flux density 

E Eccentric anomaly (a scalar) 
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A coefficient in the efficiency law; the base of the natural logarithms; 
eccentricity; subscript denoting Earth 

Spacecraft unit angular momentum vector 

Spacecraft unit radius vector 

EPS unit thrust vector 

Spacecraft unit velocity vector 

Retro stage characteristic speed exponential factor given by expression (76) 
Unit primer vector 

Auxiliary scalar function defined by (215) 

EPS instantaneous thrust magnitude ; f-function of the f and g series ; 
subscript denoting a desired value; true anomaly; auxiliary variable defined 
by equation (147) 

Retro stage thrust magnitude 

Auxiliary quantity given by expression (77) 

Auxiliary scalar functions in the coast-phase solution, equation (45) 

EPS reference thrust acceleration; g-function of the f and g series; 

BVP point-constraint geometric mean of the weighting factors 

Auxiliary quantity given by expression (97) 

Spacecraft angular momentum vector 

Magnitude of spacecraft angular momentum vector 

Spacecraft unit angular momentum vector ~ 

Variational Hamiltonian 

Cartesian components of spacecraft angular momentum vector 



max 

J 

j 

j ps 


drop 


rt 
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samp 


k t 

L 
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Thrust-switching step-function 

Subscript pertaining to an intermediate target; inclination to ecliptic; 
general subscript or running index; inclination of parking orbit about 
Earth 

Unit vector along x-axis 

Parking orbit inclination associated with range safety limit 

Index-set of the BVP dependent variables 

Unit vector along y-axis 

Unspecified-reference-power indicator 

EPS propulsion system jettison indicator (retro maneuver) 

Retro stage existence indicator 

EPS tankage jettison indicator (retro maneuver) 

Arbitrary positive constant associated with performance index; temporary 
variable ultimately equated to inverse of. the characteristic degradation 
time 

Unit vector along z-axis 

Auxiliary quantity given by expression (75) 

Intermediate -target drop-mass factor defined by expression (6) 

Retro stage tankage mass factor defined by expression (11) 

EPS structure mass factor defined by expression (8) 

Intermediate-target sample-mass factor defined by expression (6) 

EPS tankage mass factor defined by expression (7) 

Launch site latitude (scalar) 

Mean anomaly (scalar) 
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drop 
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net 
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ps 
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m 


rp 


m 


rs 


m 


rst 


m 


m 


samp 


m. 


Am 


Coefficients used in computing nuclear and total magnitudes of a celestial 
body (scalars) 

Nuclear magnitude (scalar) 

Total magnitude (scalar) 

Spacecraft total mass variable 

Auxiliary unit vector given by expression (32) 

Intermediate-target drop-mass given by expression (6) 

Net spacecraft mass 

Initial spacecraft mass (payload of launch vehicle) given by expression (2) 
EPS propellant mass 

Electric propulsion system mass given by expression (4) 

Retro stage mass 

Retro stage propellant mass given by expression (9) 

Retro stage structure mass defined by expression (11) 

Retro stage structure and tankage mass given by expression (11) 

EPS structure mass 

Intermediate -target sample-mass given by expression (6) 

EPS tankage mass 

Propellant mass increment due to primary-target spiral maneuver 

Exponent in step-size law, expression (39); subscript denoting time at 
the primary target; number of BVP dependent variables 
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Unit vector normal to the solar arrays 

Unit vector directed along a planet's north pole 

Subscript denoting launch time; subscript denoting the beginning of a 
computation step 

A celestial body's position vector; BVP partial derivative matrix 

EPS instantaneous power; subscript denoting a perturbed, or neighboring, 
parameter; auxiliary variable in equations (79) 

Ratio of housekeeping to reference power, p^/p ^ 

Total instantaneous power developed by arrays 

Housekeeping power 

EPS reference power 

Auxiliary quantities in coast-phase solution, expressions (54) and (55) 

Auxiliary variable in equations (79); solar array radiation damage factor 
Spacecraft position vector 
Magnitude of R 

Primary-target capture-orbit apocenter distance 
Earth-to-space craft communication distance 
Unit vector along line of ascending node 

Primary-target capture-orbit pericenter distance; primary-target swingby 
passage-distance 

Swingby passage-distance unit vector 

Value of r for which y-curve is at a maximum 

Abbreviation for sine function; auxiliary variable used in equations (79); 
degradation time 


xtli 



s Unit vector directed toward Canopus 

t Time 

Retro maneuver burn time given by expression (12) 

At Time-increment due to primary-target spiral maneuver 

u Generalized universal anomaly during thrust phases 

Au Generalized universal anomaly increment, equivalent to the computation 
step-size during numerical integration 

v Magnitude of spacecraft velocity 

v^ Characteristic speed of a rocket maneuver 

v Escape speed from launch parking orbit 

6 

Minimum velocity impulse required for non-coplanar injection from a 
circular orbit to a given excess velocity 

Speed of a spacecraft in a circular orbit 

Planetocentric speed at primary-target swingby closest-approach point; 
auxiliary speed given by equation (72) 

Hyperbolic excess velocity (or encounter velocity) 

Swingby planet arrival hyperbolic excess velocity 

Swingby planet departure hyperbolic excess velocity 

v b Hyperbolic excess speed (or encounter speed) 

Av Retro stage impulsive velocity increment magnitude; characteristic velocity 
associated with primary-target spiral maneuver; incremental speed required 
at powered swingby 


g 


°°A 


*D 


Av* Retro stage total velocity increment magnitude 

Av Minimum incremental velocity (magnitude) for coplanar boost out of 
circular orbit 

Av Velocity penalty due to noncoplanar boost out of circular orbit 

g 


xlv 


Av. 

1 

w 

X 


y 



a 


Ap 

y 


Velocity penalty due to launch azimuth 
Auxiliary variable in equations (79) 

First Cartesian component of position; a general variable; a general 
state variable; auxiliary variable in equations (79) 

Second Cartesian component of position; auxiliary variable in equations (79) 
Third Cartesian component of position 

EPS specific mass; geocentric right ascension of launch excess velocity 

Auxiliary parameters defined by equations (211) and (212) 

Specific mass of the solar arrays 
Communication angle (Sun-Earth-spacecraft) 

Specific mass of the power conditioning and thruster subsystem 

Arbitrary, independent angles defining orientation of excess velocity in 
( 132 ) and (139) 

Independent variable of coast-phase solution, also generalized to be the 
independent variable on the entire trajectory 

Value of P at the beginning of a computation step 

Computation step size (increment of trajectory independent variable) 

Normalized power function 


y* dy/dr 
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By/Bd, where d is the solar flux density 

Launch hyperbolic-excess-velocity asymptote declination; BVP dependent- 
variable tolerance 

Bend angles of hyperbolic arrival and departure trajectories, expression (213) 


xv 


Total bend angle given by expression (214) 

Kronecker delta function 

Auxiliary quantity in the coast-phase solution; obliquity of the Earth 
equator to the ecliptic 

EPS efficiency 

drj/dc 

In-plane thrust angle 

< 

Travel angle increment 
Travel angle 

Primer vector (adjoint to spacecraft velocity) 

Magnitude of A; a general adjoint variable; the iterator inhibitor 

Adjoint variable associated with jet exhaust speed 

Adjoint variable associated with reference thrust acceleration 

Adjoint variable associated with degradation time 

Thrust cone angle Lagrange multiplier 

Adjoint variable associated with mass ratio 

Adjoint variable associated with propulsion time 

Adjoint variable associated with thrust cone angle 

Gravitational constant of the sun; a general gravitational constant 

Gravitational constant of the primary target 

Mass ratio 

Mass ratio increment at an intermediate target 
Performance index; ratio of circle circumference to diameter 
Partial derivative of it with respect to arbitrary variable x. 



p Auxiliary variable used in equations (79) 
cr Thrust switch function 

CT* Special form of thrust switch function, given by equation (186) 

CT Portion of total thrust switch function, given by (193) 

r 

Act Propulsion-corner-proximity tolerance-interval 

T EPS propulsion time 

T, Characteristic degradation time 

d 

$ Transformation matrix for rotating from ecliptic to equatorial coordinate 
system 

0 Thrust cone angle (between thrust and radius) 

X Angle between normal to solar arrays and the spacecraft-sunline 
ij) Out-of -plane thrust angle 

Si Longitude of ascending node of an orbit 

to Angular position from the ascending node of an orbit to the spacecraft; 
argument of perifocus of an orbit 
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HILTOP is a Phase A performance analysis computer program that has 
been developed explicitly for the purpose of generating optimum electric pro- 
pulsion trajectories that are of interest for the exploration of the solar system. 
The program contains a propulsion system model that assumes jet exhaust speed 
(specific impulse) is held constant throughout the flight, while power and thrust 
may either be held constant, to simulate nuclear electric, or varied with solar 
distance, to simulate solar electric propulsion. The thrust direction time history 
and the engine on-off times are determined by the program to extremize a pre- 
selected performance index. The indirect optimization method is employed for 
this purpose. A basic knowledge of optimization theory is required to understand 
this report. For a discussion of optimization theory, see Reference [l] . 





The electric propulsion system is permitted to operate throughout the 
heliocentric flight. The initial conditions of the heliocentric flight are assumed 
to be established impulsively by a conventional launch vehicle. Two distinct 
problem formulations are available. In the launch vehicle dependent formulation, 
the initial spacecraft mass (i.e., the payload of the launch vehicle) is a specified 
function of the launch hyperbolic excess speed and launch parking orbit inclination, 
and these two parameters are available as independent parameters to optimize the 
distribution of performance between the launch vehicle and the electric propulsion 
system. In the launch vehicle independent formulation, the reference power is 
specified and the initial spacecraft mass is optimized for a given launch excess 
speed to maximize net spacecraft mass. In either formulation the hyperbolic 
geocentric trajectory ultimately established by the launch vehicle is then con- 
sidered from the viewpoint of asymptotic expansion theory and the spacecraft's 
motion is added to the motion of the Earth relative to the Sun. Only zero order 
terms of the expansion are retained such that the initial heliocentric position of 
the spacecraft is equal to the position of the Earth on the prescribed launch date, 
and the initial heliocentric velocity of the spacecraft is equal to the vector sum 
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of the Earth's heliocentric velocity and the launch hyperbolic excess velocity. 

The direction of the launch excess velocity is determined as part of the optimi- 
zation problem. 

HILTOP may simulate missions having as many as four targets (celestial 
objects) along a trajectory. One target is designated as the "primary" target 
(which may be absent in such missions as, for example, extra-ecliptic probes). 

In addition, up to three "intermediate" targets are permissible, such that the 
intermediate targets are visited prior to the primary target. At present, inter- 
mediate targets must be relatively massless objects such as comets and asteroids, 
because the gravitational perturbing effect of an intermediate target on the space- 
craft is not coded into the program. In addition, the spacecraft is optionally per- 
mitted to swing-by the primary target and continue onward ballistically to up to 5 
’^>ost-swingby" targets. A multiple -target mission may consist of any combina- 
tion of flybys, rendezvous, and stopovers involving the intermediate targets in- 
cluding possible sample retrievals and dropoffs of instrument packages. 

A power degradation model has been incorporated into the formulation of 
HILTOP. The model allows a single parameter (denoted "characteristic degrada- 
tion time") to describe the power degradation behavior of an electric propulsion 
spacecraft to a degree which fundamentally affects the solution to the trajectory 
optimization problem. 

The option of simulating spacecraft housekeeping power applies to solar 
electric propulsion with specified reference power. The housekeeping power 
is a specified constant power generated by the solar arrays and shunted away 
from the thruster power-conditioners and directly to the spacecraft payload for 
"housekeeping" purposes. A solution to the problem of optimizing electric 
propulsion heliocentric trajectories, including the effects of geocentric launch 
asymptote declination on launch vehicle performance capability, is developed 
using variational calculus techniques. The model of the launch vehicle per- 


formance includes a penalty associated with a non-easterly launch plus another 

* penalty arising from a non-coplanar launch from the parking orbit. Provisions 
for range safety constraints are included. Optimal trajectories will generally 
have the launch excess velocity offset from the initial primer vector. 

* 

The program’s ability to simulate all -ballistic missions includes powered 
and unpowered multiple swingby missions with an optional deep space burn. This 
capability renders HILTOP a powerful tool for ballistic mission design and optimi- 
zation, with tremendous flexibility for creating imaginative multi-target mission 
profiles. 

For planetary orbiter or asteroid or comet rendezvous missions, a com- 
putational technique similar to that for launch is provided to include a high-thrust 
retro maneuver to achieve the desired arrival conditions at the primary target. 

The arrival hyperbolic excess speed at the primary target is available as an in- 
dependent parameter to optimize the distribution of performance between the retro 
stage, of specified thrust and specific impulse, and the electric propulsion system. 

The use of the zero order asymptotic matching to account for high-thrust 
maneuvers in the vicinity of a planet permits all trajectory computations to be 
carried out in the heliocentric reference frame where a central, inverse square 
gravitational field is assumed. The location and motion of the planets within this 
system are defined as functions of time through an analytic ephemeris routine 
which contains osculating planetary orbital elements which are quadratic functions 
of time. Comets and asteroids are assumed to have constant orbital elements. 

The motion of the spacecraft in the heliocentric reference frame under the influence 
of the electric propulsion system is governed by a set of first and second order 
ordinary differential equations, known as the state equations. Associated with the 
state equations is another set of first and second order ordinary differential equa- 
tions, known as the adjoint equations, which are inherent in the application of the 
indirect method of optimization. These two sets of equations must be solved 
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simultaneously. During thrust phases the solution is obtained by numerical inte- 
gration; however, during coast phases, an analytic solution of the differential 
equations is known and is employed. 

The program is coded to yield a complete optimum solution by the indirect 
method only for the problem of maximizing net spacecraft mass. The reason for 
the restriction to a single performance index for the complete solution is that the 
transversality conditions, which comprise a portion of the Necessary Conditions 
of the solution, are dependent upon the choice of performance index and only one 
set of these conditions is presently coded. Coincidentally, however, a complete 
set of Necessary Conditions is also available for the problem of minimum flight 
time with fixed net spacecraft mass because the transversality conditions for this 
problem are identical to those for maximum net spacecraft mass with fixed flight 
time. A large variety of specific problems may be posed within the framework 
of these two general problems. 

The program iterator routine, whose primary function is to solve the 
boundary value problem that arises in all trajectory optimization problems, is, 
in fact, a generalized parameter optimization package. The iterator has two 
basic modes of operation: (1) the satisfy mode in which the sole purpose is to 
satisfy all specified boundary conditions, and (2) the improve mode in which a 
gradient technique is employed to improve a selected performance index while 
maintaining satisfaction of the specified boundary conditions. This iterator may 
be used in the following manner. If the specific problem to be solved is contained 
in one of the two general classes of problems indicated in the preceding paragraph, 
then all Necessary Conditions, including transversality conditions, are available 
for evaluation in the program. These transversality conditions are included as 
part of the boundary conditions, and there results a two point boundary value 
problem with an equal number of boundary conditions and unknown parameters. 

A solution of this boundary value problem is assumed to be a solution of the 
original optimization problem (specified problem to be solved) ; consequently. 


the complete solution to the original problem is assumed to be obtained upon the 
successful convergence in the satisfy mode. If, however, all of the required 
transversality conditions are not available for evaluation within the program, as 
in the case of a performance index other than as indicated in the preceding para- 
graph, then there will exist more independent parameters than boundary conditions 
that can be invoked. Thus, an infinite number of solutions exists to the boundary 
value problem posed, and all degrees of freedom may be used to improve the 
selected performance index. For such a problem, the iterator first solves the 
overdetermined boundary value problem in the satisfy mode, and then proceeds 
to the improve mode, in which the performance index is successively improved 
while maintaining approximate satisfaction of the boundary conditions. When 
the performance index can no longer be improved on successive iterations, the 
iterator assumes that the extremum has been successfully isolated and terminates 
the search. An important feature of the improve mode is that any program para- 
meter that is available for specification as a boundary condition may be chosen as 
a performance index. Consequently, the program provides considerable flexibility 
in treating a variety of problems. The primary drawback of the improve mode is 
that it is costly and time consuming to use extensively because of the normally 
slow convergence rates concomitant with direct techniques , and therefore it is 
highly recommended that the problem of interest be posed within the framework 
of the satisfy mode whenever possible. 
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n. 


FORMULATION 


A. 


SPACECRAFT AND TRAJECTORY MODELS. 


The following discussion is oriented toward the programming logic aspects 
of the HILTOP computer program. For the sake of simplicity and clarity, the dis- 
cussion of solar array radiation degradation is discussed separately in Section C. 


1. Spacecraft Mass Components. The spacecraft is composed of 
an electric propulsion system and associated tankage and propellant masses, a 
structure mass component, a retro propulsion component (for maneuvers about a 
primary target), a set of instrument package masses to be dropped at intermediate 
targets and a net spacecraft mass as follows: 


n-1 

m=m +m+m+m+m+)m , + m , 

o ps p t s r L drop i net 

i=l 


( 1 ) 


where m is the initial spacecraft mass; m , m , m. are the electric pro- 
o ps p t 

pulsion engine and powerplant, propellant and tankage masses, respectively; m 

s 

is the structure component; m is the retro propulsion mass; m , . is the in- 

r drop i 

strument package mass left at the i“* target; and m is net spacecraft mass 
(payload). In the analysis to follow, the subscript o denotes the launch body and 
n the primary (final) target. The net spacecraft mass consists of the scientific 
instruments, communications, navigation, and other engineering hardware, shield- 
ing, and any other mass components required to carry out the mission of interest. 


Under one program operating mode, known as the launch vehicle dependent 
mode, the initial spacecraft mass is governed by the performance capability of a 
launch vehicle according to the relation. 


where 






m . [ - t A ] . 

m =b e - b , 

o 1 3’ 


r 2 2T 
= v + v , 
L ®o e J 


( 2 ) 

(3) 
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provided that the declination of the departure asymptote at Earth escape does 
not exceed the latitude of the launch site (in magnitude). Whenever the departure 
asymptote declination exceeds the launch site latitude in magnitude, the above ex- 
pression for the characteristic velocity, v , should be expanded to allow non- 

c 

coplanar bums by the launch vehicle and/or non-due-East launches from the launch 

site. The equations for this more complicated case of a non-coplanar launch are 

treated in a subsequent paragraph. In expressions (2) and (3), v is the departure 

® o 

hyperbolic excess speed, v^ is an internal constant equal to the escape speed at 
185 km altitude above the Earth's surface, and b , b , b are constants which 

1 ' Z u 

define the performance capability of a given launch vehicle. These constants may 
either be input directly or taken from an internal table of constants that have been 
pre-computed for selected launch vehicles. The constants in the internal table 
were generated using a least squares algorithm and data from Reference [2] . The 
payload capabilities of the specified launch vehicles were curve fit to the above 
equation for initial spacecraft mass. 

Under a second operating mode, known as the launch vehicle independent 
mode, the initial spacecraft mass is completely unconstrained. The program then 
determines the optimum initial mass which results in the maximum net spacecraft 
mass for given reference power. The launch excess speed may be either fixed or 
optimized. 

The electric propulsion engine mass is given by, 

m ps = p ref [ “t + + Ap) “a 3 ’ (4) 

where p . is the reference power (see Electric Propulsion System) ol is the 
ref t ' 

specific mass of the thruster and power conditioning subsystem, a is the specific 
mass of the arrays, and Ap is the ratio of housekeeping to reference power, an 
input constant. The electric propulsion propellant mass m is obtained by inte- 
grating the derivative of mass ratio over all thrusting arcs and employing the 
equation. 


( 5 ) 


n-1 

m = m (1 - iM + Y* (m 
pn o n (-• samp 
i=l 


-m 


drop i\ 


where v is the mass ratio at the primary target (i.e. , prior to the optional 

retro maneuver) and m . is the sample mass picked up at the i^ target. 

samp i 


The sample masses and drop masses are specified as linear functions of 
the initial mass, as follows: 


m . = m k m . = m k (6) 

samp i o samp i drop i o drop i 


where k . and k , . are inputs and are available as independent para- 

samp i drop i 

meters of the boundary value problem. Both m . and m , . are available 

samp i drop i 

as dependent parameters. 


The electric propulsion propellant tankage mass and structure mass are 
computed, respectively. 


= k + ™ » 
t t p 


( 7 ) 


m = k m , 
s so 


( 8 ) 


where k, and k are input constants, 
t s 

The total retro mass array is given by the sum m = m + m , where 

r rp rst 

the retro propellant requirement, m , is given by the following: 

rp 


m = (m V - j m - i in ) e , 
rp on ps ps t t x 


( 9 ) 


where m and m, are the electric propulsion system and tankage masses, 
ps t 

respectively; j and j are input jettison indicators set equal to one if the 

PS L 

electric propulsion engine mass and tankage mass components are to be jettisoned 
prior to the retro maneuver and equal to zero otherwise; and e is given by 

■ IS 
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( 10 ) 


[-Av’/c ] 

- 1 r 

e — 1 - e 
x 

where Av' is the magnitude of the required velocity increment of the retro stage 

and c is the retro jet exhaust speed, evaluated as the product of the input retro 

r 2 
specific impulse and the reference acceleration of gravity (9.80665 m/sec ). The 

computation of Av' is discussed in Section B, Part 3, Target Conditions and 

Spacecraft Constraints. 

The retro structure and tankage mass component, m , , is given by 

rst 


m 


rst 


= m 


+ k , m 
rs rt rp 


( 11 ) 


where m is the input structural mass and k 
rs rt 

The bum time associated with the retro maneuver, 
of the input engine parameters. 


is the input retro tankage factor. 

t , is computed as a function 
b 


m c 
= , , r p r 

f 

r 


( 12 ) 


where f 
r 


is the input thrust level of the retro stage. 


2. Electric Propulsion System. The independent parameters 
associated with the propulsion system are reference thrust acceleration, g, 
and jet exhaust speed, c. g is the thrust magnitude at 1 AU from the sun 
divided by the initial mass, and is not to be confused with the acceleration of 
gravity at the Earth's surface, which is not used (as a symbol) in this document. 
The instantaneous thrust acceleration is computed by. 



( 13 ) 


where y is the instantaneous power ratio (described below), v is the instan- 
taneous mass ratio m/m Q , and h^ is a step function equal to one if the engine 
is operating and equal to zero otherwise. The propulsion system is assumed to 


operate at constant jet exhaust speed. The instantaneous power, p, delivered to 
the power processors, is computed by 

p = gm o cyh CT /277 = h CT yp ref , (14) 

where T) is the engine efficiency and the reference power, p is the power 
delivered to the power processors at 1 AU from the sun. 


The thrust magnitude, f, generated by the electric propulsion system is 
given by 

f=gm o yh (r =2ffl?. (15) 

The propulsion system efficiency is assumed to be a function of the jet ex- 
haust speed and is given by 

be 2 

V= -f 2 + e * a6) 

c + d 

where b, d, and e are constants specified in program input. 


The power ratio y is defined 


y= 



(17) 


The electrical power model of the electric propulsion system is specified through 
an input option indicator. For nuclear electric spacecraft, the power remains 
constant throughout the trajectory. For solar electric spacecraft, the power is 
modelled as a function of solar distance and solar array orientation relative to 
the sunline. The six options available are given below: 


Option 1 Not presently used. 

4 

Z i/4 

a. d - Ap , (18) 

i=0 
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Option 3 
Option 4 

Option 5 
Option 6 


where d represents the density of photons impinging 
the solar arrays. The definition of d is 


d = 


cos X 

2 

r 


( 19 ) 


where X is the angle between the normal to the array 
and the spacecraft-sun line. For Option 2, X = 0. 

The coefficients a^ comprise a set of internal constants 
with the following values: 

i a. 

i 

0 0.6270 

1 5.3054 

2 -10.0376 

3 7.1073 

4 -2.0021 

These values are consistent with solar distance r expressed 
in AU. Also, it is possible to input the coefficients a i to 
the program. (See HILTOP Input). The program expects 
the y-curve to peak at some distance r ^ from the sun, 

and r , is about 0.665 AU for the coefficients listed 

peak 

above. 

7 = 1 . 


Same as Option 2 for r greater than r ^ For r less 

than r . , y is held constant at the peak value by setting 
peaic 2 

“'X'I'V • 


Same as Option 2 with the side condition that y s GAMMAX 

where GAMMAX is an input constant. This is achieved by 

setting cos X “ (r/r ) 2 where r is the distance at which 
* c c 


y = GAMMAX. 


Same as Option 2 except y is held constant at the peak value 
everywhere. 
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These five available options may be described roughly as follows. Option 2 is the 
« current projected-state-of-the-art basic SEP power law. Option 3 is the nuclear 

electric propulsion (NEP) power law. Options 4, 5, and 6 are variations of 
Option 2, in which the solar panels are either shielded or tilted away from the 

* 

sun at the smaller solar distances, and in Option 6 it is assumed that reflecting 
flaps gather the required solar energy at the greater solar distances. 

3. Target Ephemerides. When it is desired to generate trajectories 
involving targets which are specific celestial objects, the program includes the 
capability of internally generating a specific target’s state (position and velocity) at 
any given time. If a target is a major celestial object in terms of mass (Mercury 
through Pluto), the orbital elements eccentricity, argument of perihelion, node, in- 
clination and mean anomaly are computed as quadratic power series in time, and 
the semi-major axis remains constant with time. If a target is a minor celestial 
object in terms of mass (comet or asteroid), the orbital elements are constant 
with time. In either case, Kepler's equation M = E - e sin E is solved iteratively 
for the eccentric anomaly E; in this equation, M is the mean anomaly and e is 
the eccentricity. Finally, the target's position, velocity and acceleration are ex- 
pressed in a Cartesian coordinate system with the x-axis directed toward the 
vernal equinox of date, the z-axis pointing toward the north ecliptic pole, and the 
y-axis completing the right-handed system. 

4. Differential Equations. The attainment of an optimum electric 
propulsion trajectory requires the repeated simultaneous solution of two sets of 
differential equations, the state equations and the adjoint, or Euler-Lagrange, 
equations. The state equations are comprised of the equations of motion of the 
spacecraft plus the equations describing the change of the spacecraft as a function 
of time. 

The motion of the spacecraft is assumed to take place in an inverse square 
heliocentric gravitational field and to be influenced by an electric propulsion system 
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with thrust directed optimally. The second order differential equations of motion 
are given by 


• • 

R = h 

a 


v 



( 20 ) 


where R is the heliocentric position of the spacecraft; r = |r| ; is the sun's 
gravitational constant; g is the reference thrust acceleration (thrust magnitude 
at 1 AU from the sun divided by initial mass); y is the normalized power varia- 
tion which may be a function of solar distance and solar array orientation; v is the 

mass ratio m/m ; e A is a unit vector in the direction of thrust; and h is the step 
o t Q 

function control variable associated with the thrust switch function. Only the sun's 
gravitational effect on the spacecraft motion is taken into account, the effects of all 
other celestial objects being neglected during heliocentric flight. The complete set 
of state equations for this problem includes the equations of motion above and the 
following three first order differential equations: 


v = -h 


gy 


cr c 


( 21 ) 


g =0, 

(22) 

c = 0 , 

(23) 


where g and c are the propulsion system parameters, defined previously, which 
are constant on a given trajectory. In addition, the following differential equations 

T = h^ , (24) 

0=0, (25) 

where T is the propulsion time and 0 is the thrust cone angle (the angle between 
the radius and thrust vectors), are considered state equations for problems with 
constrained propulsion time and fixed cone angle, respectively. The condition of 
fixed thrust cone angle is implemented through the latter differential equation in 
conjunction with the side condition. 


e • R - r cos 0 = 0. 

If 


( 26 ) 


The differential equations which govern the behavior of the adjoint variables 
are given by 


A=2£(A-r) R -M + V- 

5 ' 3 V 

r r 


y' - (A • e - - X ) 

T ^ t C V } 


— R 

+ X (e., - — cos 0 ), 
x t r 


(27) 


h g y 

V - r" (A ' v • 


(28) 


« 

X 

g 


- (A * e - - X ), 

v ' t c i/' ’ 


X 

c 


h a gy 

2 

c 


X 




(29) 


(30) 


where y' = dy/dr; A is the familiar primer vector which is adjoint to the 
velocity; X^ is the adjoint variable associated with the mass ratio; X^ is the 
adjoint variable associated with the reference thrust acceleration; X c is the ad- 
joint variable associated with the jet exhaust speed; and X x is a Lagrange 
multiplier that is identically zero if thrust cone angle is unconstrained and de- 
fined 


X 

x 




(mxe t ) 

— 9 
(mxe t ) 


(31) 


if the cone angle is fixed. The unit vector m is defined 


R x A 

m=- 


R x A 


(32) 


Also, if thrust cone angle is fixed, an additional adjoint variable X^ is introduced 
which is adjoint to the cone angle 0 , and its differential equation is 
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Finally, a constraint on propulsion time requires the introduction of X^., an 
adjoint variable associated with propulsion time. Because propulsion time does 
not appear explicitly in the state equations, \ is a constant, i.e. , 

X T = 0. (34) 

If propulsion time is optimized, the constant value of X is zero. 


5. Optimality Conditions. The control variables available for 
optimization along the trajectory include the unit thrust direction vector e and 

V 

the step function h^„ The application of the Maximum Principle of optimal control 
theory leads to the result that the proper choice of the control variables is that 
which maximizes the variational Hamiltonian, h^, at each point along the path. 
The variational Hamiltonian for the problem formulated here may be written 


='V J tr < A • ® t - * V 77 V " i (A ‘ R) - A ' * • 


(35) 


Since the quantities gyfv and h are non- negative, it is obvious that h 

O v 

is maximized with respect to e by aligning e with A, and this is the choice 

l I# 

made if the thrust cone angle is unconstrained. In the event of a constraint on cone 
angle, one simply chooses e^ as close to A as the constraint will permit. By 
considering the intersection closest to A of a circular cone about R (with half- 
angle equal to 0) and the R - A plane, one easily determines the optimum e^, 
subject to the fixed cone angle constraint, to be 



K 1 


Rcos 0+ (mx R) sin 



(36) 


with m as defined in expression (32). 


The proper choice of h^ 
a, defined by 


is seen to depend on the sign of the switch function 



( 37 ) 


. V . V . 

ff = A * e - - A + — — A . 

t c v g y r 


Again, since gy/v is non-negative, the variational Hamiltonian is maximized 
with respect to h by choosing h^ as follows: 

h =0 if a < 0 , (38a) 

Or 

h = 1 if 0- > 0 . (38b) 

Recall that the permissible values of h^ are limited to zero and one. 

6. Integration. Integration of the differential equations associated 
with the state and adjoint variables is done in two distinct modes which depend upon 
the thrust switch function, o. These two modes are discussed below. 


Thrust 


Numerical integration is required on thrust intervals. The independent 
variable of integration is the generalized universal anomaly, u. The relationship 
between derivatives with respect to time and u is given by 


du -n 
dt r 


(39) 


where r is the spacecraft's solar distance and n is an input constant (not to be 
confused with subscript n used elsewhere). Au, the independent variable interval 
used for the integration, is an input constant. Denoting derivatives with respect to 
u with the prime, the conversions from time to u derivatives are 


f “ 

x - r x , 


( 40 ) 


.. 2n ( R • R 
x" = r ^ — - — nx 


K+x). 


( 41 ) 
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A standard fourth-order Runge-Kutta numerical integration technique for 
first-order differential equations is used. A value of n equal to 1.5 will 
regularize the differential equations with respect to the mathematical singularity 
associated with the gravitational force at the point r = 0 (where the sun is located). 
This is accomplished by inherently taking smaller steps in time (which is an inte- 
grated quantity) when closer to the sun and larger steps when farther from the sun, 
assuming Au = constant. Equation (39), however, allows the removal of only the 
solar singularity from the optimal rocket problem, and has no provision for the 
other singularity, which is characterized by an extremely high thrust rotation rate 
and occurs whenever the primer vector passes relatively close to the origin of 
primer-space (A = 0) during a thrust phase. The difficulty associated with the 
primer-origin singularity is lessened by continuously cutting down the step size 
Au as the primer origin is approached. It is noted in passing that when forced- 
thrusting is invoked for an optimum flyby mission, the primer-origin singularity 
will occur at the flyby target (as will become apparent in the discussion of trans- 
versality conditions in a later section). 

Coast 

During any coast phase, the two-body equations of motion and the associated 
adjoint equations are known to possess analytic solutions obtainable in closed form. 
The particular form of the solution used in HILTOP is derived in Reference [3] 
and is simply repeated here. This solution employs a universal variable, 0, de- 
fined implicitly through the differential equation 

(42) 

where fi is the gravitational constant of the attracting body. 


For elliptic, two-body trajectories, this equation has the solution 


where E is the eccentric anomaly and a is the semi-major axis. 

The problem to be solved is that of evaluating the state and adjoint variables 
at + A)3 given the values of these variables at /J= /S . This is accomplished 

o 

as follows. Given R and R (the state at fi = p ) compute 

o o o 

i 

r = (R • R ) 2 , 
o o o 

2 

v = R • R , 
o o o 


d = R • R , 

O O 0 


(44) 


1 = 2_ 
a r 


o 


e= (A£) 2 /a = (E-E q ) 2 . 


Then, using the truncated infinite series expression 

16 k 

V^X -0m ‘ 


k=0 


(45) 


compute the functions G. for i = 5 and 4. Thereafter, the functions G., 

i i 

i = 0, 1, 2, 3 are computed from the recursion formula, 


, = (A/?/ _1 r 

r i i! a i+2 


(46) 


and are employed to evaluate the familiar f and g functions, i,e„. 


f = 1 - 


e = (r G + — — G ) , 

° 1 2 

d 


r = r G + G., +G 0 , 

° 0 Vm 1 2 


(47) 
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t-t =g + - 2 - , 
° 

& G, 


f = - 


r r 


o 

• , 

g=l-— , 


R=fR +gR , 
o & o 


R = f R +g R , 
o o 


(47) 

cont. 


which provide the state and time at the given value of (3= f$ + A/3„ g in this 

o 

section is not to be confused with the reference thrust acceleration used widely 
throughout this document. The corresponding equations for the adjoint variables 
are: 


dx. (t) dx. (t) 

X.(t) = r — \.(t )+ X.(t ), 

1 3x. (t ) i o' . l o' 

j o' J S x . (t ) J 

J° 

(t) ax.(t) . 

jo' 


(48) 


(49) 


where i, j = 1, 2, 3, and repeated subscripts in the same term imply summation 
over the range of the subscripts. The variables x.(t) represent the three 
Cartesian components of R(t) while the X.(t) represent the components of A(t). 
The partial derivatives indicated are given as follows, with 6 denoting the 
Kronecker delta function: 


dx. 


o. 

J 


P l X o. 


2G - AjSG, 


• • r i r • • 1 v \ r/ \ 

v i o.'L 3 u. J o/J ij 3 L\ 2 r / o 

i r ^ J r o 

A " A 


+ < 3G 5 - A ^G 4 , -f ] ■ 
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(50) 



-v X . -X 

ox. i o. 


dx 


J- 1 frl* o, - Vo, > 7T 1 [«»4 - 77 + < 3G 5 ' A ^ G 4> ^ ]+ • 


o. 

J 


J J 


(51) 


^ 0 n P- x 

ox . PX. _ r x o. 


I 


x . - X p.x X 

i o. _ 2 o. > o. 


dx 3 

o r r 

J o 


• * v. « V. i U . U. v _ 

-T — - J - + -(x.-x )> Lr_L /_^ G + __1 G )1 

L ^ 3 ix K j o/J r L „ 3 \ r Q o j- l/J 


r r 


K°1 


+ 


G 3 -A0G 2 


x + (2G - A0G ) x 1 + f 6 .. , 
/ o. 4 3 o.J 11 

i i * 


(52) 


dx. x X i" X o. - P 2 X o. 

— = — ~ [p x - G x 1 + — f - G x 1 

3x q r 3 Ll °j 2 V r L SH 1 °j J 

j 

x _ 

+ £ [ <G 3 ' A *V # *0, + <? G 4 - A *V * o ] ■ * i 6 iJ • 


(53) 


o i 


where 


Pl - TT [ 3G 5- A/JG 4 + J- < 2G 4' A ^ G 3 )+ r o (G 3- A *y] • 


(54) 


P 2 =2G 4 -A£G 3+ ^ (G 3 -A^G 2 )-r o A^G 1 . 


(55) 


The values of the derivatives of all other state and adjoint variables vanish 
during coast phases ; therefore, their solutions are trivial. 


7. Units. The internal units of the program are kilogram for mass, 
AU for distance and tau for time, which implies speed in EMOS. (One tau is the 
time required for a massless particle to travel one radian in a circular orbit of 
radius 1 AU about the sun). To convert from internal units to the MKS system, 
multiply by the following constants : 
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Distance - r = 1.49599 x meters = one A’J 

,2 -3,2 

Acceleration - jU/r = 5.9301282604 x 10 m/sec 

Velocity - *J\L / r = 29784. 916613 m/sec = one EMOS = one AU per tau 

Time - J t 3 /(jl/86,4C0 = 58.132440991 days = one tau 

20 3 2 

A solar gravitational constant of 1.32715445 x 10 m /sec is assumed. There 
is no conversion from input to internal units for the initial Lagrange multipliers, 
which are required to start a trajectory. Hence the input values of the multipliers 
are consistent with the internal units of the state variables; i.e. , mass in kilograms, 
distance in AU, and time in tau. 

When evaluated in internal units, the equations presented previously yield 
power in units of kg - AU /tau and thrust in units of kg - AU/tau . To obtain 
these quantities in the MKS units (i.e., power in watts, thrust in newtons), multiply 
by the following factors: 

— q — rr~ o o 

fj, /r = 176.6283757392 m /sec 
Thrust - h/t 2, = 5. 9301282604 x 10 2 m/sec 2 

Thrust in pounds is obtained by dividing the thrust in newtons by the factor 
4.4482221811. 
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B. 


BOUNDARY VALUE PROBLEM 


1. Boundary Constraints. The HILTOP program currently provides 
for nearly 100 different dependent parameters (end conditions), although certain of 
these are coded mutually exclusive. There are about half as many independent para- 
meters available in the program and therefore a subset of the dependent parameter 
array must be chosen for a well-posed problem. The independent parameters include 
the adjoint variables (i.e. , the primer vector and its time derivative at the start of 
each possible trajectory segment and the adjoint variables associated with mass ratio, 
degradation time and propulsion time at the start of the first trajectory segment), 
the reference thrust acceleration, the jet exhaust speed, the launch asymptote de- 
clination, the launch time and hyperbolic excess speed and the times and excess 
speeds at the various possible targets. Also included are the thrust cone angle (if 
constrained), the inclination of the parking orbit about Earth, the heliocentric de- 
parture velocity, and the sample-mass and drop-mass factors at each possible 
intermediate target. The dependent parameters include the spacecraft position and 
velocity differences with respect to each possible target along a trajectory, the 
launch time and hyperbolic excess speed and the times and excess speeds at the 
various possible targets, and, indeed, all of the independent variables (except the 
primer vector time derivative) are available also as dependent variables, as this 
allows the analyst to let the program iterator move one trajectory solution toward 
another trajectory solution. This provides wider flexibility and more control in 
solving the boundary value problem, particularly in instances where successive cases 
are employed to sweep a range of values of one or more of these parameters. In 
addition, total flight time, propulsion time, degradation time, net spacecraft mass, 
reference power, final solar distance, and the transversality conditions which yield 
optimum flybys, launch time and excess speed, target encounter times, reference 
power, reference thrust acceleration, jet exhaust speed, launch parking orbit in- 
clination, launch asymptote declination, thrust cone angle, and travel angle are 
available as dependent conditions to be satisfied. The complete list of options avail- 
able is found in the section Definitions of Input Parameters. 
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2. Initial Conditions for State and Adjoint Equations. The launch 

o 

planet's position P q and velocity P q are either input directly or computed for 
the specified launch date using the ephemeris routine that is available in the pro- 
gram,, Given these vector quantities, the program computes the spacecraft's 

initial heliocentric position R to coincide identically with the launch planet, 

o 

= P^o The spacecraft's initial heliocentric velocity may be input directly to 
the program; however, the program's normal operating mode is to equate the 
initial velocity to the sum of the velocity of the launch planet and the hyperbolic ex- 
cess velocity, i.e., 



+ V 


COQ 


(56) 


where V is the launch hyperbolic excess velocity. The magnitude v bo is 
. available as an independent parameter, and therefore may be either specified or 
optimized. The direction of is usually optimized , although one option does 
permit partially constraining through the geocentric declination which also is 
available as an independent parameter. The initial mass ratio V is specified to 
be one on all trajectories. 


The adjoint variables at the start of each possible trajectory segment, which- 
are contained in the list of independent parameters, are direct inputs to the program 
for the first trajectory of an iteration sequence, except that some of the adjoint 
variables associated with departing from an intermediate target may be generated 
internally by the program as being continuations from the previous trajectory seg- 
ment. For all subsequent trajectories in an iteration sequence, these adjoint vari- 
ables either remain as input or are varied by the iterator. A value for the initial 
adjoint variable associated with mass ratio must be input to the program, and a 
value for the initial adjoint variable associated with propulsion time may be input 
(the default-value of zero yields optimum propulsion time). The values of the 
initial adjoint variables not included in the list of independent parameters (i.e. , 

Ag, A c> and A^) are set to zero because the specific initial values assigned to 
these variables are arbitrary. 
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3. Target Conditions and Spacecraft Constraints. Three basic 
options are available for specifying the target conditions with each option being 
designed for a particular type of mission. Within each option, considerable 
flexibility is available for treating many variations of a mission type. The three 
options will be referred to as the ephemeris, the open angle, and the extra-ecliptic 
options. 

As the name implies, the ephemeris option is employed for missions in 
which the primary and/or intermediate targets are celestial bodies for which the 
positions and velocities (of the bodies) are either input directly or computed using 
the analytic ephemeris. Denoting a specific target’s position and velocity as P. 
and P., respectively, a constraint on the spacecraft position R^ is imposed by 
nulling the position error; i.e. , by forcing the satisfaction of the equation 

AR. = R. - P. = 0. (57) 

ill 

Similarly, a constraint on the spacecraft velocity R. is imposed by nulling the 
velocity error, 

AR.=R. - P. - V . = 0 , (58) 

i i i “i 

where V . is the arrival excess velocity at the i^ target. These equations 
00 1 

apply both to intermediate and the primary targets. 

The open angle option is restricted to problems of two-dimensional motion 
in the x-y plane. The option is designed specifically for the problem of open angle 
transfer from a given point to a specified solar distance r^. This target condition 
is written, simply, 

|R n l=r f , (59) 

where the subscript n denotes the primary, or final, target. The capability of 
imposing circular orbit conditions at this solar distance is also available. This 
vector target condition is written in the form of a velocity error as follows: 
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( 60 ) 


R 

AR = R - ,/ fkx — S- ) - V =0, 

n V \ j | / “n 

1 n 

where fi is the gravitational constant of the sun and k is a unit vector along the 
z-axis. 

The extra-ecliptic option provides the capability of targeting to a final 
perihelion distance r^ , inclination to the ecliptic i^ , and orbital eccentricity 
e f starting from a launch planet of specified position and velocity. The formula- 
tion imposes the condition that the spacecraft be at perihelion at the final time. 
Thus, four of the six conditions required to completely define the final position 
and velocity are specified. The two open degrees of freedom lead to a like number 
of transversality conditions which are 

H • C = 0, (61) 

k • C = 0 , (62) 

where H is the angular momentum vector of the final orbit, k is the unit vector 
normal to the ecliptic, and C is the vector constant of the motion defined by 


• . 

C = (Rx A) - (Rx A). (63) 

These equations may be solved for the osculating orbital elements, £2, the longi- 
tude of ascending node and, U), the angular position from the ascending node, 


evaluated at the final time. This leads to the relations 



O =tan _:L [{C • j )/(C • i )] , 

(64) 


co=tan 1 ^(-h* A)r f /(h» 

(65) 

where 

h = cos i^k + slnij (kxC)/|kxC | , 

(66) 


v f = JnQ.+e f )/r f , 

(67) 


26 


and where A and A are evaluated at the final time. The target conditions are 
then written as position and velocity errors as follows: 


c to c Ci - s to s ft c i. 


AR = r„ 


c to s £2 + s to cCi c L 


- R = 0 , 


( 68 ) 


s to S 1, 


AR = v 


-s W cfl - c to s c i^ 
-stosO+ctocficij 
c to s i^ 


- R = 0 , 


(69) 


where s and c denote sine and cosine, respectively, and R and R are the 
final integrated heliocentric position and velocity vectors, respectively. 


An alternate set of end conditions for the extra ecliptic mission is also 
available. The alternate set relaxes the requirement that the final position be at 
perihelion. The final state is specified in terms of the semi-major axis, 
eccentricity, and inclination. The other three end conditions are transversality 
conditions associated with open longitude of node, argument of periapse and true 
anomaly. The first two of these conditions are given by (61) and (62), and the 
third is: 

£ (A • R) + r 2 (A • R ) = 0. (70) 

For either of the first two options above, the capability is provided for 

evaluating the propellant and structural mass requirements of a high-thrust 

retro stage which brakes the spacecraft, approaching the primary target along 

a hyperbolic orbit of excess speed v , into an elliptical capture orbit of peri- 

°° n 

center distance r and apocenter distance r . The retro stage thrust f and 
jet exhaust speed c^ are specified by input. The retro maneuver is assumed 
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to take place at the periapsis of the approach hyperbola; therefore, the impulsive 
change in velocity is given by 


where 


A , 2 , 2i f 2r a v p 

- ( v „_ + 2v n )-(, 


5 n 


(r +r ) 
a p' 



2 

v 

P 



> 


(71) 

(72) 


and n is the gravitational constant of the primary target. Provision for includ- 
ing the finite thrust velocity penalty is optionally available. This feature employs 
the theory developed by Robbins (Reference [4]). The total velocity required, in- 
cluding the velocity penalty, is given by 


Av' = Av + c„ e f , 
1 x x 


(73) 


where Av' is solved iteratively with the following equations. 


(v 2 + v 2 )/(v 2 + 2v 2 ), 
°°n p ' °°n p ' 

(74) 

/ v r, < m o v n ■ La " K \ 2 

r p r on psps tt> 

i 2rf /’ 

(75) 

p r 


(-Av'/c ) 

■ e' r , 

(76) 

2 c 


+ Av' ] ®x ’ 

(77) 


and where j and j are the jettison flags or indicators defined previously in 
ps t 

Section A. 

The usual definition of characteristic speed v , used in the computation 

c 

of initial spacecraft mass, is given by equation (3). As stated previously, this 
definition is valid only if the magnitude of the geocentric declination 6 of the 
launch excess velocity is less than or equal to the latitude of the launch site. If 
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this condition is violated, then the launch vehicle payload (i.e., the initial spacecraft 
mass) becomes a function of the direction of V also. High geocentric declinations 
of can be achieved by establishing a high inclination launch parking orbit, by 

employing a non-coplanar injection from the parking orbit, or by a combination of 
both. The usual method of defining launch vehicle performance corresponds 
to a due -East launch from the ETR. Such a launch yields a parking orbit in- 
clination equal to the launch site latitude, or about 28.5 degrees. The HILTOP 
model of initial spacecraft mass for high declinations consists of adding to the normal 

definition of v a velocity penalty associated with a non due-East launch, to achieve 
c 

a greater parking orbit inclination, and a penalty due to the non-coplanar injection 
maneuver. 

The velocity penalty incurred with non due-East launches from the ETR is 
shown graphically in Reference [2 ] as a function of the parking orbit inclination. 

This velocity penalty Av. is adequately approximated with a quadratic curve fit 
of the form 

AV i = C l i2 + C 2 i + C 3‘ C 78 ) 


Given a reference orbit inclination i and a circular orbit speed v , the 

o 

velocity penalty Av associated with a non-coplanar departure from this circular 

& 

orbit to the desired hyperbolic excess velocity at a declination 6 is defined as 

follows. Assuming the line of nodes of this reference orbit is an open variable, one 

may choose this variable to minimize the angle between the excess velocity and the 

[ 5 ] 

orbital plane. This minimum angle is 6 - i. Gunther has shown that the mini- 
mum incremental velocity required to achieve a given v along an asymptote not 
lying in the orbital plane from a specified circular orbit is obtained from the solution 
to a quartic equation in the sine of the out-of-plane angle. Defining 


s =sin(6-i); p =v /v ; 

®o o 

2 2 

p = s (p +4); 


( 79 ) 
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2 n 2 2 
q = s (1-s )p ; 


= [y(q/2) 2 + (p/3) 3 + q/2] l/3 - \J (q/2) 2 + (p/3) 3 - q^]^; 


1/3 


w 


y = yp 2 /4 - X ; 

= |[p/ 2 + y+ y(P/ 2 + y ) 2 + 4 (x/2+ y^A + s 2 ) ].’ 


(79) 

(cont) 


then Gunther’s solution for the magnitude of the minimum velocity impulse required 
to accomplish this maneuver is 


g 


“ v q Jp 2 + 3-2 y (1 + pw - w 2 )(2 + pw) , 


(80) 


and the penalty Av is the difference between v and the velocity increment re- 

g g 

quired if the out-of-plane angle were zero, i.e. , 

Av =v - (Jv 2 + 2v 2 - v ) . (81) 

g g W ®o o of 

Thus, the definition of the characteristic speed for those cases in which the 
asymptote declination lies outside the interval [-i, i] is 


7 2 2 

v +2v + Av. + Av = v +v + Av. . (82) 

°°o o i g o g i 


The transversality conditions to be developed later define the optimum 

choices of 6, i, and the geocentric right ascension & of the excess velocity. 

Once these are known it is necessary to evaluate V in the ecliptic Cartesian 

oo o 

coordinate system. This vector is evaluated 


V =v 

0OQ OOQ 



COS CL COS 6 
sinO! cos 6 


sin 6 


9 


(83) 
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where $ is the transformation matrix, 




10 0 
0 cos e -sin € , 

0 sin e cos c 


(84) 


and e is the obliquity of the ecliptic. A more complete development of this 
formulation, with examples, is given in Reference [6] . 


At intermediate targets, the mass ratio may be discontinuous to account 
for mass drops and sample pick-ups. Since these mass components are defined 
to be proportional to the initial mass, the mass ratio discontinuity is 


Av. = vt - vT = k . - k, 
ill samp i drop i 


(85) 


Position and velocity at intermediate targets are continuous. The intercept of an 
intermediate target is achieved by imposing a constraint on the position error. At 
the i** 1 target, this constraint may be written 

AR. = R.-P.= 0, (86) 

ill 

One may optionally constrain the passage speed v . at an intermediate target. 

The constraint equation for this is 

AR.=R. - P. - V . = 0. (87) 

i i i ®i 

The direction of V . is optimized. 

OOj 

In simulations of trajectories which are all-ballistic, the program is capable 
of simulating a single deep-space burn, or impulsive velocity-change, at any point 
prior to arrival at the primary target. The three components of the incremental 
velocity AV are independent variables of the boundary value problem, such that, 
at a specified time, the spacecraft velocity is incremented: 
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R = R + AV . 


( 88 ) 


The use of this program option is described in the Sample Problems and Results 
section under Sample Case H, Multiple Ballistic Swingby Mission. 

In the HILTOP program, stopover missions having optimum stopover time 
are simulated simply by forcing the spacecraft to rendezvous with the desired inter- 
mediate target. If the trajectory segment immediately following the intermediate- 
target arrival-time begins with a coast phase, then the duration of that coast phase 
is the optimum stopover time. If that trajectory segment begins with a thrust phase, 
then the optimum stopover time is zero. To simulate a stopover mission having a 
specified stopover time, as in Sample Case E displayed later in this document, 
the same intermediate target should be specified twice consecutively, and of course 
the spacecraft should be forced to rendezvous with the intermediate target at the 
first encounter. Then inputting values for A and A at the start of the stopover 
trajectory segment (as boundary value problem independent variables) to be re- 
latively small with respect to the mass ratio multiplier X ^ will force the thrust 
switch function to be negative and cause the spacecraft to coast along with the inter- 
mediate target until the desired departure time is encountered* . In this manner the 
trajectory block print and extrema of selected functions are available during the 
stopover phase. 

The capability is provided for constraining the total flight time of the mission 
(not including the flight time of the optional ballistic swingby-continuation trajectory 
segment). This feature is particularly useful in problems for which the analytic 
ephemeris is employed and for which the launch and (final) arrival dates are opti- 
mized but flight time is specified. The end condition is written 



( 89 ) 


where t, is the specified total flight time and t and t are the launch and (final) 
I on 

arrival dates, respectively. 

*Or, alternately, an imposed coast phase maybe employed. 
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The capability for constraining propulsion time is implemented through 
satisfaction of the optional end condition, 

T f"/ n ' h a dt (90) 

O 

where h is the step function defined in Section A and T. is the desired pro- 
a f 

pulsion time. It should be noted that this constraint is an equality constraint; that 
is, if invoked, the propulsion time is forced to be equal to the input value (after 
convergence is achieved) regardless of whether this input value is less than or 
greater than the optimum propulsion time. Given a solution for constrained pro- 
pulsion time, one can easily determine whether this propulsion time is less than or 

greater than the optimum value by the sign of the associated adjoint variable, X . 

T 

A negative value of X indicates the constrained propulsion time is less than 

T 

optimum whereas a positive value of X indicates the constrained propulsion time 

T 

is greater than optimum. 

Several spacecraft parameters are available as constraints. These are the 

net spacecraft mass, m , the intermediate target drop and sample masses, 

net 

m drop i anc * m sa mp i’ res P ec ti v ely, and the reference power, P re £» The equations 
employed in computing these two parameters are given in preceding sections, 
Spacecraft Mass Components and Electric Propulsion System, respectively. 

All trajectories are integrated forward from the launch date t and are 

terminated at the time (t -t ) later. Primary-target conditions are computed 

n o 

using the variables evaluated at this trajectory termination time. For multiple- 
target missions in which intermediate targets are present, the integration is 
interrupted at intermediate times t. and appropriate quantities are stored for 
the computation of dependent conditions. 

4. Transversality Conditions . The application of the indirect method 
of optimization leads to a set of necessary conditions, some of which are known as 
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transversality conditions, that must be satisfied by the solution. In essence, for 
every boundary condition left open* in the problem posed, the indirect method pro- 
vides a transversality condition. For a given performance index n which is to be 

minimized, the general equation for the transversality conditions is written 
n 


kd?r + 


It* 


dR - A • dR+A 


v 


dl/' + A dg + A dc+A d0 

g C 0 


+ A 


dr-h dt 
T v 




(91) 


The convenient choice is made whereby A , A , and A , are forced to be con- 

g c 0 

tinuous at each intermediate target, which means that, for example, only A (t^) 
need appear in the derived transversality expressions rather than the cumbersome 
expression 

n-1 


vv-l vv-vv'-yv- 

i=l 


This is because A (t ) alone, with A (t ) = 0 and A (t.) = A (t.) for each i, 
g n' g o' g i g i 

has the same value as the cumbersome expression cited above if A^ (t Q ) were not 

zero and A (t.) were not continuous, and this is due to the absence of A in the 
g i g 

differential equations, the same being true for and A^. The scalar k is an 
arbitrary positive constant which expresses the arbitrariness of the performance 
index; in other words, the minimization of n is equivalent to the minimization of 
2 IT, 3ir , . . . , etc. k effectively renders the general transversality condition linear 
and homogeneous in the adjoint variables, thus allowing the elimination of one ter- 
minal condition from the problem by appropriate choice of a value for k. Due to 
the independence of an unconstrained boundary condition, its differential is 
arbitrary, forcing its coefficient to vanish independently of all other terms in the 
equation and thereby yielding the desired transversality condition. 


In the program HILTOP, all available transversality conditions are derived 

for the problem of maximizing net spacecraft mass, i.e. , V— - m . From the 
1 net 

♦The terminology "open" is synonymous with unspecified and "fixed" is synonymous 

with specified. An open parameter is generally optimized. Also, final time refers 

to time at the primary target. 
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earlier definition of m , , one may write 

net J 

ff= j m +m jk +k -(l + k)l/ +j (1 + k Je |(l + i kjy 
J r rs ols t v r n J r v rr xL v J t V n 


n-1 n-1 

-iAa+Y (k. k 


t t i—i ' samp i drop i 
i=l 


L samp i 
i=l 


n-1 


-k V k .)+m fl-j j (1+k )e 1, 
t L dropiJ ps L J r J ps rt xJ 
i=l 


(92) 


where j is a constant equal to one if a retro stage is employed and equal to zero 

otherwise. For the launch vehicle dependent formulation, m^ is a function only 

of the launch excess speed v and possibly the geocentric declination 6 and 

launch parking orbit inclination i, and m is evaluated 

ps 


where 


m =0!p =agcm /277 , 
ps ref & o' 


a = a + (l + Ap) a . 

X cl 


(93) 


Of course, if reference power is constrained in this formulation, then m is 

ps 

fixed, and only two of the three parameters g, c, and m are independent. On 

o 

the other hand, for the launch vehicle independent formulation the reference power 
(and therefore m^ g ) is always fixed and m^ is obtained directly from the equation. 


m o = 2 ”P r ef /gc • 


(94) 


such that m Q is functionally dependent upon g and c. Thus ir may be written 
functionally in its most general form. 


Tr=ff < v m«» g, c, 6, i) . 

®o «n n 


(95) 


Consequently, using the notation ir = 3 rr/3x, one obtains 

3C 
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dn = n dv + n dv +n dv + n dg + n dc + 7r..d6 + n di . (96) 

v 00 o v °°n n g c o i ' ' 


oo o 


oo n 


To write the partial derivatives indicated, it is convenient to first define a factor, 
j , which is equal to zero if reference power is fixed and equal to one otherwise. 
Also, employing the notation of earlier sections, define 


_ 4c f (1 - e ) _ 

+ k rt >[ 1 + . * _ * T ] • 


2c -c (f -e ) 
r l x x 


(97) 


(where g is a new symbol similar to e and f ) and note that 

X X X 


IT 


bn 


n-l 


m dm 
o o 


= k + k + j -- ff£ c (1+k ) u + (l + k ) Y k 
s t J p 2 7) t. n t i L sampi 


n-l 

“ k t Z k drop i + g x + ^t k t^n ~ -^t k t ^ + L ^samp i k dr op i^ 

i=l i=l 


i=l 

n-l 


“U 




P PS 27) 


]• 


(98) 


Then one may write the indicated partial derivatives of n as follows for the launch 
vehicle dependent formulation: 

3 m 

n = Tt ° , (99) 

v m ov 
»o o 00 o 


2(m v -j m -j m )(1 +k ,)(1 -e ) v _ — , 

jr =i ■■-_<>-* ) t Jps ps M x ; -n r 1 

V r " 2 . 7 2ir„ ~ L . 2 . 

«n 


2c. e f v 
x x c 


(v 2 +2v 2 ) z [2c — c. (f -e ) 2 ] 
®n c r lx x 


(v " + v 2 )(v 2 +2 v 2 ) 3 
00 n c 00 n c 


?] 


\ =m o [B x (1+i t k t ) ' (1+k t )] 


( 101 ) 


m 

n = j — - (1 -i g ) , 
g > .g J psV 


( 102 ) 


36 


, ( 100 ) 


* 



„ „j !e»(i. salVi _ j g ) , 

c p C \ t? A ps x/ 


'P C \ t? 

3m 

it, = v -r— , 
o m 3 o 
o 


( 103 ) 


(104) 


3m 


o 


ff . = ff .. , 

i m 3i 
o 


aos) 


where 77 * =dr)/dc and 


3m 

0 

dm 

0 

3v 

c 

3v 

dv 

3v 

00 0 

c 

oop 

3m 

dm 

3v 

0 

0 

c 

36 

dv 

c 

36 : 

3m 

dm 

3v 

0 

0 

c 

31 

dv 

c 

31 

dm 

0 

dv 

= - W e " 


“(v /bj 


(106) 


The forms of the derivatives 3v /3v , 3v /36 and 3v /3i depend on the de- 

c “o c c 

finition of v . For relatively small geocentric declinations of V , v is given 
c 09 o c 

by equation (3). Then 


3v v 
c_ _ ”0 

v 

°°o c 


3v 

c 

36 


= 0 , 


( 107 ) 


3v 

rr~- = 0 . 
3 v 
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In the high-declination case where includes the penalties due to launch azimuth 
and non-coplanar injection, the derivatives are considerably more complex. The 
formulas for the derivatives of v^ with respect to the three independent’ parameters 
v . , 6 and i are 

00Q 


dv 


3v 


3v 3v 
os 0 «OQ 


Sv 3v 

36 ’ 


(108) 


where, from (78) 


3v c 3v 3Av. 

3i 3i 3i 


3 Av./B i = 2 c 1 i + c 2 , 


(109) 


and, from (79) and (80), 


3v /3 i = - 3 v /36 . 
g g 

The derivation of the partial derivatives 3v /3v and 3v /3 6 is straight- 

g o g 

forward although somewhat cumbersome. The equations for these partial 
derivatives .employing the notation defined earlier, are as follows: 

2 .. . 2 2 _ 2 


( 110 ) 


3v 


3v 


v r a w(3+2pw-w )(3p/3v )+(3p+2p w-3pw -4w)(3 w/3v )^ 

•f-kiP- 7== F ■ -b < 111 > 

®o 2^(l + pw-w )(2 + pw) 


00 o g 


3v 

36 

where 


2 2 2 
V q (3P +2p w - 3pw - 4w) 

ST ’ 


2v . /(I + pw - w*) (2+p w) 

g v 


( 112 ) 


3p/3v = 1/v , 

®o o 


(113) 
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X 


±!L ,i/i »£_ + + 

3v 2 1-2 5y 3v L\ _ 
®o «°o ®o 2 




)^ + (f + y )ei^ 


2 / Bv 

S “0 


+ g^- )] / (2w - p/2 - y)} , 


( 114 ) 


If - HU ♦ [H + +(f +y) U J/^w-P/2- y,} . am 


Bs/B6 = cos (6~i), 

B y/ Bv oo 0 = [(p /2 > s p/ Sv =o 0 - dx/Sv <» 0 ] /2 y» 

By/d 6= - (dx/B 6 )/2y, 


(116) 

(117) 

(118) 


§* = I j- (q/2Ma q/au ) Mp/ _ 3} .l 8£ /S . uj . + |3 (q / 2 / + (p/3) 3 +q / 2 ] 

6y y<q/2) 2 +(p/3> 3 


-2/3 


-2/3 


_ - u ) + fe/6) @ p/ 5 u ) - f^][\A q//2)2 + (p/ 3 ) 3 -q/ 2 ] . ( 11 9 ) 

6 /(q/2) 2 + (p/3) 3 


with u = v or 6, 


coq 


Sq/Sv^g 2ps 2 (i -s 2 )(Sp/5 v ooo ) » 

Bq/36 = 2p 2 s(l -2s 2 )(Bs/3 5), 


dp/Bv =2ps <3p/3v ), 

eo o ® O 


Bp/36 = 2s (p + 4 )(Ss/ 36 ). 


( 120 ) 

( 121 ) 

(122) 

(123) 


For the launch vehicle independent formulation, the partials V and 

aojj 

77^ remain unchanged from expressions (100) and (101). The remaining partials 
n 

become. 



rr 

v 


(124) 


IT* = IT. = 0 , 

0 i 

t»0 


TT 

g 


m 

— — TT 
g in 


9 



a25> 


(126) 


Of course, for the launch vehicle independent formulation, the factor j is always 
taken to be zero, since a condition of the formulation is that reference power is 
fixed. 


Now, consider the differentials dR and dR at the initial and final times. 

At launch 

dR = dP = P dt , (127) 

o o o o / 


dR = dP + dV =P dt +dV . 
o o “o o o “o 


(128) 


The differential dV may be written in terms of differentials in its magnitude 

and any two arbitrary, but independent, angles which would uniquely define the 

orientation of V . Because of its importance in the formulation, one of these angles 
00 o 

will be chosen to be the geocentric declination, 6 . For convenience, the other 
is chosen to be the geocentric right ascension, a. In terms of these angles, dV 

°°o 

may be written 


- dV = 

CD O 


e°0 


dv +(n xV )da+[(V xn )xV /v cos6]d6, (129) 
v °°0 P »o 00 O p “0 00 o 


c°o 


where n is a unit vector in the direction of the Earth’s North Pole. 

P 

At the primary target the form of the differentials will vary according to 
the target condition option chosen. For the ephemeris option 
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(130) 


: dR = P dt , 
n n n 

and, providing the velocity at arrival is not left completely open, 


dR = P dt + dV 
n n n 08 n 


(131) 


As above for dV , we may choose any two independent angles, say or and a , 
® o X 2 

and the magnitude as a set of three parameters uniquely defining V , and 

dV may then be written 
®n 


V 

oorj — — 

dV = dv + (a x V ) da, + fa x V ) da„ , 

®n v ®n ' 1 =°n 1 ' 2 ®n 2 

°°n 


(132) 


where a and a are unit vectors normal to the planes in which a and a , 

1 a 1 6 

respectively, are measured. 


For the open angle transfer option, the fixed final solar distance gives rise 
to the scalar equation. 


R • dR = 0 , (133) 

n n 


and, if the final orbit is constrained to be circular. 


^n = J ~ ^ X dR n ) + dV «n ’ a34) 

r 

n 


with the expression for dV being as given above. For the extra-ecliptic option, 

°°n 

it may be recalled that the final radius, speed, inclination to the ecliptic, and flight 

path angle are fixed. Optionally, the semi-major axis, eccentricity and inclination 

may be fixed. It is therefore convenient to write the differentials of position and 

velocity in terms of the two or three non-vanishing differentials dto , dw , and, 

n n 

when applicable, df^, where f is the final true anomaly. These are, 


dR =(kxR)dft +(hxR)(dco + df ) , ( 135 ) 

n n n n n n 
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(136) 


dR - (kxR ) dfl + (hxR )(dco +df ), 
n n n n n n 


where h is a unit vector along the final angular momentum R x R „ 

n n 

At intermediate targets, the differential of the position vector is 


dR. = dP. = P.dt. , 
1 1 li 


( 137 ) 


while, for the velocity vector, 


dR. = dP. + dV . = P. dt. + dV . 

i i t i “i 


( 138 ) 


As above for dV and dV , we write dV as follows 
OO o 00 n oat 

v »t 


dV .= dv .+ (a.,xV .) da. + (a xV J da, . 

00 1 v_. °°i 1 il 00 r il i2 »i i2 


( 139 ) 


Other optional constraints which result in variations in the form of the 
transversality conditions are the flight time, which results in 


dt = dt , 
n o 


( 140 ) 


and the reference power, p ^ , which yields , 


g 


(i-2) 

\c t) / 


dc + 


dm 


m 


= 0 , 


( 141 ) 


where, in its most general form, 

d m . dv 

dm = -r — — ( 
o d v v 


dv 
c °°o 


dv + > A 
00 o 0 0 


3v dv 
= d6 + — c dl). 


( 142 ) 


Finally, the fact that the initial mass ratio is always taken to be unity yields the 
result, 


du = 0 . 

o 


( 143 ) 


42 



The constraint equations (in differential form) are employed to eliminate 
from the general transversality equation (91) a like number of differentials. Re- 
quiring that the coefficients of the remaining differentials vanish provides the 
appropriate set of transversality conditions for a specific problem. However, be- 
cause the adjoint equations are linear and homogeneous in the adjoint variables, it 
is possible to fix the initial value of one of the adjoint variables. It is convenient 
to choose the initial mass ratio multiplier X as the independent parameter to 
be held constant, and choosing a value of unity 0 will generally result in initial values 
of the other adjoint variables also of order one. Also, the arbitrary positive per- 
formance index constant k in expression (91) will be assigned a value which causes 
the transversality condition associated with the final mass ratio to be satisfied; 

k = X /it . (144) 

v v ' ’ 

n n 

The choice of eliminating X from the boundary value problem independent para- 

o 

meters is made by the program user when the program inputs for a particular pro- 
blem are composed. Indeed, it is sometimes, but not often, advantageous to not 

hold X fixed when attempting to obtain a converged solution to a particular pro- 
o 

blem. The choice of eliminating the final mass ratio transversality condition from 
the boundary value problem dependent conditions is made arbitrarily and is accom- 
plished automatically by the internal coding within the program. These two choices 
reduce the order of the boundary value problem by one, which is considered to be 
advantageous in this particular instance (in terms of conserving computation time), 
even though it is not always advantageous to reduce the order of a boundary value 
problem, thereby forcing the boundary value problem onto a possibly highly non- . 
linear constraint subspace. 

Once these choices are made, the appropriate transversality conditions may 
be derived and written as follows : 

For open launch excess speed with m Q being independent of 6 and i the 
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transversality condition is, 


kTT 


gX 

“o e o 

— - (1-jJ 


J p' X m Bv 0 


-1 = 0 


(145) 


o o 


For open launch excess velocity direction with initial mass independent of 6 and i: 


A x V = 0. 
o “o 


(146) 


This implies V || A ; V aligned with A is usually the proper choice. 

oo O O “O O 

For cases in which m Q is a function of 6 and/or i, the transversality 
conditions are modified as follows, where 


f = [k +k A -(l+kJf -gX /km ] dm /dv . 
st t n g o oc 


(147) 


Open excess speed: 


f(B v /Bv ) - (A * V )/v = 0 . 

g “o o “o °°o 


(148) 


Open geocentric declination of V : 

oo o 


f(Bv /8 6)-A ' [(V xn )xV /v_ cos6] = 0 
g o ' “o p °°o °°o 


a49) 


Open right ascension of V : 

“O 


-A • (n xV ) = 0 
o p ®o 


(150) 


Open launch parking orbit inclination: 


f(d Av /Si - 9 v /B 6) - 0. 

i g 


(151) 


Note that these last two equations are functions only of initial conditions*; hence, they 
may be solved directly for a and i, which therefore need not be included explicitly 
as independent parameters. Equation (150) dictates that a be equated to the geocentric 


♦Assuming f / 0. 
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right ascension of A , or 180 degrees therefrom. Equation (151) may be solved 
for i; however, the expressions are sufficiently complex that the solution must be 
obtained iteratively. 

For open excess velocity direction at an intermediate target, the condition is. 


(A. - A. - ) x V . = 0. 

i i “i 


(152) 


For open excess speed at an intermediate target, the primer is continuous, i.e., 


A. - A." = 0. 

i i 


(153) 


For open excess velocity direction at the final target, the transversality condition 
is 

A x V =0. (154) 

n «n 

The proper direction of is usually opposite A^. 

For open arrival excess speed (at the primary target) in problems where a retro 
stage is employed, the transversality condition is. 


k?T 


5 n 


-1 = 0 . 


(155) 


n 


For open thrust acceleration with unspecified reference power in the launch vehicle 
dependent formulation, the transversality condition is, 

k?r 


+ 1 - 0 . 


(156) 


g 


For open jet exhaust speed with unspecified reference power in the launch vehicle 
dependent formulation, the transversality condition is, 


kff 


+ 1 = 0 . 


(157) 
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These last two equations also apply for open thrust acceleration and open jet ex- 
haust speed, respectively, in the launch vehicle independent formulation with 
fixed reference power. Of course, the appropriate expressions for the partials 
of ff must be used, however. 


If the reference power is specified using the launch vehicle dependent 
formulation, but both reference thrust acceleration and jet exhaust speed are 
left open, the last two transversality conditions are replaced in favor of the one 
condition. 



(158) 


In the preceding equations X and X are evaluated at time t , the time at 

g c n 

which the spacecraft is to be at the primary target. When using the open angle 
transfer option, the transversality condition associated with the open angle is 

TA xR -A xR ]• k = 0. (159) 

L n n n nJ 

For either the ephemeris or open angle options, if the final velocity is completely 
unspecified, as in the case of flyby missions, the appropriate vector transversality 
condition is, 

A r = 0. (160) 

This causes the primer locus (trajectory in primer-space) to home-in on the 
origin of primer-space, which is a singularity of the optimal rocket problem 
(as previously mentioned) if the spacecraft should be thrusting at that point. 

However, if forced-thrusting is not invoked via the propulsion-time multiplier 
X T , condition (160) tends to cause a coast phase to occur at the time of flyby, 
thus avoiding the singularity, A similar situation can occur at the initial time 
for the launch vehicle independent mode if launch excess speed is to be optimized. 
This leads to the condition 

A q = 0. (161) 
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For the extra -ecliptic option, the fact that Ci and to are left open in specifying the 
final state gives rise to the two conditions, 


° k = 0 , 

(162) 

• h = 0, 

(163) 


where C = RxA-RxA is the vector constant of the motion on a given 
o o o o 

trajectory segment and h is a unit vector along the angular momentum of the 
final heliocentric orbit. If the option is employed which additionally leaves final 
true anomaly open, the extra condition, 

~ (A* R) + r 2 (A • R) = 0 , (164) 

applies. The transversality condition associated with open launch date for the 
ephemeris option is , 

.. 

-A “P+A • P +h =0, (165) 

o o o o vo 

the condition for open encounter date at an intermediate target is, 

- (A. + - A.“) • P. + (A. + - A.”) • P.+ h + - h 7=0, (166) 

l 111 i ivivi 

while the appropriate condition for open arrival date (at the primary target) with 
the ephemeris option is, 

A • P - A * P - h =0. (167) 

n n n n vn 

Since the variational Hamiltonian h is a constant of the motion on a given trajectory 

v 

segment, the time at which it is evaluated on that segment is arbitrary. The preced- 
ing conditions pertaining to initial and final time are applicable if the total flight time 

is unconstrained. In the event that the total flight time is fixed while both t and t 

on 

are left open, the two preceding conditions are replaced with the single condition re- 
presented by the sum of the two, i.e. , 
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(168) 


A • P - A ’ P -h -A ' P + A • P +h =0. 
nn nnvn oo oovo 

For either the open angle or extra ecliptic options, the appropriate transversality 
condition associated with open arrival date (i.e. , open flight time) is 

- h y = 0 , . (169) 

where the signs (+ and -) in all the preceding expressions are displayed as they 
are coded in the program,, 

Finally, the transversality condition associated with unspecified, but con- 
stant, thrust cone angle is, 


X, (t ) = 0. 
0 n 


(170) 


5. Partial Derivatives. In solving the two-point boundary-value 
problem, the finite difference method is used for computing the required matrix 
of partial derivatives of the dependent variables (i.e., end conditions or constraints) 
with respect to the independent variables. The finite difference method involves 
individually perturbing each independent variable x. to be optimized and computing 
an associated perturbed dependent parameter vector Y . The amount each inde- 
pendent parameter is perturbed, Ax^ , is specified by program input. The desired 
partial derivatives are then approximated as follows: 



(171) 


where Y represents the nominal dependent parameter vector. The finite difference 
method is completely general and is applicable to all end condition options available 
in the program. 
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C. EXTENSION OF SOLUTION FOR POWER DEGRADATION. 

The assumed spacecraft and trajectory models are as described earlier in 
Section A and are not repeated here ; the nomenclature used in the following 
analysis is also described in Section A, except for the introduction of new symbols 
which are described in the text as they appear. 

Historically, this electric propulsion power degradation model first appeared 
in the literature in [ 7 ] , and then soon afterward an improved discussion appeared in 
[8], in which several of the ramifications and consequences of the theory are also 
discussed. For deeper insight into the analysis, the reader is therefore referred 
to [8] , from which much of the analysis below is extracted. 

The model discussed here is concerned with the manner in which the per- 
formance of a solar array degrades due to high energy particle damage. It is 
assumed that one can define a damage factor q which has a value of unity at the 
initial time and which decreases in value with time during the course of the mission 
such that the power output, p, of the arrays at any time may be written 

p = q 7 p ref» (172) 

where p is the reference power and y is the power factor which is a function 
of solar distance and array orientation relative to the sun. The damage factor q 
may also be thought of as a time-dependent efficiency. The derivative q is nega- 
tive, and is assumed to be linearly proportional to q and to the density of particles 
impinging on the solar cells. For simplicity, it is assumed that the particle density 
d is of the same form as the density of photons striking the surface of the array, 
i.e., 

e • n 

d = “V" ’ (173) 

r 

where e is a unit vector along the sun-spacecraft line, n is a unit vector normal 

r 

to the arrays such that e^ • n > 0 implies the front of the panels faces the sun, and 
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r is the solar distance of the spacecraft. Then 


and 



- kqd, 


(174) 



( 175 ) 


where k is the constant of proportionality. It is convenient to introduce a para- 
meter s, called degradation time, 


or 



ddt'. 


s = d . 


(176) 

(177) 


Thus, under the above assumptions the degradation time accumulates at a faster 
rate when the spacecraft is nearer the sun, which is a characteristic one might 
expect. Note that when the spacecraft is maintained at 1 AU with the panels 
normal to the sun line, s = 1 and the degradation time is equal to the flight time. 
Then k may conveniently be thought of as the inverse of a reference time, called 
the characteristic degradation time denoted by T^, and the degradation factor q 
becomes 

-s/t. __„ l 

q = e d. (178) 


Actually, there is little reason to allow the degradation time s to continue to 
accumulate during coast phases since the arrays may be turned edgewise to the 
sun and the degradation process may be effectively halted. Therefore, we adjust 
equation (177) to read 

s* = h^d , (179) 

where h^ is a step function equal to one during thrust phases and equal to zero 
during coast phases. 

The characteristic degradation time Tj is an engineering parameter that 
must be determined experimentally. For example, by exposing a solar cell to 
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the particle emission of a solar simulator and measuring the performance of the 
cell over a period of time, one should be able to estimate a reasonable value of 
Tj. Another source of information would be measurements from actual spacecraft 
which employ solar cells for power supply. 

The assumed exponential form of the degradation factor, although intended 
for use with SEP systems, is applicable for NEP systems as well. The principal 
difference is in the definition of s ; for example, s = h^. The exponential form 
also permits one to evaluate radio-isotope systems by defining s = 1 and letting 
T\j represent the time for the radioactivity to dissipate to l/e of its initial level. 

In the development which follows, the formulation applicable to SEP is used 
exclusively. The equations of motion and adjoint equations are given in Section A; 
the equations affected by the inclusion of power degradation in the model are given 
below. 


Power degradation affects the problem in a very fundamental sense, beginning 
with the rocket-thrust term in the equations of motion: 


V = R = a e t - R , 
r 


« 

R = V , 


( 180 ) 


where R is the position vector, r is the magnitude of R, V is the velocity 
vector, M is the gravitational constant of the sun, a is the magnitude of the thrust 
acceleration and e^ is a unit vector in the direction of thrust. The thrust accelera- 
tion a is a function of several variables and may be written as follows: 


a = h 


gyq 


( 181 ) 


where g is a reference thrust acceleration evaluated under a prescribed set of con- 
ditions, v is the ratio of current to initial mass, and q is the degradation factor 
defined above. The power factor y is assumed to be a function of the density, d. 
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of photons incident on the arrays, where d is as written in equation (173). The 
general form of y is, 

4 

y= d£ a.d l/4 . (182) 

i=0 

The coefficients a^ are chosen so that this equation will adequately describe the 
power output of a given array. The only restriction of the a. is that their sum 
is equal to one. Then at r = 1 AU, with the arrays normal to the sun line, 
y= d = 1. The summation term in (182) represents the temperature effect. 


The mass ratio satisfies the differential equation, 


v = - h 


gyq 


a c 


(183) 


using V = 1 as an initial condition, where c is the jet exhaust speed which is 
assumed to be constant over the trajectory. 

The variational Hamiltonian becomes, 

h=A <, R-A-R + Xt'+XT+Xs’+Xg+Xc 
v v r s g c 

= h o[^ <A- Vc V + V + X t]- ^(A- R)-A -R, 

r 

and the adjoint equations are 

A.-JL A+ aB p . A)R+h [«Ca (A . ;t .K v + x ]|±. 

r r 


(184) 


(185) 


X* = - h -2 s - (A • e - - X ), 
g a v t c v 
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where 


X = 


c 



gyq 

2 

c 




= h 


ZXSl. 
a i /r, 
d 


(A 


vfv. 


4 

y* = dy/dd=V a. (i+“J-)d 1 ^ 4 , 
i=0 



ass) 

cont. 


The control variables are the thrust direction e , the switch step function 
h , and, providing the array orientation is not constrained to yield maximum 
power, the normal direction n_ According to the Maximum Principle, these con- 
trols are chosen to maximize the variational Hamiltonian (184). The maximum of 
h^ with respect to h is seen to depend totally on the sign of the term in square 
brackets. That is, denoting 


= ^Z£_ ( a. e t -|x^) + X s d + X r , 

then choose 

_fl if a* >0 
a to if cr* < 0 


(186) 


(187) 


Maximizing h^ with respect to e t is also accomplished by inspection. Since 
appears only in the dot product with A and since the coefficient of that dot product 
is non-negative, i.e., 

h gyq/V* 0 , (188) 

a 

then h^ is maximized with respect to by making A • e as large as possible, 
which is the same result as when there is no degradation. The control vector n 
appears explicitly in (184) through the density d as given by equation (173). In 
fact, since n appears in h^ only through d, n affects h^ only through its dot 
product with e . Letting the angle between n and e r be denoted X such that 
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e * n = cos X > then it is clear that there will be a "best" angle X between n 
r 

and e^ to maximize h^, but that n may lie along any element of a right circular 
cone of half angle X about e^. For the moment, we will put aside the question of 
the explicit direction of n and concentrate on defining the optimum x or » alterna- 
tively, the optimum d„ The optimum value of d is determined by maximizing h^ 
with respect to d, i.e. , by solving for the root of the equation 


3h 

v 

dd 


h 

C 7 



(189) 


Performing the indicated differentiation yields 


or 


v 


cr + X 
r s 


= 0 , 


X V 
s 

* 


(190) 

(191) 


where, using (182), 


and 


4 

y*=f 

i=0 


a. 

i 






(192) 


(193) 


i/4 

Because of the form of (192) equation (190) is a quartic in the variable d , and 
is solved by iteration in the program. A more detailed discussion of the solution 
to (190) is given in Reference [8] . for a specific set of coefficients a f . For now, 
assume that the optimum value of d is found from (190). Then the optimum angle 
X is immediately obtained 


cos X = d r 



(194) 


Of course, equation (190) does not take into consideration the fact that d can never 
exceed the inverse square of the solar distance. Consequently, the right hand side 
of (194) may exceed unity, under which condition the program sets cos X = 1 (i.e. , 


54 


X = 0) implying that n is directed along e^. 

If the d that represents the solution to (190) exceeds the upper limit allowed 
for d whether that limit is imposed by the problem or by nature (i.e. , l/r ), the 
correct choice for d is that upper -limiting value. Likewise, on the lower side, d 
is physically limited to be non-negative. Therefore, a negative solution to (190) is 
disregarded, and d is set to zero which corresponds to X = tt/ 2 (panels oriented 
edgewise to the sun), and the engines are shut down. 


The precise definition of n has no bearing on the solution of the problem, 

except as it affects d as defined in (173). The appearance of n in the state and 

adjoint equations and the variational Hamiltonian is solely through the density d 

except in the equation for A where n appears explicitly as part of the partial 

dd/dR, defined following equations (185). Actually, this partial is valid only if 

d is permitted to vary with R. That is, if either d = 0 or d = constant is imposed 

then dd/dR becomes the null vector and the entire term drops from the equation 
♦ • •• 
for A . Furthermore, if d is the solution to (190), then the last term of A in 

(185) again drops out because the term in square brackets is the left side of (190). 

• • 

Therefore, the only time the term in question remains in the equation for A is 
when d = l/r^ which corresponds to cos X = * n = 1 and implies n = e^. 

Under this condition 


dd = _2_ ~ 

SR " r 3 6 r * 


(195) 


The boundary condition pertaining to the initial degradation time is 

s(t ) = 0 . (196) 

o 


If the final degradation time, s q , is unspecified, the transversality condition 
associated with degradation time is 

A =0. (197) 

s_ 
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The initial value of A ^ is unknown and therefore becomes one of the independent 
parameters of the boundary value problem. From equations (185), it follows that 


x* =-- s - 

S 7 d 



(198) 


and, therefore, 





Using the boundary conditions X = X = 0, it follows that 

s a 
n 6 o 


(199) 




( 200 ) 


Thus, if an approximate value of A is available from a trajectory similar to the 

g n 

one of interest, a reasonable guess of A is easily approximated. 

s 

o 
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D. AUXILIARY COMPUTATIONS 

The following paragraphs present the equations employed in a variety of 
auxiliary computations throughout the program. The computations are termed 
auxiliary because they do not influence the optimization results or procedures, 
with the possible exception of the Spiral Capture computations. Rather, the com- 
putations are made for printout purposes. 


1. Standard Block Print Variables . A standard print block is 
employed for printing information at various points along a trajectory. Each 
standard block contains a total of 40 parameters as follows : 


TIME 

SEMI-MAJOR AXIS 


ECCENTRICITY 


Time since departure, in days. 

Semi-major axis of the osculating heliocentric trajectory. 

a = 1/ (f'7r)' lnAU - 

Eccentricity of the osculating heliocentric trajectory. 


e = Jl - h 2 //ia , where h is the magnitude of the 
angular momentum. 


INCLINATION 


NODE 


Inclination of the osculating heliocentric trajectory re- 
lative to the ecliptic. 

i = cos *(k • H/h), in degrees where H is the angular 

© 

momentum vector, H = R x R . 

Longitude of ascending node of the osculating heliocentric 
trajectory plane on the ecliptic, measured eastward from 
the vernal equinox. 

sin 1 £j • (kxH/h)J=cos *£i • (kxH/h)J, in degrees. 
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ARG POS 


RMAG 


TRAVEL 


Rl, R2, R3 
VI, V2, V3 
MASS RATIO 
THRUST ACC 

LI, L2, L3 
L4, L5, L6 

L7 

HAM 

LG 

LC 

LPHI 


Angular position in osculating orbit plane from ascending 
node measured in direction of motion. 


C0= cos 


IJ-Rj. (k x H) j _ s . Q -l j~ R ♦ [Hx (kxH)] 


L rh sin 
in degrees. 


rh sini 


]• 


Magnitude of the instantaneous heliocentric spacecraft posi- 
tion vector, r = |r | , in AU. 

Approximate travel angle since launch, in degrees. 

- . R. • R . 

0 = E. 0 . where 0 . = cos ( — — — ) 

t i i i \ r. r. , / 

i i-l 

tVi 

and R. is the spacecraft position vector at the i 
computation step. 

Cartesian components of current position vector, in AU. 
Cartesian components of current velocity vector, in EMOS. 
Ratio of current mass to initial mass. 


Ratio of current thrust acceleration to current solar 

2 . 

gravity acceleration (=h^gyqr /\iv), 

Cartesian components of primer vector, A . 

Cartesian components of time derivative of primer vector, 

A. 

Variable adjoint to the mass ratio, A^. 

Variational Hamiltonian, h . 

v 

Variable adjoint to the thrust acceleration, A^. 

Variable adjoint to the jet exhaust speed, A q . 

Variable adjoint to the thrust cone angle, A^. 
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CONE, CLOCK 


HMAG 

POWER FNCT 
SWITCH FNCT 
PSI, THETA 


PHI 

LATITUDE 

LONGITUDE 


Two angles, in degrees, defining the angular position of 
Canopus relative to a spacecraft fixed coordinate system. 
Employing a built-in unit vector s in the direction of 
Canopus, the two angles are defined 


CONE = cos * (s • e ) 

s 


CLOCK = cos 


_l f (e xe t ) * 


L lz ~ 


VV 


<v 5) 1- sin -i r (e s xg t )x(i s xg) 

|e g xii J L |e g xe t | |i g xi| 



where e is a unit vector pointing from the sun to the spacecraft, 
s 

2 

Magnitude of the angular momentum vector, in AU /tau. 

♦ 

h = |Rxr|. 

The power function y q defined in an earlier section. 

The thrust switch function a . 

Thrust angles relative to the instantaneous osculating 
trajectory plane. PSI is the out-of-plane angle and THETA 
is the in-plane angle, in degrees. 

0 = sin 1 (e t * H/h) 

„ .-If- HxR \ -1 - 

8= sm e^ • r J = cos (e, • R/r cos W) 

\ t hr cos ip / t 

Thrust angle relative to the sun-spacecraft line 

0 = cos 1 (e^ • R/r), in degrees. 

Ecliptic latitude of spacecraft, in degrees 
= sin 1 (z/r). 


Ecliptic longitude of spacecraft, in degrees 


= sin 1 (y/Jx 2 +y 2 ) = cos 1 (x/Jx 2 + y 2 ) . 
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FLT PTH ANGLE 


Heliocentric flight path angle, in degrees 

-1 

= sin (R • R/rv). 

VMAG Heliocentric speed, in EMOS 

v = |r| . 

PROP TIME Total time propulsion system has operated, r, in days. 

The above standard block may be augmented in two ways. When power degra- 
dation, as indicated by the input variable TPOWER, is simulated, a single line of 
information is automatically added to each block, as displayed in the Sample Pro- 
blems and Results section, Case G, Comet Rendezvous Mission. When the input 
variable MPRINT is 2 or 3, three extra lines of information are generated per 
block, as displayed in the Sample Problems and Results section, Case H, Multiple 
Ballistic Swingby Mission. These two types of additionally printed lines may appear 
simultaneously. 

(a) . Power Degradation. The single line of power degradation in- 
formation contains eight parameters as follows: 

S Degradation time, s, since departure, in days. 

LS Degradation time adjoint variable, X . 

S 

-2 

DENSITY Density parameter, d, in AU . 

DPOWR qdy/dr, in AU _1 . 

DPOWD qdy/3d. 

DEGRAD The degradation factor, q. 

CHI Solar array orientation angle X » in degrees, 

CHI REF Solar array orientation angle which the arrays would have 

if oriented for maximum power, in degrees. 

(b) . Target-Relative Coordinates and Comet Magnitudes. The 

three extra lines which may appear via using MPRINT contain the following information 
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R1 REL 
R2 REL 
R3 REL 


Cartesian components of current spacecraft position vector, 
with respect to the next astronomical body to be encountered 
along the trajectory in a moving coordinate system generated 
by that body, with the x-axis pointing outward along the 
body’s heliocentric radius vector, the y-axis in the body's 
orbit plane in the sense of the body's motion, and the 
z-axis completing the right-handed orthogonal system, in 
kilometers, with the origin of coordinates at the body. 


VI REL 
V2 REL 
V3 REL 


RMAG REL 
VMAG REL 
S/C NUC MAG 


S/C TOT MAG 


Cartesian components of current spacecraft velocity vector, 
in kilometers/second, in the target-relative coordinate system 
described directly above (see R1 REL). 

Magnitude of R1 REL, R2 REL, R3 REL, in kilometers. 

Magnitude of VI REL, V2 REL, V3 REL, in kilometers/second. 


Nuclear magnitude (of comet) of the next astronomical body to 
be encountered along the trajectory, as seen by the spacecraft. 

M = M + M, log |r -R | + M log |R | 

N o 1 6 10 1 targ 2 & 10 targ 1 


+ .03 cos 


-1 


R 


• <R* 


targ targ 


-R) 


|R | |r. -r 

1 targ targ 


] c °- 


where M , M, , and M„ are magnitude constants associated 
o 1 2 

with the target, and C° is the radians-to-degrees conversion 
factor. The arc-cosine term is the phase angle. 

Total magnitude of the next astronomical body to be en- 
countered along the trajectory, as seen by the spacecraft. 

M = M +M log |r -R | + M log |r | 

T 3 4 to 10 ' targ 1 5 10 targ 1 

where M„, M„, and M_ are magnitude constants associated 
3 4 5 

with the target. 
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GEO NUC MAG 
GEO TOT MAG 
ANG(V, R) 


ANG(V, XY) 


R1 REL ECL 
R2 REL ECL 
R3 REL ECL 

VI REL ECL 
V2 REL ECL 
V3 REL ECL 


Same as S/C NUC MAG, except as seen by the Earth. 

Same as S/C TOT MAG, except as seen by the Earth. 

Angle which (VI REL, V2 REL, V3 REL) makes with the 
positive x-axis in the target-relative coordinate system des- 
cribed under R1 REL, in degrees. 

ANG(V, R) = cos" 1 (VI REL/VMAG REL) 

Angle which (VI REL, V2 REL, V3 REL) makes with the xy 
plane in the target-relative coordinate system described under 
R1 REL, in degrees. 

ANG(V, XY) = sin" 1 (V3 REL/VMAG REL) 

Same as R1 REL, R2 REL, R3 REL except expressed in the 
ecliptic coordinate system of date. 

Same as VI REL, V2 REL, V3 REL except expressed in the 
ecliptic coordinate system of date. 


RMAG ECL Magnitude of R1 REL ECL, R2 REL ECL, R3 REL ECL, 

in kilometers. 


VMAG ECL 


Magnitude of VI REL ECL, V2 REL ECL, V3 REL ECL, in 
kilometers/second. 


2. Extremum Point Summary Print . For the final trajectory of 
each case, a table of extremum points of selected functions occurring along the 
trajectory is evaluated and printed. This table contains the following variables; 
the number of iterations required to isolate the given point, time from launch, 
ecliptic longitude, solar distance, communication angle and distance, thrust switch 
function, the three thrust angles, (!), 0, and 0 , the power input to the thrust sub- 
system, and the array orientation angle x • Extrema are evaluated for all of these 
functions except the time from launch, ecliptic longitude, and array orientation 
angle. In addition, the table includes the initial and final times, all thrust switch 
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points, and the points at which the power profile changes under Options 4 and 5 
of the power function selector. The values of all functions contained in the table 
are printed each time an extremum or special point of any function is encountered. 

Extrema points are evaluated by locating points at which the derivatives of 
the functions go to zero. The other special points (e.g. , thrust switch points) 
are obtained by defining a function which goes to zero at the special point (e.g. , 
the thrust switch function cr) and isolating these roots. A general set of routines 
is employed for determining all the extrema and special points. The procedure 
is to check at each computed point whether any root of interest has occurred 
within the current (most recent) computation step. If a root is known to exist, 
an iteration scheme is initiated to isolate the point and store the related informa- 
tion. This is undertaken on all trajectories to define the special points and the 
extrema of the thrust switch function. The extrema for all other functions are 
obtained only for information purposes, however, and therefore are evaluated 
only for the last trajectory of each case. 

The ecliptic longitude printed in the table is referenced to the longitude 
at launch; consequently, the value printed at the initial time is zero. There- 
after, the longitude is accumulated in increments so that the value will exceed 
360 degrees if the spacecraft traverses more than one revolution of the sun. 

The angle is printed in degrees. 

The solar distance is simply the magnitude of the spacecraft's helio- 
centric position vector and is printed in AU. Extrema of the function are ob- 

• 9 . 

tained by isolating those points where r = R • R/r is zero. 

The communication distance is defined to be the distance from the Earth 
to the spacecraft and is evaluated in AU as follows: 

r = |R - P 1 , (201) 

C 6 
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where P is the heliocentric position of the Earth, 
e 

evaluated by locating the points when 


The extrema in r are 
c 


• • 

(R-P )• (R-P ) = 0 . (202) 

6 6 

The communication angle is the angle in degrees subtended at the Earth between 
the Earth-spacecraft and Earth-sun lines. It is evaluated as follows: 


a 

c 


cos 


-1 


[ 


P * (P -R) 

|pj lv R l 


]• 


except at the initial time when R = P ; then, 


a 

c 


cos 


-1 


P • (P -R) 


P -R 
e 


■/ - 

ri J* 


(203a) 


(203b) 


Extrema of this communication angle correspond to the roots of the equation. 


[<R-P e ,xp-P e )] 


x (R - P ) 
e 


|pj IB-Pj 


P + 


(P X P ) X P 
e e e 


e . l . .3 ( R - p J = 0. (204) 

e |R - P | |P r e 


The above equations for communication distance and angle assume that the input 
parameter NDIST is set to its default value of 3, which corresponds to Earth. 

Any other celestial body may be chosen as the reference for the calculations, 

• 

through proper choice of NDIST, in which case the vectors P and P then 

e e 

represent the position and velocity vectors,, respectively, of the corresponding 
reference body. 


The thrust switch function is evaluated as defined earlier. Extrema of 
the thrust switch function correspond to the roots of the equation*, 

q = A » e •+ A • e ,-h ^ a - (R • R) A = 0 , (205) 

t t a vc 2 t 

& y r 

*In the absence of power degradation; q = 0. 
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with 


e = " i — ■ [(Ax A) x A], 
‘ |A I 3 L J 

if thrust direction is not constrained, and 


e =e cos0 + (mxe +mxe ) sin0 , 
t r r r 


if thrust cone angle 0 is constrained where 


e = R/r, 
r 


e r = -g C(R x R)xR] , 


e x = A/X, 


- -j- [(Ax A) x A ] , 

A 

R x A 


|R x A| 


e xe, 
r A 


x (e^ xe^ + e r xe^)xmj 


The extrema of the thrust angles 0 , 0 and 0 are defined by isolating 
the points at which their time derivatives vanish. For the case of unconstrained 
thrust angles, these derivatives are defined: 

• i - _ _ - 

$ r (e , * e + e * e , ), 

cos 0 t h t h 

• i - - - - 

. 0 = - — ■ (e , * e + e, • e ) , 
sin0 t r t r 


x 

0 = - : — x r~ (sin 0 0 - cos 0 sin 0 lb), 

sin 0 cos 0 


x r- (e • e + * e + sin 0 sin ib ft) , 

cos 0 cos 0 ' t v t v 


65 



where e ^ is as defined above for the unconstrained thrust angle case, and 

e, = (R x R)/ jR x R I = H/h, 
h 

i h = ^[(HxH) X H ] , 
h 

H = h^-^— (Rxe t ), (210) 



The two expressions for 0 are employed to avoid singularities when either 
sin 0 or cos 0 vanish. 

For the case of constrained thrust cone angle, no extrema of the cone 
angle 0 is sought since the angle is a constant. The derivatives of the other 
two angles are as defined above using the expression for e appropriate to the 

o U 

case of constrained cone angle and setting 0=0. 

The instantaneous power input to the thrust subsystem, in kilowatts, is 
also listed. Extrema of this function are isolated by locating the roots of the 
function d(yq)/dt. 


3. Swingby Continuation Analysis . Auxiliary computations are 
optionally provided, invoked by the NAMELIST input vector MOPT4, whereby 
ballistic swingbys past the primary target may be simulated. 

In one mode of program operation, invoked by MOPT4(l) > 0, single 
swingbys past the primary target may be simulated to up to ten post-swingby 
targets per case. 
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In another mode of program operation, invoked by MOPT4(l) < 0, multiple 
swingbj's along a single trajectory may be simulated, first swinging past the 
primary target and then subsequently swinging past more targets downstream 
along the trajectory. One multiple swingby trajectory may be simulated per 
case. 


In either mode of operation, the following basic assumptions are made. 
The swingby continuation computations are independent of the trajectory leg 
leading up to the swingby target, which may consist of an optimized electric 
propulsion trajectory segment (if the swingby planet is the primary target), 
except that the arrival and arrival time at the swingby planet are used in the 
determination of the swingby passage conditions. Each swingby maneuver is cal- 
culated under the assumption of the patched-conic approximation, and the swingby 
planet's sphere-of-influence is assumed to have zei 4 o radius as seen from inter- 
planetary space and infinite radius as seen from the planetary vantage point. 

The passage time in the swingby planet's sphere-of-influence is neglected (taken 
to be zero in the heliocentric frame). 

Each swingby maneuver may be either unpowered or powered, and these 
two cases are discussed in the following sections. Since the unpowered swingby 
solutions are embedded in the wider class of powered-swingby solutions, tending 
to appear in pairs which are separated by a region of braking powered swingb ys, 
the more general case of powered swingbys is discussed first. 

(a). Powered Swingbys. This type of swingby maneuver is re- 
stricted to occur at the mutual perifoci of the approach and departure hyperbolic 
arcs; the powered phase is impulsive and the thrust is colinear (pro or con) to 
the velocity at closest approach. Whether the swingby is powered or unpowered, 
the trajectory segment leading up to the swingby planet has been pre-determined, 
this being the method by which the program has been designed to obtain swingby 
solutions. Therefore the swingby time and the arrival hyperbolic excess velocity 
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are know. In the following analysis, subscript A pertains to arrival at 
the swingby planet and subscript D pertains to departure. 

A basic assumption of the powered swingby problem posed here is that 
the flight time from the swingby planet to the next target is specified. This being 
so, the program is able to converge, by iteration, on some ballistic trajectory 
from the swingby planet to the next target having the specified transfer time, im- 
plying that the departure hyperbolic excess velocity V ^ at the swingby planet 
is thereby determined. Therefore, the heliocentric trajectory before and after 
the swingby planet is determined, and it then remains to perform the required 
computations pertaining to the hyperbolic arcs within the swingby planet's sphere 
of influence. 


and 


The closest approach distance is found by iteration as follows. Let 


a = 1 + 
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where = I | , v^ = jv^ | , r^ is the (unknown) passage distance, 
and H is the swingby planet's gravitational parameter. Then the approach and 
departure hyperbolic bend angles are given by 


A -1 . -1 

— = cosec a = sin 


D -1 . -1 

— • = cosec a = sm 

2 D 


and these must sum up to the total bend angle, 


V . and V „ 
°°A “D 


a/v 

a/<y , 

which is specified in terms of 
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Therefore, using r^ as the independent variable, the zero of the quantity 


F = sin 1 (l/a ) + sin 1 (l/a )-cos 1 
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is obtained by Newton's iteration, using the derivative, 
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When the iteration is converged, the passage distance r^ is in hand, and the im- 
pulsive velocity increment is computed, 


Av = 





(217) 


where the square-root-quantities are the hyperbolic speeds at closest approach. 
The remaining parameters defining the planetocentric transfer are computed as 
follows. The inclination of the swingby orbit plane to the planet's equator is given 

by 


i = cos 1 (h • n ) 
v P 


(218) 


where h is the unit vector along the angular momentum of the hyperbolic passage 

trajectory and n is a unit vector pointing toward the swingby planet's north pole. 
P 

The ascending node angle of the swingby orbit plane is computed as, 


Cl~ tan * (-h /h ), 
.x y 


(219) 


and is placed in the proper quadrant by using the system library routine DATAN2. 
The argument of perifocus is given by, 
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where is the unit vector pointing toward the closest approach point and r 
is the unit vector lying along the line of nodes and pointing toward the ascending 
node. This is adjusted for the proper quadrant by the test, 

If h (r x r ) < 0, w -» 2ir - to. 
z n p z 

In the right-handed planetary reference frame, the z-axis is toward the planet’s 
north and the x-axis points toward the ascending node of the planet’s equator on 
the ecliptic . 


(b) Unpowered Swingbys. This type of swingby maneuver is 
considered to be a powered swingby having Av = 0. The program adjusts the 
post-swingby heliocentric trajectory segment, by iteration, until the swingby 
departure magnitude equals the given arrival magnitude. The primary 
independent variable in this iteration is the post-swingby transfer time to the 
specified target, which was held constant in the powered swingby case. Thus 
v ^ = v = v , and the swingby passage distance is obtained from the formula, 

GO I j 00 ^ CO 
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p 



2v, 





(221) 


The other orbital parameters are obtained from the same relations given above in 
the section, Powered Swingbys. 

The program can generate multiple-revolution ballistic arcs, and a parti- 
cular solution obtained by the program may not be unique, even for the same trans- 
fer time. All solutions are reachable, however, by means of inputting an approp- 
riate initial velocity guess for the trajectory segment in question. 


4. Spiral Capture . The program provides the option of computing 
approximate performance requirements of an electric propulsion spiral capture 


maneuver at the primary target planet. The approximation is based on asymptotic 
matching techniques developed by Fimple and Edelbaum (Reference [9] ) and by 
Breakwell and Rauch (Reference [lO]). The technique assumes that a helio- 
centric trajectory to a conceptually massless point with position and velocity 
of the primary target planet is available. The approximation then yields the 
additional propellant and propulsion time that would be required above that of 
the heliocentric trajectory to insert the spacecraft into an orbit of periapse r^ 
and apoapse r using the electric propulsion spiral maneuver. It should be 

H. 

noted that the additional propellant and time computed in this approximation does 
not represent the propellant and time spent performing the spiral with very high 
accuracy because the heliocentric trajectory included a trajectory segment which 
was within the geometric boundaries of the sphere of influence of the planet. The 
additional propellant and time computed is more appropriately considered a cor- 
rection to the heliocentric trajectory which, when added to the requirements of 
the heliocentric trajectory, yields a good estimate of the total performance re- 
quirement, including those of the spiral. 


Defining the semi-major axis of the capture ellipse, 


a = (r +r )/2, 
c a p 


( 222 ) 


and the thrust acceleration at the end of the heliocentric trajectory 


a - gyq/i', 

n n 


(223) 


then an incremental velocity associated with the spiral maneuver is calculated, 

j 2 

/ V. P . a a . 

Av = V — fl-l-8 *(— VI, . (224) 


which leads to the additional propellant 


a ,, _ -Av/c 

Am = m V (i-e ), 
p on 


(225) 
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and the additional time 

At = c(l - e” Av//c )/a . 

n 

This option is invoked with the input parameter MTMASS. 
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E. SPECIAL PB OGHAM FEATURES 

The following paragraphs describe several special features available in 
HILTOP which are provided to alleviate certain numerical difficulties and to in- 
crease the program's generality and flexibility* 

1* Perturbation Step Size Selector . The program generates a 
partial derivative matrix of the dependent variables with respect to the independent 
variables by integrating trajectories neighboring the current nominal trajectory. 
The perturbation step size of each independent variable, which is used to vary 
that lone variable in order to generate its associated neighboring trajectory, may 
be input to the program. Each neighboring trajectory is used to generate one 
column of the partial derivative matrix, each element of which is constructed 
by forming the simple ratio of the difference between the neighboring and nominal 
dependent variable values to the perturbation step size. Each element of the 
matrix thus constructed represents the secant-slope-approximation to the actual 
dependent variable slope, and this approximation comes closest to the true value 
for some unique perturbation step size value for each independent variable. The 
approximation becomes poor for large step sizes because the secant -slope deviates 
farther from the true slope, and becomes poor for very small step sizes because 
the numerical accuracy of the computer and also of the trajectory generation 
algorithm with its numerous iterations introduces computational noise. 

A program option controlled by the input variable KPART is available 
which attempts to determine the optimum perturbation step size for each inde- 
pendent variable. The program accomplishes this by taking a linear walk in 
the base-10 logarithm of each independent variable step size starting from the 
input or default value and not exceeding KPART steps. The program first steps 
in each direction (smaller and larger step size) to determine the proper direction 
of the walk, and each step consists of varying the step size one-half order-of- 
magnitude. For each linear walk (for each independent variable), that one column 
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of the partial derivative matrix associated with the independent variable is com- 
puted for each step of tlie walk, and each element of that column is compared be- 
tween the n tn step and the (n + l) m step of the walk. The element which has 
the largest normalized error in comparing the and the (n + 1)^ steps is 
selected as the criterion function, and this maximum-error-element is allowed 
to vary as the walk progresses. The walk continues until the criterion function is 
minimized, at which point the optimum perturbation step size is considered to be 
determined to within one-half order-of-magnitude. The process is repeated for 
each independent variable perturbation step size, to arrive at an optimum set of 
step sizes, which are then input to the program's iterator in place of the original 
values. A summary of the step size optimization is printed. 

2. Avoiding Corners in the Propulsion-Time Function. The pro- 
gram monitors the thrust switching function along each trajectory to determine 
when to optimally switch the thrust on or off (at the roots of the switch function). 

At times the program’s iterator moves the trajectory into a region where two 
successive switch function roots approach each other and become very close to- 
gether, which physically represents very rapid thrust-switching, and which mathe- 
matically represents a corner in the propulsion-time function. This means that 
the iterator is attempting to either add or eliminate a coast phase or a thrust 
phase to the trajectory. The iterator has severe difficulty when switch function 
roots become very close together, simply because this represents a crease in the 
otherwise-smooth dependent variable functions of the boundary value problem, and 
the iterator is not designed to handle such situations. To avoid this difficulty, a 
program option controlled by the input variables GAP and NHUNG is available 
which monitors the closeness of the switch function roots. When the program de- 
tects that the iterator is having difficulty with a propulsion-time corner, the pro- 
gram halts the iterator, returns logic control to the MAIN program, and intro- 
duces two cases (in addition to the user's input cases). The first additional case 
avoids the propulsion-time corner by forcing the spacecraft to thrust continuously 
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throughout the mission* After the iterator attains convergence with forced-thrust- 
ing, the program attempts to converge on the original case using the forced-thrust- 
ing trajectory as the initial guess* The program generates the forced-thrusting 
trajectory by setting X ^ = 10 X (t^) . 

3. Lagrange Multiplier Scaling . The Lagrange multipliers, or 

o 

adjoint variables A and A , at the start of each trajectory segment may be 
scaled such that the initial mass-ratio Lagrange multiplier has unit magnitude* 

The program will do this automatically if the input variable NORMAL is set 
equal to one. 


4. Generation of an ASTEA Tape. The HILTOP program is 
capable of generating a trajectory tape or sequence of punched cards suitable 
for input to the ASTEA (Arbitrary Space Trajectory Error Analysis) computer 
program^ 5 * * * * * 11 This is accomplished by means of the NAMELIST input variables 
MPUNCH and NTAPE* The contents of the HILTOP output tape or punched 
cards are described in the section, Program Output. 

5. Scratch Pack Output. The HILTOP program was developed 

for use at the computing facility at the NASA Goddard Space Flight Center. That 
facility supports conversational remote input terminals, and special provisions 

have been made in the program for effective use with the remote terminals* 
Specifically, fortran WRITE statements are included throughout the program to 

direct selected program output to unit 11 and unit 12 output devices. 

Job control cards are used to define these units to be special disk packs from 
which the data can be retrieved and printed at the remote terminals. If this 
feature is not wanted or is incompatible at the installation, IBM JCL DD DUMMY 
cards for each unit may be included in the job control cards. The run summary 

is written on unit 12, and the iterator independent and dependent parameter values 
are written on unit 11. The independent variables are written in a format com- 
patible with NAMELIST input and therefore may be used directly in a continuation 

run. 
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6. Normal Run Termination, The computing facility at NASA 
GSFC requires as input an estimate of the machine time necessary to execute 
the run. The operating system will not permit this time to be exceeded, and 
it automatically terminates the run if the input time estimate is reached. To 
avoid losing valuable information should this occur, the time remaining for 
the run is continually monitored during execution (using the GSFC library 
routine REMTIM). If the remaining time becomes less than an input number of 
seconds, control in the program is transferred to MAIN which causes the most 
recent trajectory to be integrated and all requested information printed for that 
trajectory. The run then terminates normally. REMTIM is called only from 
subroutine TIKTOK, and may be dummied by setting its two arguments to 100. 

7. Ballistic Trajectory Option . HILTOP specifically contains the 
ability to generate one- or two-impulse ballistic interplanetary transfer tra- 
jectories, i.e. , trajectories that do not use an electric propulsion system. The 
impulse at launch is provided by the indicated launch vehicle, while the impulse 
at the primary target, if required, is accomplished with a retro stage defined 
by input characteristics. Two options are available for generating ballistic tra- 
jectories. 

The first option is simply to use the program in the same manner as 
for electric propulsion trajectories but force the thrust switch function to always 
be negative such that the electric propulsion engine never turns on. This is most 
easily accomplished by setting the input parameter IBAL equal to 1. This simply 
causes the program to set a number of parameters internally to assure that a 
ballistic trajectory is generated. The program uses the coast phase solution to 
generate the trajectories and the iterator is employed to converge on the desired 
end conditions. 

The second option is more straightforward and automatic. Setting the in- 
put parameter MOPT equal to 1 causes the program to generate a ballistic 
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flyby transfer from the launch planet to the specified primary target; setting 
MOPT equal to 2 results in the generation of a ballistic rendezvous trajectory 
to the specified primary target. In this option the ballistic trajectory is com- 
puted before entering the logic in which electric propulsion trajectories are 
normally evaluated. Therefore, this procedure may be used to obtain the im- 
pulsive solution as a first guess to the electric propulsion trajectory, since 
this option automatically generates values for the initial adjoint variables 

and A and also for v 

o “O 

The latter option requires fewer inputs and less detail to obtain the 

ballistic trajectory, but is restricted to fixed launch date, fixed arrival date 

single-target cases. The former option may be used for multiple -target 

missions and permits one to optimize the launch and/or arrival dates because 

the appropriate transversality conditions are valid for ballistic as well as 

electric propulsion trajectories. It is sometimes very convenient to use the 

two options sequentially, using the second option in the first input case to get 
© 

A^, 'A q and v , and then using the first option in the second input case to 
get the full computational and printout advantages of the program. 

8. Rotation of Primer Vector with Launch Date. The mission 
analyst is frequently faced with a problem of possessing a solution for one launch 
date and desiring a solution for a different launch date< Theoretically, such a 
problem may be solved by treating launch date as both an independent parameter 
and a dependent parameter. However, unless the launch date difference is re- 
latively small, this approach to the solution often leads to numerical difficulties 
due to high sensitivities of the trajectory to arbitrary changes in the components 
of the initial primer vector and its time derivative. 

Generally, the initial primer vector and its derivative will maintain 
similar relationships to the initial position vector as launch date is varied. 
Consequently, a simple rotation of the vectors through an angle equal to the 
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displacement of two initial position vectors will often yield acceptable first guesses 

of the multipliers from which convergence can be achieved. An option is provided, 

using the input parameter IROT, which automatically accomplishes the rotation of 
© 

and A^. In using this option, a case must first be executed to initialize the 
position vector of the original case. This initialization may be accomplished by 
running a single trajectory or a regular iteration. In either case, the initial 
position of the last trajectory generated is used as the reference position. On the 
next case, the parameter IROT should be set to a non-zero value and the initial 
position vector is evaluated using the (supposedly) new launch date. The angle be- 
tween this new position vector and the reference position is evaluated, and A 

• O 

and A from the preceding case are rotated through this angle about the ecliptic 
North Pole. The rotated vectors are then used as initial guesses for the second 
case. 

9. Limitation of Power to Power Processors. An input variable, 
GAMMAX, provides an option of simulating the limitation of a maximum power 
that can be processed by the power conditioners. For example, the propulsion 
system may conceivably be designed for a maximum power equal to that developed 
by the arrays at 1 AU. If a trajectory passed through a perihelion distance less 
than 1 AU, the power developed by the arrays might exceed the design point. 

Using this option, the program would simulate the constraint by assuming the 
arrays are tilted to the sun line by an amount necessary to reduce the power factor 
y to a value equal to GAMMAX. The solar distance below which the arrays must 
be tilted is evaluated internally and depends on the coefficients of the solar power 
law. This distance is evaluated once for each case, and the array tilting is 
triggered on crossing this solar distance. The option may be used with 
the power degradation option, but the tilt angle does not take q < 1 into account. 
To invoke this option, it is necessary to set MODE = 5. 

‘ 10. Housekeeping Power. This program option applies to solar 
electric propulsion and is currently not combined with the power degradation 
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option. The program therefore does not allow the use of both options simulta- 
neously. The program input quantity controlling the housekeeping power 
simulation is DPOW, which is the ratio of housekeeping power to reference power 
and is given the symbol Ap: 

Ap = P h /P ref • (227) 

This spacecraft model assumes that the power required by the operating 
components at each instant of time exactly match the power developed by the 
solar arrays: 

P a = P + P h ’ (228) 

where p is the power developed by the arrays, p is the instantaneous power 

3 . 

delivered to the power conditioners, and p^ is the housekeeping power, which 

is constant with time. Currently, all trajectories generated by the program 

using the housekeeping power simulation must satisfy the condition that p > p , 

a h 

so that p remains positive; in other words, engine shutdown when p “* 0 is not 
coded into the program. 

11. Imposed Coast Phases. It is occasionally desirable to impose 
coast phases over specified time intervals of a trajectory. The option is provided 
through the two input arrays TOFF and TCOAST to impose up to 20 coast phases 
throughout a trajectory. To invoke this option, the elements of TOFF must con- 
tain an ascending sequence of times, in days from launch, at which the engines 
are shutdown. Corresponding elements of TCOAST define the duration, in days, 
of the associated coast phases. To work properly, one should assure that 

TCOAST (I) < TOFF (1+1) - TOFF (I). 

An imposed coast phase may partially or completely overlap an optimum coast 
phase and vice-versa. 
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PROGRAM DESCRIPTION 


This section provides a users' manual of the HILTOP program and pro- 
vides an introduction to the basic logic flow of the program. A brief description 
of the function of each routine is provided in a subroutine glossary, along with a 
subroutine calling sequence table. Subroutine and labelled common cross refer- 
ence tables are presented, followed bj ? tables of each labelled common variable 
cross referenced to every subroutine referencing the variable. A complete list 
and description of the input parameters are given and are followed by a description 
of the computer program output. 

The job control language required to execute an object program module 
residing on a user disc pack at the GSFC IBM 360, Model 91 computer facility 
is presented. Additionally, the basic machine requirements for program execu- 
tion on the GSFC 360/91 system are presented. 


AON SAlvia 


OHiMDUUd 
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A. PROGRAM STRUCTURE 

1. MAIN Program. The subprogram MAIN represents the first 
execution entry-point and the primary driver of the HILTOP program. MAIN calls 
the routines which read the program inputs and perform initialization computations 
for each case, transfers control to the iterator, and calls the summary print rou- 
tines after the iteration phase is completed. A logical flow chart of MAIN is pre- 
sented below. 



BEGIN 


INPUT 


QSTART 



GUESS 


MINMX3 


FINISH 
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2. Subroutine Glossary. 


AEINWT 

ALBEDO 

ANSTEP 

BEGIN 

BOOSTR 

CARKEP 

CDERIV 

CHECK 

CHECKI 

CHKINT 

CONVER 

CONVRT 

CORNER 

DATE 1 
DECLIN 
DERIV 


Assigns target body oi*bit elements for selected comets and asteroids. 

Computes photometric magnitude of a given astronomical body as 
seen by both the spacecraft, and an observer on Earth. 

Takes analytic coast-phase step in state and adjoint variables. 

Assigns default values prior to reading input. 

Initializes launch vehicles coefficients. Entry point in OMASS. 

Converts from Cartesian position and velocity vectors to Keplerian 
orbital elements. 

Computes functions monitored by subroutine CHECK. 

Monitors trajectory for engine switch points, target times and 
other extremum-table entries. 

Initializes CHECK for each trajectory. Entry point in CHECK. 

Initializes CHECK for each trajectory segment. Entry point in 
CHECK. 

Converts vector from ecliptic to equatorial coordinate system. 

Converts derivatives between time and generalized universal 
anomaly. 

Initializes quantities for additional cases which avoid propulsion- 
time corners. 

Computes Julian Date from calendar date. 

Computes launch asymptote declination. 

Evaluates derivatives of state and adjoint variables for numerical 
integration. 
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EFM 

EFMPRT 

ETA 

ETAINT 

EXTAB 

FINISH 

FUNCT 
GET I 
GET Q 
GET RV 

GUESS 

GUNTHR 

IMPRNT 

IMPULS 

INCOND 

INPUT 

INTERP 


Computes target position, velocity and acceleration using stored 
or input osculating elements. 

Prints target and spacecraft position and velocity at target-inter- 
cept times. 

Computes engine efficiency and first derivative with respect to 
jet exhaust speed. 

Initializes constant coefficients for engine efficiency. Entry point 
in ETA. 

Prints extremum table of selected functions. 

Recomputes final trajectory for full program printout, case 
summary, optional punched cards and trajectory tape generation. 

Computes auxiliary functions required at each computation step. 

Computes launch parking orbit inclination. 

Computes iterator dependent variables. 

Computes optimum final position and velocity for extra-ecliptic 
missions. 

Main control subroutine for iteration to find ballistic two-body 
trajectory and associated adjoint variables between specified 
end conditions in given flight time. 

Computes minimum impulsive speed between circular orbit and 
hyperbolic orbit and derivatives with respect to v^ and i . 

Prints initial and final positions and velocities each iteration in 
GUESS. 

Trajectory generator used by GUESS. 

Transforms position and velocity from polar to Cartesian coordinates. 
Reads and writes program inputs (NAMELIST) for each case. 
Interpolates for roots of functions monitored by CHECK. 
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LOAD 

MAIN 

MINMX3 

MORE 

OMASS 

PAR INC 

PDATE 

PMPINT 

PMPRNT 

PRINT 

PRINTR 

PRIOR 

PUNCH 

QPRINT 

QSTART 

RADAR 

READER 

RETINJ 

RIDGE 


Computes quantities for extremum table. 

Program entry and master control. 

Generalized iteration and parameter optimization routine. 

Controls ballistic swingby continuation analysis and computes 
ballistic trajectory extension beyond primary target. 

Computes for a specified launch vehicle, the initial spacecraft 

mass and derivatives with respect to v and i. 

00 

Adjusts perturbation step sizes to maximize accuracy of partial 
derivative matrix. 

Computes year, month, day, hour from Julian date. 

Establishes dimensions for partial derivative matrix print routine. 

Partial derivative matrix print routine. Entry point in PMPINT. 

Prints iterator independent and dependent variables on scratch 
packs (units 11 and 12). 

Iterator summary print routine. 

Saves integrated variables and restores their derivatives at each 
computation step. 

Punched card and magnetic tape output routine. 

Case setup and summary print routine. 

Performs initialization computations for each case. 

Evaluates communication distance and angle. 

Reads special A-format input cards containing independent para- 
meters. Entry point in PUNCH. 

Computes retro engine or electric propulsion spiral performance 
requirements at primary target. 

Monitors propulsion-time corner proximity. Entry point in 
CORNER. 
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RKSTEP 

SCOMP 

SETUP 

SIMEQ 

SMQENT 

SOLAR 

SOLINT 

SPRINT 

STEP 

STORE 

STOREI 

SUMMRY 

SWING 

SWPRNT 

SWSTO 

SWTRAJ 

TAP 

TAPSET 

TFORM 


Takes thrust-phase step (fourth-order Runge-Kutta integration). 

Series computation for f and g series solution of trajectory 
during coast phases. 

Computes iterator logical variables, primarily for use in GET Q. 

Solves a system of simultaneous linear equations. Entry point in 
SMQINT. 

Initializes for SIMEQ. 

Computes spacecraft power ratio yq and related parameters. 

Initializes for SOLAR. Computes solar power law critical radii 
by iteration when coefficients are input. Entry point in SOLAR. 

Point by point trajectory print routine. 

Performs a computation step along a trajectory. 

Reorders and removes multiple entries in extremum table. 

Initializes for STORE. Entry point in STORE. 

Stores specific quantities and prints run summary. 

Computes ballistic swingby continuation trajectory to specified 
post-swingby target. 

Dummy print routine passed to MENMX3 through its argument list. 
Entry point in SWTRAJ . 

Stores output data at switch points. Monitors Hamiltonian accuracy. 
Trajectory generator used by SWING. 

Trajectory integration supervisor. Computes one trajectory seg- 
ment. 

Initializes parameters which force TAP to generate a ballistic tra- 
jectory segment. 

Vector transformation routine. Entry point in VSCAL. 
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THANG 

THANGD 

TIKTOK 

TRAJ 

TRAJI 

TRAVEL 

TRJINT 

TWINKL 

UNITD 

VADD 

VCR OSS 

VDOT 

VMAG 

VPRINT 

VSCAL 

VSUB 


Computes thrust unit vector in fixed cone angle case. Entry point 
in THANGD. 

Computes time derivative of thrust unit vector in fixed cone angle 
case. 

Monitors remaining execution time on computer, by calling 
REMTIM. Provides normal run termination. 

Basic mapping routine for the MINMX3 iterator. Supervises 
trajectory-segment computation and initializes for GET Q. 

Initialization for the beginning of each trajectory. 

Computes incremental travel angle. 

Initialization for TRAJ. Entry point in TRAJ. 

Computes unit vector in direction of Canopus. 

Computes time derivative of a unit vector along an arbitrary 
vector, given the arbitrary vector and its time derivative. Entry 
point in VSCAL. 

Vector addition routine. Entry point in VSCAL. 

Vector cross .product routine. Entry point in VSCAL. 

Vector dot product routine. Entry point in VMAG. 

Vector magnitude routine. 

Print routine for optimum COV Earth departure (e.g. , NERVA). 
Computes product of a scalar and a vector. 

Vector subtraction routine. Entry point in VSCAL. 


Detailed descriptions of all subroutines are presented at the back of this 
document. 
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3. Subroutine Calling Sequence . The table on the following two 
pages displays the hierarchy of calls to the various subroutines in the program. 
The order of subroutine names shown is approximately the order in which the 
calls occur in the listings, but not necessarily the order they occur in the logic 
flow. Multiple calls to any given routine from another routine are noted only 
once in the table. The table on the third page following presents the same infor- 
mation in a different format. Each subroutine is listed alphabetically and is 
followed by a list of all sub-programs that reference the subroutine. 


As noted in the sequence table following, the program contains a number 
of secondary entry points within selected sub-programs. A complete list of these 
entry points follows. 


Entry Point 

Sub-program 

CHECKI 

CHECK 

CHKINT 

CHECK 

RIDGE 

CORNER 

ETAINT 

ETA 

BOOSTR 

OMASS. 

PMPRNT 

PMPINT 

READER 

PUNCH 

SIMEQ 

SMQINT 

SOLINT 

SOLAR 

STORE I 

STORE 

SWPRNT 

SWTRAJ 

THANG 

THANGD 

TRJINT 

TRAJ 

VDOT 

VMAG 

VCR OSS 

VSCAL 

VADD 

VSCAL 

VSUB 

VSCAL 

UNITD 

VSCAL 

TFORM 

VSCAL 


I 

I 

I 


I 

I 

I 
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I 

I 


MAIN 


r-TIKTOK- 


- RRMTIM 
h FINISH — 
*- SUMMRY 


h BEGIN 
QSTART— | 


DATE1 

- AEINWT- 

- SO LI NT* 

- EFM 

- Q PRINT — 

- SETUP 

- BOOSTR* 

- ETAINT* 

- TWINKL 

- TRJINT* - 


-INPUT 

-READER* 
[-CORNER 
rGUESS — 


-TKTOK- 


MINMX3 


PD5 


MINMX3 ■ 
IMPULS - 
f- SMQINT 
SIMEQ* 
SMQINT 


FINISH - 


TRAJ 

- PRINTR 

- PRINT 

- EXTAB 

Q PRINT - 
I- SUMMRY 
h CORNER 


( SEE BOTTOM OF PAGE) 


DATE1 


(SEE CALL BY FINISH) 


r ETA 

OMASS 

RETINJ 

l EFM 


DECLIN 
GUNTHR 
L GET I — 


GUNTHR 


(SEE TOP OF PAGE) 


(SEE CALL BY MAIN) 


-i-SCOMP 

1-ANSTEPjPRIOR 

l-SCOMP 


- PMPINT 


-TIKTOK 

— (SEE TOP OF PAGE) 

f - TRAJ 

(SEE NEXT PAGE) 

- SWTRAJ 

(SEE CALL BY SWING BELOW) 

[- IMPULS 

— (SEE CALL BY GUESS ABOVE) 

[-PIUNTR 


j- SWPRNT 


(- IMPRNT 


- PAR INC 

TRAJ (SEE NEXT PAGE) 

- PMPRNT* 


- SIMEQ* 



(SEE NEXT PAGE) 


EFMPRT-PDATE 

(-BOOSTR* 

(-PDATE 

SOLAR 


PUNCH - 


LMORE • 


EFM 


SWING — rAEINWT — DATE1 
|- EFM 

MINMX3 — (SEE CALL BY MAIN) 
SWTKAJ— j-EFM 
C ON VER -TAPSET 

TAPSET L tap (SEE CALL BY TRAJ) 

(-TAP (SEE CALL BY TRAJ) 

EXTAB 

L EFMPRT- PDAT E 


♦Entrj' Point 
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TRAJ 


1-TAP 


hEFM 


rTRAJI 


1-ETA 


-RIDGE* 
-EFM 
-RADAR — 
CHECKI*- 
SOLAR 
-THANG* 
THANGD 
SOLAR 
CIIKINT*- 


-EFM 

■STOREI* 


•CDERIV 


• LOAD - 
STORE 


THANGD 
RADAR — r-EFM 
LAE 

TRAVEL 


AEINWT 


CONVRT 
DERIV — 
CDERIV - 
FUNCT— 
STEP 


-SOLAR 
-(SEE ABOVE) 

-(SEE CALL FROM STEP) 
-RKSTEP r-PRIOR 

L] 


T 


[-CHECK- 


- DERIV SOLAR 

PRIOR 
SCOMP 

LFUNCT p-THANGD 

-THANG* 

-SOLAR 
-(SEE ABOVE) 

STEP — (SEE ABOVE) 

ODER IV (SEE ABOVE) 

TRAVEL 


STOREI* 

1-TRAVEL 1-SOLAR 

FUNCT— 


-ODER IV 
■INTERP 




-LOAD - 
■ STORE 
-STEP — 


(SEE ABOVE) 


SPRINT ■ 


-THANGD 


LSWSTO — 
■ OMASS -yDECLIN 

Lgunthr 
Lget i — 

RETINJ 

•GET RV— INCOND 
GET Q 
V PRINT 
PRINT 
RIDGE* 

PRINTR 


P THANG* 

-ALBEDO LsOLAR 

■SOLAR 

■CARKEP 

•EFM 

-PUNCH EFM 

• SOLAR 


GUNTHR 


♦Entry Point 
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SUBROUTINE CROSS REFERENCE TABLE 


NAME 
AEINWT 
ALBEOO 
ANSTEP 
BEGI N 
BOOS TR 
CARK EP 
CCHECK 
CDER IV 
CHECK 
CFECK I 
CHKI NT 
CONVFR 
CONVPT 
C OPN ER 
DATE 1 
DECL IN 
DER I V 
FFM 

EFMPRT 

ETA 

ETAI NT 

EXTAB 

FINISH 

FUNCT 

GET I 

GETG 

GETRV 

GUESS 

GUNTHR 

IMPULS 

I NCOND 

INPUT 

INTERP 

LOAD 

MI NMX3 

MORE 

OMASS' 

PAR I NC 

PDAT E 

PD5 

FMPI NT 
PMPRNT 
PRINT 
PR I N TR 
PRIOR 
PUNCH 
QPRI NT 
C START 
RACAR 


OSTAR 

SPR IN 

IMPEL 

MAIN 

GPP IN 

SPR IN 

MINMX 

CHECK 

TAP 

TRA J! 

TAP 

SW ING 

TAP 

FINIS 

AF INW 

CM ASS 

RK STEP 

ALBEDO 

TRAJ 

MCRE 

TRAJ 

OSTART 

FINISH 

MAIN 

CHECK 

OMAS5 

TRAJ 

TRAJ 

MAIN 

GETI 

GUESS 

GETRV 

MAIN 

CHECK 

CHECK 

GUESS 

OPR INT 

TRAJ 

MINMX3 

EFMPRT 

M INMX3 

MINMX3 

MINMX3 

FIN ISH 

FINISH 

ANSTEP 

FINISH 

FINISH 

MAIN 

CDER T V 


RA CAR 
STEP 
QS T ART 


INTERP 


MA I N 
QSTART 
PR INT 
TAP 
PUNCH 
TRAJ I 
QPRI NT 


MORE 
T I K T OK 
STEP 


CMASS 


MAIN 


QPRINT 


SOLAR 
SOLAR 
RK STEP 
SPR I NT 
QSTART 

TRAJ I 


ubrcutines 

T 
T 

S 

T 
T 
3 


H 

T 


EFERENCING MEMBER 
SW ING 


TAP 


QSTART RADAR SPRINT 


TAP 


SWING 


TR AJ 
TRAJ 


SWI NG 


S WT RAJ 
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CROSS REFERENCE TABLE (CONTINUED) 


NA ME 

SUBROUTINES REFERENCING MEMBER 

READER 

MAIN 





REMT IM 

T I KTOK 





RET I NJ 

TRAJ 





RIDGE 

TRAJ 

TRAJ I 




RKSTEP 

STEP 





SCOMP 

ANSTEP 

I M PULS 




SETUP 

C S T ART 





£ IMEG 

GUESS 

MINMX3 




SMQ1 NT 

GUESS 

MINMX3 




SOLAR 

CHECK 

TRAJI 

OERI V 

FJNCT 

Q PR I NT 

SPR I NT 

SOLI NT 

GST ART 





5PRI NT 

INTEPP 

TAP 




STEP 

CHECK 

INTERP 

TAP 



STOP E 

CHECK 





S TOR El 

CHECK 





SUMMFY 

FINISH 

T I KTOK 




SWING 

MORE 





SWSTC 

TAP 





S WTP A J 

SWING 





TAP 

VCPE 

SWTR AJ 

TR AJ 



TAPSET 

MOPE 

SWTRAJ 




TFORM 

SPRINT 





THANG 

FUNCT 

TAP 




THAN GD 

CDERT V 

FUNCT 

TAP 



TI KTCK 

INPUT 

MA I N 

MI NMX3 



TRAJ 

FTN ISH 

PARING 




TRAJ I 

TRAJ 





TRAV EL 

LOAD 

TAP 




TRJI NT 

GST ART 





TWIN Kl. 

GST ART 





UNIT D 

CCER I V 

THAN GD 




VADD 

CDERIV 

TAP 




VCRO SS 

CARKEP 

CDERIV 

getrv 

INCCND 

CPR I NT 

VOOT 

ALBEDO 

CARKEP 

CDERIV 

RADAR 

SPRINT 

VM AG 

ALBEDO 

TAP 

CARKEP 

CDER I V 

QPR INT 

RADAR 

VPR I NT 

TRAJ 





VSCAL 

CDERIV 

INCOND 

SPRI NT 

SWING 

TAP 

VSUB 

ALBEDO 

CDERIV 

GETRV 

SWING 



SWSTG 


SPRINT 
SWING 
SPR INT 


TAP 


SWING 

SWING 
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4.' Common Array Information. Throughout the HILTOP program 
a total of 8 labelled common arrays are employed. In the tables to follow are 
presented cross-reference information sufficient to inform the user as to the 
occurrence of specific common arrays and of specific common variables through- 
out the program. The first table presents for each common array a list of all 
subroutines in which the named common appears. This is followed by separate 
tables for each common, containing subroutine references to each common vari- 
able. These tables list only variables that are actually referenced one or more 
times throughout the program. For each variable are listed the Fortran name, 
the type of variable (i.e., R*8 for double precision real, 1*4 for integer, etc.), 
the address of the variable in decimal bytes relative to the start of the common, 
and the name of the subroutine in which the variable appears. A separate line 
appears for each subroutine in which the variable is referenced. For the second 
and subsequent subroutine references, no other information is repeated unless 
the Fortran name of the variable is different, in which case the new name and the 
relative address are repeated. The definition of any referenced common variable 
is given in the External Variables Table of the subroutine referencing the variable. 
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COMMON CROSS REFERENCE TAELE 
NAME SUBROUTINES F.'EFERENCING MEMBER 


•EXTREM 

BEGIN 

COER IV 

CHECK 

EXTAE 

INTERP 

LOAD 

PRINT 


STORE 

TAP 






GUNC CM 

GET I 

GET Q 

GUNTHR 

C MA S S 




INTG F<4 

BEGIN 

ODER IV 

CHECK 

C CRNEP 

DER IV 

EFMPRT 

EXTAB 


FIN ISH 

GETQ 

INPUT 

I NT ERF 

LOAD 

MAIN 

MINMX3 


MORE 

PARING 

PR INT 

PRINTR 

PRIOR 

PUNCH 

GPR INT 


QST ART 

RADAR 

RK STEP 

SETUP 

SOLAR 

SPRINT 

STORP 


SUM MR Y 

SWING 

SWSTO 

S wT R A J 

TAP 

TAPSET 

TIKTOK 


TR AJ 

TRAJ I 






I TER AT 

BEGIN 

CDER I V 

CORNER 

FINISH 

GET G 

GETRV 

GUESS 


INPUT 

MI NMX3 

□ MASS 

PARI NC 

PR I NT 

PR I NTP. 

PUNCH 


QPRINT 

QST ART 

RADAR 

SETUP 

SOLAR 

SPR I NT 

SWING 


TAP 

TRAJ 

TRAJ I 

VPR I N T 




I TER 2 

BEG IN 

GETQ 

GUESS 

I MPULS 

MI NMX3 

PAR INC 

PRINTR 


QST ART 

SWING 

SWTRAJ 

TIKTCK 

TRAJI 



LOGI C4 

BEGIN 

CDER IV 

CHECK 

CORNER 

DECLIN 

DER I V 

EXTAB 


FINISH 

FUNCT 

GETI 

GETC 

GETRV 

INTERP 

LOAD 


MA I N 

MI NMX3 

MORE 

C MAS S 

PAR INC 

PR I NT 

PRINTR 


PRIOR 

PUNCH 

QPRINT 

Q STAR T 

RADAR 

RET INJ 

RK5TEP 


SETUP 

SOLAR 

SPRINT 

STEP 

STORE 

S'JMMRY 

SWING 


SWSTO 

TAP 

TAPSET 

TFANGD 

TIKTCK 

TRA J 

TRAJI 

real e 

ALBEDO 

AN STEP 

BEGIN 

CCER IV 

CHECK 

CONVRT 

CORNER 


DECLIN 

DERI V 

EFM 

EFMPRT 

ETA 

EX T AB 

F IN I SH 


FUNCT 

GETI 

GETQ 

GETR V 

GUESS 

I MPRNT 

I MPULS 


INPUT 

interp 

LOAD 

MAIN 

MINMX3 

MORE 

C MASS 


PRINT 

PR I NTP 

PR I OR 

PUNCH 

QPRINT 

QST ART 

RADAR 


RFT INJ 

pkstep 

SOLAR 

S PR I N T 

STORE 

SUMMRY 

SWING 


SWSTO 

SWTRAJ 

TAP 

TAPSET 

THANGD 

TI KTOK 

TRAJ 


TRAJI 

TRAVEL 

TW INKL 

V PR I N T 




SOLS YS 

BEGIN 

EFM 

EFMPRT 

M CR E 

QPRINT 

QST A RT 

SPR INT 


SWING 

TRAJ 
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T A 3 L 


COMMON VARI /• 31. £r CROSS Kt'FRcNCc 


SUMMON IThR s 


LONG TH 


Vv.Y,~, 


OF POOR QUALITY 


var I a:-;l t 

i yn 

ADDS 

8 : J.3 R our 1 N 

i 

R A 8 


G UE S 3 
1 N G 
TRA JI 
I I'VUlS 
M 1 NM X 3 
PA PI NS 
0 5 T A o r 
S a TP A J 
T I KTO< 

p i 

r< r ' -15 

r 

3 EG I M 

0 

R 6 

?P.C: 

Gk i 3 
5 W I •' v G 
1 M PUL 5 
M I NM X 3 
PA PI NS 
5 V. TP A J 

ns 

R * » 

56 C 

G UE 5 S 
S \v 1 M G 
MI NM X 3 
OS TART 

JW 

R><fc 

84 r 

G US 5 5 
5 V. I .N G 
MI NM X 3 
QG TART 

0 M 1 N 

w * a 

1 \ 2 0 

6UES 5 
S Wl N 3 
M I NM X 3 
PR INTR 
05 TART 

G M A X 

R <r 8 

3 408 

G UE 3 f. 
S V. I N G 
M I NM X 3 
PR I NT R 

o j i a r r 

OB 

R -r a 

1 680 

PAP I NS 

BBS 

R * a 

1 680 

GUESS 
S W I N G 
M I NM X ?. 
OS TART 

YYY 

a * a 

3. 960 

M l NM X 3 

8 NO MX 

r * a 

2 24 0 

MI NM X 3 
T I KT □ < 

R M 

r a a 

2 520 

G UE S S 
M I NM X 3 
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C OF M ON VA9 I A BL £ 


RE^LQtNfCE T A 3 L 


CROSS 


COMMON REALS 


L ENj T H IS DO 


VARIABLE TYPE 

R ? I P-? 8 

SAVE R*S 

PAYLOD R*6 


AM p»g 

X MAS S P A a 


CTANK P 8 



FT R * 8 



VJ R«8 


VIMP R*8 


ADDR 

0 


3 

p 


7 £ 


es 


96 


104 


S UBR cu T I N 
BEGIN 
SU M.M.P Y 
GE TO 
TP A J 
P P I N T 
PUNCH 
PR INTP 

op pi nt 

DMA 3 S 

geto 

TP A J 
E XT A 3 
PRINT 
PUNCH 
PP INT P 
OP PI NT 
RETIN J 
TAPSET 
TP A J 
3 EG I N 
I NPJT 
PUNCH 
QP PINT 
TP A J 
I NPUT 
P UNC H 
OPRI NT 
T A D 
GET 3 
OERIV 
F UNC r 
S OL A P 
SWSTD 
T PA J I 
CD EP 1 V 
QSTAPT 
RETINJ 
SPRI NT 
GET 3 
TRA J 
SQL A 3 
TRA JI 
CDERI V 
QSTAPT 
RETI NJ 
VPRI NT 
TRA J 
TRA JI 
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C 0 M M (j N «EAL2 ( C t J f J 7 I N U E D ) 


VAR I A3L -I 

TYPE 

A DDR 

S U ft ft OUT I 
DE CL I M 
OP PINT 
OS TA-7T 

OT:V AS". 

P-'8 

1 3. 2 

TR A J 
P UNC H 
GPP.IMT 
RE TINJ 

L>M P £ T ft 

P- -j 

152 

I NPUT 
P UNC 8 
RE TI NJ 

CTfttT.T 

ft * S 

1 6C 

TR A J 
-3 EG I N 
I iNPJT 
P UNC H 

ft PER 

ft * 8 

1 6 8 

U c G l N 
I NPUT 
PUNCH 
CP PI NT 
RE TINJ 

PAP 

p * y 

1 76 

3 EG I N 
I NPUT 
P UNC H 
OP ft] NT 
RETI NJ 

THRET 

k >- a 

IS* 

8 EG I N 
I NPUT 
PUNCH 
OP ft I IN T 
RE TINJ 

SP I RET 

R * 8 

i 02 

3 EG I N 
I NP JT 
P UN C H 
OPRI NT 
PE TINJ 

DVEL 

ft tS 

200 

P UN C H 
OPRI NT 
RETI NJ 

V L OSS 

ft * 3 

20 Fs 

TR A J 
PUNCH 
OPRI NT 
RE TINJ 

V CP B 

ft 4- 8 

216 

PUNCH 
QPPI NT 
RE TI NJ 

ASCL 

Rj- 8 

22 A 

I NPJT 
SOLAR 

XJLD 

Pr- 6 

2 6 A 

OPRI NT 

TR I P 

R * 8 

272 

D UNCH 
(PRI NT 

PT 

ft * 8 

288 

P UNCH 
QPPI NT 
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COMMON! 


REALM (CONTINU'D) 



VAR I ABLE TYPE 

F R * 8 

FFF R*-0 

P tU M R»>e 

fivlAX A! A a 

A N G D R v ft 

HOUR R*8 

TPM AX 9 >' ft 

A A I R ~ 6 

STATE rah 


XO R ■* 6 


P C. R * 8 


6 1 RT6 


O I R -t G 


El R* Q 


01 R * 8 


ADDR 

SUBROUTI 

29 6 

P UNC 8 

296 

OP PI NT 

30 4 

PUNCH 


OPR I NT 

312 

PUNCH 
OP PI NT 

32C 

PUNCH 
0 ' J P I N T 

328 

a EG I N 
i nput 

PUNCH 


qstapt 

33 6 

0 ETG i N 

3*4 

1 nput 
PUNCH 


os tar. t 

3S? 

MAIN 
BEGIN 
I NP JT 
P UNCH 


OS TART 

4C 0 

GETO 
TP A J 
BEGIN 
GETRV 
GUESS 
I NPUT 
PUNCH 
RADAR 
TRAJI 


OPRI NT 

ostar r 

456 

GET3 
TRA J 


G UE 5 S 
TRAJI 

512 

B EG I N 
I NPUT 
PUNCH 


OPRI NT 
QSTART 

52 C 

BEGIN 
1 NPUT 
PUNCH 


QPRI NT 
QSTART 

523 

I NPUT 
PUNCH 


Q3RI NT 
QSTART 

53 6 

I NPUT 
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COMMON 

RE ALE- (COW 

T I M .UE D ) 


VAR 1 A Ol. !>' 

TYPE 

A OOP 

EUCiP OUT 1 
R UN C H 
CP PI NT 
0 S T A R T 

8 2 

P •- 3 

54 4 

I NPUT 
P UN C 5 
OP PI NT 
OSTAPT 

:3 3 

• t<*b 

55 2 

I NJRUT 
P UN C H 
GP PINT 
os tap r 

PO A F I X 

r< * s 

5 6 0 

T R A J 
B EG I N 
I N^ JT 
PUNCH 
OSTAPT 

AR 

;P <• £} 

.56 8 

B EG I N 
GETPV 
I IMP J r 
P UN C H 

cs tap r 

AV 

,".^C 

R a c 

5 76 

GETPV 

OSTAPT 


P A 6 

58 4 

G ET P V 
OS TART 

F. I 

R * £ 

552 

G ET P V 
OSTAPT 

C I 

Si 8 

6C 0 

GETS/ 

OSTAPT 

AE 

Pi 8 

608 

I NP U T 
PUNCH 
GSTAPT 

2 A I 

IP X fc 

61 6 

E r M 
BEGIN 
INP'JT 
P UN C H 
OP RI NT 
08 TAP T 

EC l 

8 

62 4 

EFM 
I NP JT 
PUNCH 
OSTAPT 

CM I 

P i; Pj 

632 

E r v, 
I NOJF 
PUNCH 
OSTAPT 

CM I 

P* 3 

64 0 

EFM 
I NPUT 
PUNCH 
OSTAPT 

SO I 

P * 6 

64 a 

EFM 
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COMMON 

REAL?. { CON 

T I NUGC > 


V AR I ABL B 

TYPE 

ADOP 

SU8R OUT 1 

i input 

PUNCH 
G £ T A R r 

TP I 

R ci 

65 6 

t •- Y, 

1 NPJT 
PUNCH 
OP R l NT 
03TART 

EMU ODD 

p xi Aj 

66 A 

E”M 
T P. A J 

i kput 

PUNCH 
S Wl No 
OSTART 

PA DODD 

O t fi 

67 2 

TP A J 
8 EG I N 
1 NPJT 
PUNCH 
. S W I N 3 
OS TART 

0 BET A 

or 

a 

680 

T A 3 
C HE C < 
I NTF..VP 

STEP1 

R * Q 

683 

TA 3 
BEGIN 
I NPJT 
P UNC H 

STEP 2 

R*e , 

696 

TA 3 
BEGIN 
I NPJT 
PUNCH 

T M A X 

r* e 

70 A 

T A 3 
GE T, 0 
TRAJ 
GUESS 
TRAJI 
I MPULS 
QSTART 
SWTRA J 
TAPSEI 

T 2 

R*'fi 

71 2 

MORE 
BEGIN 
I NPUT 

GOO 

R Y 6 

792 

S A I N 3 
S v. TR A J 

TDV 

R*3 

80 C 

TA 3 
GET 3 
BEGIN 
I NPJT 
OSTAST 

TOP.LV 

R*8 

80 8 

TA 3 
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COMM ON 


ft H ALT: ( C OKT I MUc D ) 


VAR 1 A L r 

• VP E 

GAP 

ft*rt 

T G 

R'v 8 

S V: HUNG 

H *- 8 

HO HUNG 

ft ~ 8 

XCOM 

R a 8 

GL A:-. E 

R <: Q 

5 : f; r V. A 

R 6 

c: r m n 

r * 8 

5. t£PMC 

ft* 8 

Sfc'FMO 

P*e 

ANGLE 

R* 8 


T BURN 

e 

DATE’ 

ft * e 


a on ft 

SUftft OUT 1 N 
0; TART 

S 1 6 

BEGIN 
C HE Z < 


I NPUT 
FI Nl 3 H 
OS TART 
T I K T 0 < 

8? 4 

C HEC < 
CORNER 

SO 4 

CHECK. 

CORNER 

904 

C HE C < 
CORNER 

944 

LOAD 
CmEC < 
RADAR 
COER 1 V 

99? 

LOAD 
ft AD A R 

1 ccc 

GET D 
TRA J 
TRAJI 
QPRINI 
OS TAR r 

1 0 58 

GFTD 
TRA J 
T RA JI 
OPR I NT 
GS TART 

1 3 3 a 

GET D 
TRA J 
TRAJI 
OS TART 

3 ?. 6 3 

GETD 
TRA J 
TRAJI 
OS TART 

1 22 4 

TA 3 

LOAD 

PRINT 

TRAJI 


OPR I NT 
SPRI NT 
TAPSET 

1 232 

PUNCH 
QPRI NT 

3. 24 0 

TRAJ 
PUNCH 
T RA J I 
GP R I N T 
GSTAR7 
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COMMON 


PEALS (CONTINUED) 


VAP 1 ABLE 

TYPE 

A DDR 

SUBRCi'JT 1 N 

TO ATE 2 

R*8 

3 24 e 

MO R E 
TRAJ 
3 A I N j 
TRAJI 
OPR I or 
OSTAP r 

ID ATE X 

R* >3 

1256 

EK vi 
TRAJ 
PUN CP 
P AD A 3 
S Wl No 
TRAJI 
OPRI NT 
qstart 


£PRI NT 


r\V 


SKA 

r t a 

1 264 

T W I N< _ 

•SDC 

R*fc 

1 272 

Ttf I NKL 

SPV 

R » a 

128 0 

CPF: I NT 
T ,v 1 N.< _ 

TANGLE 

p -*< e 

1 30 A 

LOAD 
TRAJI 
TAPSET 
TP AVE_ 

A A A 

n*8 

5 31 2 

GETI 

GMASS 

6B0 

r» a 

1 320 

GETI 

OMASS 

ccc 

P *•- 8 

1 328 

0 MA S S 

-A A N G 1 

r ^ a 

1 335 

GETI 
0 MASS 
PUNCH 
T PA J I 
GPP I NT 

AN 02 

R * 8 

1 34 A 

0 MASS 
PUNC H 
OP PI NT 

E 

R « 6 

1 35? 

. BEGIN 

SE 

R v- 8 

1 360 

GET 3 
BEGIN 
TRAJI 
DECL IN 
T «i IN<_ 

CE 

R * S 

1 366 

GETD 
8 EG I si 
TRAJI 
DECL IN 
TW INK. 

OfcCL 

R* 6 

1 376 

GET3 

OMASS 

PRINT 

PUNCH 
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reals (cGNTitwe:oi 


COMMON 
V AR I ALL E. 

A i 3 

X INCL 

A L T A !J 

T Cj 0 
R M A X 

fc V, I M 

P M A X 
COMO I S 
CO ANG 
RTS w I T 

P -5 1 
t ml-:ta 

PHI 


TYP'rl 
P ■* S 

u * a 

R 3 


f< 5 ^. 

R •--- >j 

k 8 

K 6 

r v- a 
Rv 8 

R-S 

R*-e 


AO DR 


1 38 A 



.1 40 8 
1416 

i 4 ? A 


1 4 4 C 


144 8 


I 486 


1464 
1478 
1 48 0 


SUP.R CUT I \! 
Ijr- CL l vi 
OPR l NT 
T A 3 

r uncf 
solar 
DMA ss 

PUNCH 

or pi Nr 

T A 3 
CC T 0 
F UNC r 
S V ST 3 

tra j i 
sp w i n r 

TARGET 
MOHS 
BEG I N 
I NP J T 
C H £ C. < 
PRINT 
P UN C H 
STOR£ 
OPR I NT 
CHEC< 
PRINT 
P UNC H 
STORE 
QPR1 NT 
CHEC< 

S TORE 
OPR I NT 
£ XT AS 
PUNCH 
TARGET 
c. XT A 3 
PUNCH 
TAP3ET 
TAP 
£ XT A 3 
COER I »/ 
OS TART 
.SP RI NT 
LOAD 
COER I V 
SPRINT 
LOAD 
COER I V 
SP R-I NT 
LO A3 
COER I 
SPRI NT 
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COMM CN 


REAL 6 (CONTIWUdO) 


VAR I A9L C- 

type 

A DDR 

SUBR CUT I V 

G PLAN 

ft* 0 

1 469 

(jcTI 
TRA J 
OETIMJ 

PPL AN 

ft =! c 

1 496 

TRA J 
RETINJ 

CP I SPY 

R -/ 8 

I ft 04 

OPRINT 
RE TINS 

TH3PY 

R * 8 

1512 

OPR I NT 
RE TIN J 

T I MSP Y 

R - 8 

1 52 C 

OPRINT 

RETINJ 

X MSP Y 

i< -/ 8 

.1 52 3 

&PRIMT 
ft E T I N J 

fi V C 

ft v ft 

1 53o 

R AD A ft 
CDE.R I \/ 

X ANG 1 

R •• S 

1 584 

I NPUT 
DMA 5 5 

X ANG 2 

p * e 

1 592 

I NPUT 
Cf/ASS 

ft 

p* 8 

3 60 C 

T A 3 
GETO 
LOAD 
TRA J 

kunct 
PUNCH 
CD Eft I V 
CONVRT 
OPR I NT 
SPRINT 

PP 

r >• e 

1 61 6 

TA 3 
C. HEC < 
FUNCT 
COER I V 

5 W ITCH 

P * 8 

1 632 

TA 3 
GET 3 
LOAO 
C HE C < 
FUNC r 
CDERI V 
SPRI NT 

QVt ITCH 

RT 8 

1 64e 

TA 3 
FUNC T 
CDERI ✓ 
VPSINT 

S A I X 

R A 5 

1 800 

EFM 
I NPUT 
RADAR 
S WI NS 
OSTART 

EC I < 

R A 8 

1 640 

EF vl 
I NPUT 
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C CM M (>l 


P EAL8 (CONTINUED) 


V A «? I A B l_ r. T CPE A D D R £ U B R O' J T 1 N 





R ADA 3 
S W I N j 
G S TAP r 

x 

G 

P.y S 

1 0 8 0 

E F v < 
i np-j r 

R ADA ? 

s i*. no 
ostar r 

C M I X 

X <"■ 8 ' 

5 c - ° 0 

E F x, 

i npjt 

RADA? 
b I'l N j 
0 START 

£ Cl I X 

P * 8 

*. 9 6 0 

E = vl 
I NPJT 
f< ADA P 
S rt I N 3 
o. c , tapt 

TP I X 

R *■ 8 

200 c 

Mr VI 
I INPUT 
R AJAR 
SWING 
qsta?t 

MM 1 1 OO X 

R * 3 

EGA A 

E F X 

i np j r 

RADA? 
5 W 1 N 3 
OSTAP, r 

PAD COX 

P» b 

£ one 

I NPor 
R ADA? 
b WI NS 
0 ~> TART 

AL PHA A 

P :Sc 8 

2 ?. 20 

TP A J 
BEG I N 
I NP U T 
P UN C H 
OP PI NT 

ALPHAT 

PA & 

2 1 2 8 

TR A J 
BEG I N 
I NP JT 
PUNCH 
QPP1 NT 

AN 

^ n* & 

2 ACC 

TA? 
BEGIN 
DEP I V 
I WPUT 
PUNCH 
CONVRT 

R «. N 

Rvfi 

240 8 

TAP 
DEP I V 
COER I V 
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COMMON 

PEAL ft 

(CONTI NUF.P) 


V AR I A3L E 

TYPE 

AD DP 

SU3R CUT I ME 
CClNVP.T 

P?N 

R * 8 

£4*6 

T A D 
DERI V 

V DOO 

R * 8 

2 4 2 4 

GE TO 


GETRi/ 
T RA J I 
OE CL i N 
OP PI MT 


P I MO 

P. ft 

5 44 8 

GETO 
TRA J 

P I GN 

r * a 

£4ft6 

GE TO 
TRA J 
BEGIN 
I NPUT 
PRI NT 
TRA JI 

TSC ALE 

PA 6 

2 464 

BEG I N 
I NR J T 
QSTART 

DE CLAM 

R ■% g 

£ 472 

PRINT 

TRAJI 

REVS 

R« 8 

2 4 80 

G UE S S 
I NPUT 

XA.MODA 

R + 6 

2 4 8 8 

PRINT 
MI NM X 3 




PR INTP 

DEG 

R v- ft 

2 4 96 

EFM 


GET I 
GETO 
LOAD 
BEG I M 
O MA E 5 
PRINT 
P UN C N 
S Wl M3 
ALBEDO 

ef mpr t 

QPRI NT 
03TAPT 
SPRINT 
TW INK- 


FP SIN MM 

R * ft 

2 50 4 

GETO 
BEGIN 
QSTART 
VPRI NT 

SUN-MU 

R*8 

2512 

E=M 
TRAJ 
BEGIN 
Q START 

C 0 N T M 

R4 g 

2 52 C 

TAP 

GETO 
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{ C r? M T I NUE D) 


COM WON ft F ALE 

V Aft I A 31. E TYPE ADDR S J C RCUTI N 

LOA5 

T r - A J 

o eg i m 

C HE C P 

p pi n r 

P UNC H 
P. ADA ? 

S QL A 5 
S W I N 3 
5 Vs S T 3! 

T RA J I 

corner 

I N TE ? 3 
OP PI NT 
OSTAP r 
EP R I NT 
S W TP A J 
T A P3 E T 

CON jP R*8 2528 £ T A 

T A 3 
GET 3 
REGIN 
C M A 3 5 
SV.IN3 
E r MP ft r 
OP PI NT 
P F T I N J 
SPRINT 


CGNAO 

P * 3 

2 536 

GET 0 
BEGIN 
P E T I N J 

C.ONPW 

R * H 

254 4 

TP A J 
B EG 1 N 

P I 

P * 8. 

2 55? 

0 EG I N 
R ADA 3 

T v. G* 3 I 

R > - b 

2 560 

Er 'A 
3 EG I N 

CO NOS 

ft *8 

£56 5 

BEGIN 

3PPINI 

C.ONLBO 

R * 6 

2 576 

BEGIN 
PE TI N J 

CONG 

R>' 3 

£ 58 4 

BCG I N 
PE T I N J 

DCLTAV 

R ; i : 8 

£ 59 ? 

TA 3 
GET 3 

XS a ING 

R 6 

260C 

MOP = 
I NPUT 

YE t. ING 

K * 6 

3C-4 0 

MOP S 
5 V. I N 3 

TCIFF 

ft « 8 

4 2 C 0 

JP.AJ 


BEG I N 
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COMMON 


KEAL8 (CONTINUED) 


1 AR I ABLE 

TYPE 

A DDR 

5UBROUTI \i 
I NP J T 
OPR I NT 
GSTAPT 

TCOAST 


A 36 0 

TRA J 
1 NP'JT 
OP PI NT 

T C HECK 

r* e 

A 52 0 

T A 3 
TRA J 
C HEC < 
CDL-.R I Y 
TAPSET 

F 

R ~ £ 

5 5. 6 C 

G UE 3 3 
ANSTE 3 

G 

n « v 

516 8 

AN STEP 

' FX 

R<! b 

517 6 

G UE 3 5 
ANETE 3 

o K 

Ri-a 

5224 

GUESS 

ANSTEP 

ER 

Ft 3 

544 3 

ostapt 

EPD 

tt*8 

5472 

OSTAPT 

FN 

R *• 3 

5 4 96 

Ostap r 

END 

R -••• 8 

5 52 C 

OSTAP r 

AXIS 

R«3 

5 54 4 

DEP I V 
OSTAPT 

A X I f, D 

R A 

5 568 

OSTAPT 

4L T IT U 

R -* 6 

5 592 

GET0 
I NP JT 
QSTARI 
VPRI NT 

T 5 W 

RtQ 

56C0 

5 V ST 3 
PRINT! 

H j W 

RT Q 

6 0 DC 

SVST3 

PVFLUC 

R A 6 

64 00 

E te I N 3 
Sv.'TR A J 

T ii A 3 E 

R*6 

6 4 24 

S tel N3 
S'.v TRA J 

T S ■ J M 

R * eJ 

6432 

MO P E 
S tel N 3 
S W TP A J 

zetate 

R * a 

6440 

SWI N 3 
S w TR A J 

ETH 

R * 8 

6 60 0 

T A 3 

DEP I Y 
FUNCT 
S i’. S T 3 
CD ER I Y 
OSTAPT 
SPRINT 
T H AN 3 3 

L7HD 

R*8 

6624 

T A 3 
F UNCT 
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C. GMM ON R Al. «•; ( C ON T I MUZ U ) 


V AR I ABL “ 

TYP f. 

A 0 0 ft 

SUB ft CjUT 1 
CDtS I 7 
OSTAP r 
sp pi mt 
THAN GO 

S I T 

R ~-8 

6 72G 

T A 

CHECK 
D EP I 7 
F UN C T 
CDE.ft I 7 

S I PH I 

ft ± a 

6 7 2 8 

D 2 ft I V 
CO 2 ft! I 7 
TH AN GO 

Clip HI 

Rr. .3 

L ~7 3 6 

T A A 

D fc‘ ft I 7 
CD HP I 7 
TH AN GO 

PLC 

ft * 0 

£ 744 

TAR 
10 E R 1 7 
F UN C T 
S 'AS r 0 
s p r< i n r 

HAM 

ft ci 

6 760 

TAR 
TP A J 

HA MX 

R . a 

6763 

' GST 0 
TP A J 

F THS W 

ft. * C. 

6 300 

S V- 5 T D 

X JL A 

r<? a 

SCO 0 

OP PI NT 

HC'JRU 

ft * 8 

0 00 ! 3 

OP PI NT 

HOUR A 

R >'• b 

601 6 

CP PI NT 

r 3 

P .£/ 

6 02 A 

OP SI NT 

C A 

K ^ C 

f 032 

OP PI NT 

!;£p 

P r: ri 

6 04 0 

OP ft I N T 

ARP 

iM 8 

804 8 

0 *■> PINT 

ASP 

r ~ a 

f 056 

0 R ft I N V 

DPOW 

ft ? a 

6056 

TkAJ 
I NPUT 
■ s ola ft 
QRRI NT 

GAMMA X 

ft A 0 

6064 

BEGIN 
I NPUT 
SOLA ft 

ONi I N 

R * E 

9C72 

CHEC < 
SOL A ft 

DP DMA X 

R*8 

6 06 0 

TAR 
SOLA ft 

CHFNC2 

Rj C 

90 8 8 

SOLA ft 
CDEPI 7 

CHFNC 

ft v y 

6112 

C HEC < 
SCLAft 
TftAJ I 



C CM M O M P E A L 3 (CONTI N'JC D I 


V A.-; I A3L 2 

TYRE 

A DDR 

S'JDRGUTl 
C D L R I V 

SAMP?, 

R 8 

9180 

TRAJ 
OPR I NT 

DROPS 

R v Q 

o 1 2 a 

TRA J 
OP PINT 

PC DOT 

9 * 8 

9136 

TAP 
F UN or 
CD ER I Y 

PMN 

R* 0 

916 0 

TA3 
GETO 
TP A J 
0 tip I V 
F UNCT 
SOLA? 
T PA J I 

RMS 

• R k- 8 

91 66 

TRA J 
O E R I y/ 

func r 

TRAJ I 

RT 

R*6 

9176 

TAR 
D ER 1 </ 

funct 

SOLAR 
Sv.STD 
TRAJ [ 
AN STEP 
OPRI NT 
SPRINT 

PS 

R * 8 

91 84 

OtRIV 
S W S T D 

PC 

r * a 

919? 

TAP 
D EP I V 
F UNC r 
S A STD 
A N ST E p 
SPRI NT 

ACC 

R *- 8 

<2200 

D ER I Y 
FUNCT 

AMO 

r* e 

920 3 

Of/ ASS 

F M 0 

R*8 

9203 

CETD 

TRAJ 

POV/R 

R*6 

9232 

TAP 
LOAD 
DER I V 
FUNCT 
SOLA? 
S VST D 
CDERI Y 
OPRI NT 
P E T I N J 
SPRI NT 
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COMMON 

PE ALP 

(CONT IMIJED) 


VAR I AOL E 

type 

A DDR 

SO BP OUT I N 

D P 0 1! ft 

fj * g 

9 24 0 

TA? 
D ER I V 
FUNCT 
SOLAR 
CD BP I V 
89 R I NT 

T A UP GW 

R a 6 

9 £4 8 

T A 3 
D EP I V 
F UN c r 
6 OL A ft 
C D ER. I </ 

dma:-,s 

R * 8 

9256 

0 Eft 1 V 

0 M OV C 

p / 0 

9 264 

GET 3 
0 MA S S 

AM 

St'K 

9 27 2 

GET 3 
RETINJ 

EX 

R <- 3 

9 28 9 

TP A J 
PE 7 I NJ 

FP 

r a a 

9288 

TR AJ 
RE 7 1 NJ 

GCST 

r c e 

9 296 

TRA J 

P.ETI NJ 

VH 

p* e 

9 3 C 4 

GET 0 
TRA J 
F-ETINJ 

VHS 


9 31 2 

TR A J 
PETI NJ 

V INF 

ft* e 

9 520 

GE TO 
TP A J 

petinj 

V J P F T 

R* 8 

9 326 

TRAJ 

RETINJ 

VS 

R+ 8 

9 336 

TRAJ 

RETINJ 

XU.NG 

R ~ u 

9344 

03 TART 

Y L ON G 

R * 8 

934 4 

SPPI NT 

AV J 

i<* : 8 

9352 

T A 9 

TRAJ 
DERI V 
F UN Z T 
SV.STD 
TRA JI 
SPRI NT 

FTV J 

R *• 3 

9360 

TRAJ 
TRA JI 

FTCVJ 

R-*6 

9363 

DSP I V 
TRA JI 

r-MS I 

R *b 

5376 

GET 3 
TRAJ 
G M ASS 

SKOUNT 

r* e 

9384 

LOAD 


Ill 



COMMON 


REALS (CONTINUED) 


VAR I ABLE 

TYPE 

ADDS 

SUBROUT I M z 
C HEC < 
INTERS “ 

DM AX 

R * ft 

9 39 2 

SOLAR 
SPRI NT 

AK 

r* e 

9 4CG 

GET 3 
TP. A J . 

D 

i« 6 

940 6 

CD ER I V 

DENS I T 

r* e 

940 e 

T A 3 


LOAD 
D£R I v/ 
F UN Z r 
SOLAR 
S '//ST 3 
SPRI NT 



PMO 

R 6 

9A1 6 

GET 3 





TP A J 

A* 

P M 9 D 

ft *8 

942.4 

GETO 
TRA J 
TRAJI 

TEST 

R>f 8 

9432 

TRAJ 

TRAJI 

GSUOX 

ft* 8 

944 0 

GETG 

T EMP 8 

P. * H 

94 4 8 

TRAJ 

GETO 

TRAJ 




D PC ‘ADD 

R * 6 

9 4 56 

T A 3 





EUNC r 
SOLAR 





CDcR I Y 


TEMP« 

R* a 

94 64 

GET3 

TRAJ 


A JT 

R v o 

94 72 

TRAJ 


A JPP 

R * 8 

9480 

GET 3 
TRAJ 


FETA 

R *8 

9488 

GET 3 
TRAJ 


SCALE 

R* 6 

9496 

GET 3 
TRAJ 

- 

TP OWE R 

R* 6 

9504 

BEGIN 

i np j r 

SOLAR 





QPRI NT 


DEGRAD 

R* 3 

95X2 

FUNC T 
SOLAR 





SPRINT 


P M A X C 

R» 3 

952 0 

S OLA R 


RPNAXO 

R* 3 

952 8 

SOLAR 

OSTART 


RPC a C 

Rt6 

9 53 6 

SOLAR 





QSTART 


DP OWD 

R *8 

9544 

SOLAR 
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COMMON |<EAL-, (CONTINUED) 
VAR 1 ABI £ TYPE ADD« 

V; p : ? I M R * B 955 2: 

SX SR*P. 960 0 


x k * s ? c c c :• 


^ V/'Oi > 

... sq 


a, 

- ■ %* . 


XD P-=3 iC4CC 


S'JB 9 OUT 1 Nl 
CD CP I M 
SP PINT 
GS T 3 
TRAJ 
TA 3 
LOAD 
G UE S S 
PRIOR 
AN ST E 3 
i m pp n r 
I MPU_ S 
INTER 3 
OS TAP r 
RK STEP 
T A 3 
GET 3 
LOAD 
MOP S 
TKIA J 
C HEC < 
DtR I V 
F UNO T 
G PT R V 
G UE S S 
PRIOR 
PUNCH 
RADAR 
SOLAR 
SWING 
SWOT 3 
7 PA J I 
ANSTE 3 
CDER I V 
CON VRT 
I MPRNT 
I M PULS 
I N T C 3 3 
OPPI NT 
05 TART 
PK.STE? 
SPPI NT 
SWTRA J 
TARGET 
VP PI NT 
TRAJ 
DERIV 
FUNCT 
PRIOR 
ANGIE 3 
CDER r V 
OSTART 
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COMMON 


REAL 3 (CONTINUED) 


VARIAOLE TYPE A DDR SUER CUT IN 





F- K S” 3 3 
SP PINT 

S BET A 

R~-3 

I C8CC 

PRIOR 
AN5TEP 
I NTERP 
RK5T= 3 

BET A 

Rv S 

1 0 30 6 

CHEC.< 
PRIOR 
TPA J I 

AN STEP 
COcR I Y 
RKSTE=> 
TAPS 3 T 

r T A u 

R " £> 

1 C 6 1 6 

PRIOR 
PK STEP 

TAU 

R '~ 8 

1002 4 

GET 3 
PRINT 
PRIOR 
P UN C -i 
TRAJI 
PR INT R 
OPRI NT 
RK ST E 3 
SPRINT 

X I NT 

R v S 

1 J COO 

GET3 
TRA J 
UP PI NT 

XD INT 

R* 8 

1 3000 

TR A J 

X T INT 

R * e 

15 0CC 

GETS 
TR A J 
QPRI NT 

X T D I N T 

R * C 

1 5240 

GET 3 
TRA J 
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COMMON 


variable cross reference T A 3 l 


COMMON 

F XTPEM 

L E NS 

T H 217 5 

VAR I A C: L 

TYPE 

AO DR 

£ OOP CUT 1 N 

r ~ « 

R: F d 

o 

B EG I N 

TIME 

O is o 

r 

E XT A 3 
STORE 

T 04 V 

P* 3 

30 0 

E XT A 3 
S TORE 

D I 3 T 

; P - 0 

1 6 CO 

E XT A 3 
STORE 

ON OFF 

P - 0 

2 .2 0 C 

c XT A 3 
P R 1 N r 
S TORE 

01 S COM 

E* 3 

4 80 0 

E XT A 3 
STORE 

AN G C 0 M 

R * e 

64 OS 

E XT A 3 
S TO R E 

AN GPS I 

R *■ 8 

S C 0 o 

E XT A 3 
STORE 

AN CT HE 

IT •< 8 

0 60 0 

E XT A 3 
STORE 

AN GO HI 

R A d 

1120 0 

E XT A 3 
STORE 

P OWE X 

R*8 

1 2 80 C 

E XT A 3 
STORE 

CM I X 

R ~ c 

1 4 4 C C 

C HEC< 
E XT A 3 
STORE 

AK O JN'T 

r -- e 

16 0 0 0 

E XT A3 
S TORE 

COPS 

P. £j 

.1 6800 

LOAD 
C. HE C < 
STORE 

R 

R * 8 

2 3. 2 60 

T A 3 
LOAD 
C HE C < 
C D EP I V 
I N TER 3 
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COMMON VARIABLE CROSS RE^E ?ENC£ TA3LZ 


COMMON GUNCuiV LENGTH 6 4 


VAR I A 01. F 

TVPE 

A DDR 

SU8.R OUT I NZ 

VO 

k*-8 

C 

0MAS5 

GUNT.hr 

V 0 D I V 

R? S 

6 

DMA 3 5 
GUN'THR 

voo 

R=ce 

16 

GET I 
GRASS 
GUNTUR 

ICC 

Rtf 6 


GET I 
OMASS 
GUNTHR 

ov 

R ~ 8 

22 

OMASS 

GUNTHR 

0 V V C G 

ft? 0 

4 0 

GRASS 
GU NT HR 

0 V I c 0 

R'-'b 

48 

GET! 
GE T D 
0 R A S 5 
GUNTH? 

GL I 

L A 

56 

GET 3 

L D V V A 0 

L* 4 

56 

GET I 
OMASS 
GUNTH < 

L D V I C 0 

L* 4 

6C 

GET I 
GRASS 
G'J NT H R 
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COMMON VARIABLE CPOf.S REFERENCE TA3L 


COM MON 

I NT GR4 

l_ 

EN3TH 4 f- 

V AS-; I Af.d_ £ 

1 YP E 

ADDP 

SUBRCSUTI 

I f 'L 

1 * 4 


rJtGIM 
I N : P J 7 
QSTAy T 

I CAVE 

I >' 4 

r 

5U MMR Y 

MODE 

1 >: 4 

4 

TAP 

LOAD 

0 EG I N 
CHECK 
E XT A 3 

1 t iOjT 
P UN C H 
E CL A R 

GO TART 

IPK 

I -- 4 

0 

3 EG I N 
I NPUT 
PUNC 3 

1TF 

I * 4 

1 Z 

b EG I N 

I HP J T 
T I KT 0 < 

L I N E 

I * 4 

16 

PRINT 
E. CL A R 
MI NM X 2 
0 S TART 

JP A 1 NT 

i - ’ 4 

SC- 

TP A J 
I NP J r 

M BOOST 

14-4 

24 

I NPUT 
PUNCH 
OP PI NT 
GC TAPT 

MOP 7 

I»4 

23 

WA I vl 
I NPUT 
PUNCH 
CORNER 

NT APE 

1*4 

3 ? 

3 EG I 'J 
I NPUT 
P UN C H 
QS1 AR r 

M DAY 

I * 4 

36 

B EG I N 
I NPUT 
P UNC H 
OSTART 

MONTH 

1*4 

4 C 

B EG I N 
I NP JT 
PUNC H 
OSTART 

M YEAR 

I a 4 

44 

BEGIN 
I NPUT 
PUNCH 
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COMMON INTGR4 ( C ON T I NUEC ) 


V AP I Ai.jL E 

TYPE 

ADDR 

SUBROUT I i 
0 S T A k T 

M R r- A D 

I * 4 

4 8 

MA I N 
I NP-JT 
CORNER 

MUPl.AT 

I * 4 

5? 

B EG I N 
I NPUT 

corner 

FI N I S 4 

f«PP 1 NT 

I * 4 

56 

1 NP JT 
FI NISH 
OS TA R r 

NSET 

I v- 4 

60 

MA 1 N 
BEGIN 
I NP J T 
PRINT 
P R I N T R 
OP R I NT 
OS TART 

i nor 

I * 4 

80 

I NPUT 
PUNCH 
CORNER 
QSTA R T 

jpp 

I * 4 

84 

TRAJ 
I NPUT 
PUNCH 
OP PI NT 

JT 

1*4 

88 

TRAJ 
I NPUT 
PUNCH 
OP PINT 

MOPT2 

I » 4 

9 ? 

TRAJ 
I NPUT 
PUNCH 
TRAJI 
EF MPRT 
QPRI NT 
OSTART 
SPRI NT 

M OPT 3 

I * 4 

96 

MORE 
TRAJ 
8 EG I N 
I NP JT 
PUNCH 
TRAJI 
QPRI NT 
Q START 
SPRI NT 
TAPSET 

MOPT 4 

1*4 

ICO 

MORE 
I NPUT 

NTAi-iG 

I * 4 

140 

SWI N 2 
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I igT GP.'t (CONTINUED) 


COMMON 


VAR i A“!L r 

TYPE 

AC- DR 

SUER CUT I 
5 * TR A J 

MT M ASS 

I 4 

t 4 A 

I NPJT 
P UN C 3 
C P K l N r 
Q START 

NSPFX 

I « 4 

1.48 

CHEC < 
H XT A 3 
PRINT 
.3 TG R z 
TRA JI 
TAPS FT 

MPUMCH 

1*4 

. 1 5o 

I NPJT 


P UNC H 
F 1 N I S3 
I NTF.P=> 
0 .START 
SPRINT 


L COUNT 

I * 4 

7.6C 

ostart 
SPRI NT 

NR R I NT 

I * 4 

1 64 

T A 3 


TP A J 
5 LG I M 
I NPUT 
PRINTS 

CL TART 


NiPS 

I ft 4 

3 e. b 

HI Nt SH 
MI NMX 3 
OPR] NT 
QS7ART 

NS VI PAR 

l >;■ 4 

1 8 4 

3 EG I N' 
I NPJT 
MI NMX 3 

NORMAL 

I • 4 

168 

I NPJT 
OSTAk t 

LI MPH I 

I x-4 

19 ? 

GET a 
SETUP 
TRAJl 

KP ART 

I v A 

1 96 

MA I N 
I NPUT 
CORNER 

I H'J MG 

I ft 4 

POO 

TRA J 
C HEC < 
CORNER 

NHUNG 

I V 4 

204 

3 EGI N 
I NPJT 
CORNER 

K DUiMT 

I ft 4 

21 ?: 

MA I N 
E XTA3 




PRINT 
P UNO 3 
CORNER 
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COMMON I NT GW 4 ( COM! I NUel) ) 


VARI A3L. r. 

7 yp a 

A DDR 

SUDR OUT I 
EF MPRT 
FI NISH 
OP RI NT 
GSTAST 
SPRINT 
T I KTO< 

1 OAL 

14-4 

21 6 

I NPJT 
PUNCH 
OS TAR T 

MV C ARD 

I * 4 

220 

OPR I NT 

.MYEARA 

I * 4 

224 

QP PINT 

MONO 

1*4 

22 8 

OPR I NT 

MON A 

I* 4 

222 

QP R 1 NT 

Ml A VO 

I * 4 

23 6 

OPR I NT 

M DAY A 

1*4 

24 0 

OP PI NT 

L X X 

I v: 4 

244 

TPA J I 
C P R I N T 
OSTAP r 
T I KTO< 

MX X 

1*4 

24 8 

GE T 0 
OPRINT 

ostart 

J J 

1*4 

252 

TAR 
C HE C < 

KF 

1*4 

25.4 

T A 3 
CHECK 

J H'JNG 

1*4 

2C0 

MA I N 
CORNER 
FINISH 
71 K TOC 

MAJOR 

1*4 

264 

TRAJ 
P PINT 
MI NMX 2 
PR INI R 
T I KTOC 

MAJORS 

I * 4 

2 68 

MI NMX 3 
PP INTR 
71 KTOC 

M I NOR 

I X- 4 

272 

MI NMX 2 
PR I NT R 

LAUNCH 

I * 4- 

276 

I NPJT 
PUNCH 
S ET IP 
OS TART 

. I STAR 

1*4 

28 0 

PUNCH 

I END 

I *4 

284 

P UN Z H 
SPRINT 

1 SP I N 

1*4 

286 

I NPJT 
S CL A R 

I OUT 

1*4 

292 

GE T 2 
I NPUT 
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( C Ci NT I NUF. D ) 


COMMON 1 NT Ci<* 


V AR 1 AGL 

T YP E 

4U0!-' 

5 GBP OUT I 
PR I N T 6 
OS TAP T 

N Vj'.V 

I * 4 

296 

P P 1 N T 
S WS T J 
T PA JI 
M I KM X 3 
PR I N T 2 
T A PS E r 

I S a 

I -• 4 

A 0 c 

S IV 5 T j 
P P INI 2 

1 nTPR 

I - 4 

550 0 

I NPIJT 
I NTS 2 = 

LEG 

I v 4 

55 * 4 

; 2 A J 

S P K I NT 
TA PSET 

l.L. GMA X 

I 4 4 

5 OS 

GET 2 
TP A J 
SE.TJP 
E F MPR r 
Pk I N T 2 
OPPI MT 
0 START 
SF PINT 
T A PS E r 

MOOT X 

1 4 

51 2 

TP A J 
I NP J T 
OP RI NT 
GS TART 

sprint 

7 A PS £ r 

I N7 E r< 

1*4 

52.2 

TP A J 
GSTA2T 
SPPI NT 

LOAOX 

1*4 

55? 

I NP 0 T 
0STA2T 

MPiJvV 

1*4 

55 6 

E XT A 3 
I NPUT 
GPRI NT 
0 START 

N S V, I N G 

I * 4 

560 

MOP z 
I NPJT 

NDIST 

1*4 

564 

3 cG I vi 
I NPUT 
R AO A 2 

N P ii f< F 

1*4 

56 3 

I NPUT 
OS 7A2T 

M P n ° F 

i -4 

572 

MA 1 7 
OSTART 

f 7 PR' 17 

1*4 

57 6 

I NP U i" 
MI NMX: 

INTERX 

1*4 

560 

MOP E 
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COMMON -1N7GP:4 (CONTINUED) 


VAkIAl-LS: TYPE ADDR SUBROUTINE 





swing 

SPRINT 
S W TR A J 


iM A XH AM 

17 4 

564 

0 EG I N 

1 NP J r 
SwSTO 


NLEAVc 

r~4 

560 

£ V. I N 3 
SWTRA J 


MS '/< IN6 

I x A 

59? 

MORE 
I NPUT 


LL 

I *' A 

74 0 

TRAJI 
PARI NC 
G START 


MM 

I + A 

I 02 C 

GET0 
M I NM X 3 
OS TAP r 


MPH I 

I *4 

1 50C 

TA3 


NP H 12 C 

l *• 4 

1 30 4 

T A 3 
DtP I Y 
PRIOR 
RKSTc 3 


NCHK 

I a 4 

i 20 6 

C HE C < 


NCOP 

I * 4 

1 31 2 

LOAD 
C HEC < 
STORE 


NSTtP 1 

I ~ 4 

1 316 

TA 3 
P R I N r 

TRAJI 


NST&P2 

I f 4 

1 320 

T A 3 
PRINT 
T RA J I 


JC 

1 + 4 

1 324 

TA=> 
LOAD 
TRA J 
CHECK 
CDEPI Y 
TAPSE7 


JCMAX 

I 'r 4 

1 32 0 

TA 3 
LOAD 
TRA J 
CHECK 
TAP3ET 

- 

NC.'/X 

1 + 4 

1 400 

PRI NT 
M 1 NM X 3 

4 
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COMMON VAR I AGLn CK CSS- 


REF Lf< ENCE 


A 5 L 


COMMON 
V AR 1 A ML L 

rj K 


OY 


CL NX 


CONY 

O 


TEFAT 


L E N 3 T H 7 A A 


TYPE ADDR 

r * n c 


PC 


>5 


peoo 


4480 


Pvfe 5C4C 


« 


560 0 


K* 6 6 J 60 


SUbP OUT I N 
TRA J 

3 E G I 
G UE 5 0 
I NPUT 

p (> r n r 

P UN C H 
SET J 3 
CORNER 
FI NI S H 
PA PI NS 
GPP I NT 
G START 
5 EG l N 
I NP JT 
5ETU 3 
CORNER 
PR INTR 
QPRI NT 
05 TART 
GETS 
TP A J 
B EG 1 N 
G LIE ;> i 
S CL A r 
FINIS i 
PAP I NC 
PR IN TP 

os tar r 
S EG I N 

p « i n r 

PP I NTS 
05, TART 
T A 3 
GF7T0 
TP A J 
GCTRV 
O MA S S 
F: ADAR 
SOLAR 
T RA J I 
CDERI V 
FINISH 
PRINTS 
QPPI NT 
OSTART 
SP PINT 
VP PI NT 
TP A J 
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COM N’C'i'l I TE PAT ( C.ONT 1 NUF.D ) 


V,U'I AflLE TYPE 


FXL R-'-8 


F XL 1 R*8 


A DDR SUBROUTIN 
PRINT 
PUNCH 
5 OL A ^ 

S W I Ni j 
F I NISH 
OP PI NT 
6 7?0 OF TO 

PR! NT 
PR I NTT 
OSTAP T 
728“ MINMX3 

PP INTP 


124 



COMMON 1 

C OM i-‘C : J 

VAR I ABL 
E RROH 


l ft n or x 

CO NVRG 


f- ALEE 

F I X P Li .V 
iJ' J R 

KONOtR 

HUNG 

PARNS W 
QOECL 

TUDFLG 


VARIABLE CROSS REFERENCE TABLE 


_ GG 1 C4 


LENGTH SOOT 


TYPE 
L *■ A 


L *4 


L'-A 


L 4 
L >■ A 


LA A 
LA A 


L w A 


L * A 
LV A 


L* A 


ADDA SUBROUTINE 
G T A 3 

GET 1 
MA I N 

a eg i m 

C HE C < 
PUNCH 
ft OL A R 
5 tiS T D 
FI N I SH 
OP PI NT 
ft.U MM R Y 
T I KTCX 
C TK A J 

7 RAJ I 
A MA I M 

3 EG I N 
PUNCH 
FI NISH 
C)P PI NT 
O START 
SU MM A Y 
T I KTCX 

ft begin 

ft GE T 0 

TP A J 
PRI NT • 



OP RI NT 
OS TART 

1 2 

MA I N 
OSTAPT 

16 

TR A J 
CHECK 
MI NMX 3 
PA RI NC 
OSTAP r 

£ C 

CORNER 
FI NI SH 
Ml NMX 3 

£4 

M I NM X 3 

2 G 

O MA S S 
PRI NT 
TRAJI 
OP RI NT 
OS TAR r 

” O 
*-> C. 

T A 3 


CHECK 
0 EF? I \> 
PRI 3P 


125 



COMMON LOGICS (CONTINUED) 


V AG I AbL H TYPE 


CUTE CL L ■' 4 


FI XT HR L v 4 


A DDR SUBROUTINE 

R AO A R 
UE CL IN 
G S T A R T 
thangd 
36 GET 2 

TR A J 
R AO A R 
TPA j I 
QPRI NT 
OS TAR T 
SPfil NT 
AC T A 3 

GE T D 
CHE C< 
DER W 
r UN G T 
RRI OR 
TRAJI 
CDtRIV 
QPRI NIT 
OSTART 
R K ST t£ 3 
SPRI NT 


QVLOSS 

L : - A 

A A 

TRAJ 

OSTART 

FL YBY 

L * 4 

A 3 

GETD 
TRAJ 
GPRI NT 
OSTART 
RETINJ 

BP. AK E 

L< A 

52 

QPRI NT 
OSTART 

RENCEZ 

LA A 

56 

OSTART 
RET I NJ 

POSVEL 

L- r A 

60 

OSTART 
RF. TI NJ 

SP IR AL 

L* A 

64 

TRAJ 
QPRI NT 
0 S TA R T 
RETINJ 

VELOSS 

L* A 

6 6 

TRAJ 
OSTART 
RETI NJ 

REG ION 

LA A 

72 

TA 3 
C HE C < 




DER I Y 
FUNCT 
SOLAR 
OSTART 

ERODE 

LA A 

76 

TA 3 


LOAD 
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CUM I*: )M 
VA f- I ACL ti 


F L 4 ° L ! 4 

i-P'IN L*4 

TILT L * 4 

HCUSE L ■' 4 

TRACK L 4 


PL AN FT L - 4 


0 J E X L > v 4 


FT XT A U L T 4 


SUSP OJT I N 
CHECK 
D ER I </ 

[-: xr a a 

F UN I I 
P PI OR 
S OL A R 
tpaji 
COER IV 
O 3 r<I NT 
F K ST E 3 
SP PINT 
E T 
C< 

otp r v 

S CL A R 
S OL A R 
C HE C .< 

S OL A R 
SCLAR 
T A 3 
S V INS 
FIM5H 
OS TART 
SP PINT 

ge r y 

7R A J 
RADAR 
T PA J 1 
FI Im I S -I 
OPR I NT 
US TART 
SP F? I NT 
TAP 
TP AJ 
C HE C < 
PRINT 
SOLAR 
STORE 
S W I NS 
TPAJI 
COOP I V 
CORNER 
FI N1 SH 
INTERS 
PR IN TR 
GSTARI 
10 8 TA 3 

CHECK 
F UN S T 
T IRA J I 
TAPSET 


L 3 G I C 4 ( C O N T I N U E 0 ) 

TYPE A DOR 


00 


TA PS 
CHS 


84 

88 

9 2 
9 6 


1 0 r 


1 C4 
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COMMON LOGIC* (CONTINUED) 


VAC I A3L K 

type 

A DDR 

SUGPCUTI 

He AT 

L * A 

1 12 

TAP 
LOAD 
C HE C < 
F UN C I 
SOLAR 
T RAJ I 

L OG5 E 

L i( t 

.1 J. 6 

TP A J 
T RAJ I 
OP PI \'T 

F I P 5 T 

L*4 

: 2 •: 

TAP 
LO AO 
COER 1 

CO ATT 

L»4 

*24 

TAP 
Lu A 0 




ST E 3 
CHECK 
r UN c r 
S *STO 
COER I v/ 
INTER 3 
SPPI NT 

OF ROD IE 

L v 4 

128 

TAP 
DER 1 Y 
PRIOR 
RKSTE? 

RR Z C R 0 

Lv4 

13? 

GET 0 
TRA J 
T RA J I 

X L 0 A D 

L* 4 

1 3£ 

TRA J 
05 TART 

'A I PL 

L * 4 

140 

TAP 
OS TAR r 

P C'JRV 

LV.-4 

144 

TAP 
LO A 0 
C HE 2 < 
DEP I V 
fun:i 

SOLA R 
TRA JI 




COER l v' 
05* TART 
SPRINT 

MA XPOW 

L * 4 

146 

LO A 0 




CHECK 
S OLAR 
TRA J I 
QS TAR T 

EDGE 

L v 4 

152 

TAP 
CHECK 
S OL A R 
£PRI NT 
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COM FLU LOGICS (CONTINUED) 


AN I AL L H 

t Y ?" 

ADD;' 

su hb out i 

PLUS 

L * A 

1 S C 

T A ? 
CH8C< 
0 66 I V 
TH AN GO 

JUMPED 

L. ••• A 

: 6 c 

T A 3 

5 V;S T 3 

Al. V.-AYS 

L - A 

1 6 4 

GC 7 3 
TP A J 

L DC CL 

L r ft 

1 f: 3 

0 K 1 F? 7 

1 k i n r 

7 PA J I 

BALL I G 

4 

172 

T A 3 
C Ha C < 

GUI f 

L * 4 

1 7 1 

S V; I N O 
T I K 7 0 < 

F XTO A 

L. *• 4 

\ 8 0 

MOPE 
FI NI SH 
GPP I NT 

0 f/i Qi? F 

L >■ 4 

3. 3 4 

MO P E 
GPP I NT 

f A N Dr M 

L 4 

i b a 

MOP E 
5 >« I N3 

A * A 

L * 4 

6C C 

OS T 3 
SETUP 

A ‘i 0 

L 4 

60 4 

GET 3 
TRA J 
.6 F T U 3 

Al C 

L* 4 

60 8 

GET 3 
3ETU 3 

a :> a 

L v 4 

61 2 

GF.T 3 
S ET J 3 

A > B 

1.4 4 

6 5 6 

Gt: T 3 
TP A J 
SET J 3 

ABC 

L'- 4 

62 C 

GE T 3 
S ET U 3 

A A 

L* 4 

624 

OF T 3 
SE T J 3 

A 3 C 

L'.4 

LES- 

GET 3 
SET J 3 

A4 A 

L. - 4 

OVS 

GET 3 
S ET J 3 

A A 3 

L * 4 

636 

GET 3 




SETU 3 

A AC 

L---4 

64 C 

GET 3 
S ET J 3 

A 5 A 

L=> 4 

64 4 

GET 3 
SETU 3 

A. SB 

L * 4 

64 ft 

GET 3 
SETS 3 

ASC 

L ^ 

652 

GET 3 
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COMMON LOGICS { C ONT I RUE D ) 


VAR T A'iL 5 

T Y P .3 

AD Ok 

LUBROUr I N 
3 E T U 3 

A 6 A 

L 4 

656 

Gt r T 3 
5ETJ3 

A 6 -3 

L •- 4 

66 C 

GET 3 
SET U 3 

ACC 

L -■ a 

664 

GET 3 
SET J 3 

A'/A 

L ~ 4 

665 

GET 3 
S ET U 3 

a y i-j 

L* : 4 

67? 

GET 3 
S E T U 3 

A 7 C 

L * 4 

67 6 

GET 3 
GET J 3 

- A <3 A 

L x- 4 

68 C 

GET 3 
SET J 3 

A B 3 

L * 4 

684 

GET 3 
SET J 3 

A • 1 A 

L v 4 

68 B 

GET 3 
S E T U 3 

A * ", B 

L - 4 

69? 

get 3 

SET J 3 

All C 

LA 4 

696 

GET 3 
TRA J 
S ETU D 

At C A 

L ' A 4 

7CC 

GET 3 
SET 3 = 

AT 5 3 

L • - 4 

704 

GET3 
S E.T J 3 

A i 3 A 

L*A 

708 

GET 3 
TRA J 
S E T U 3 

A l 3 B 

Ls 4 

71? 

GET 3 
SET J 3 

A! A A 

L * 4 

716 

GE T 3 
5ETJ 3 

At AB 

L*4 

72 0 

GET 3 
3 ET J 3 

A 1 4 C ’ 

L v 4 

7? 4 

GET 3 
SET J 3 

A 1 5 A 

L*4 

72 . 6 

GET3 
S ET U 3 

A*. :5B 

L-T4 

73? 

GET 3 
5 ET U 3 

A 1 6 A 

L4 4 

726 

GE T3 
S ET U 3 

A “ 6 3 

L Y 4 

74 0 

GET 3 
SET J 3 

A 1 6 C 

L ~ 4 

744 

GET 3 
S E T J 3 

APMI 

L * 4 

130 

74 B 

GET 3 
SET J 3 
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CG> 
\J A 1 


■i 0 U LOG JC4 ( C O N T I NUB 0 ) 


i a,-’l f: 

TVPt 

AODK 

SUB 9 CUT I N 

A i 7 A 

L - A 

A? S 

BE T 0 
SETiJ=> 

A ;■ 7l.i 

L “ 4 

3 ? 

GE T o 
SET U 3 

A * • • A 

L *" 4 

636 

GETO 
TEA J 
SET 3 

A 1 •'} 0 

L 4 

34 0 

Gt.- TO 
E f-Tija 

A ; 9 A 

L “ 4 

34 4 

GE T 0 
TP A j 
SETa-"* 

A .i U 3 

L* 4 

«4 6 

GET 3 
SET J 3 

A £ A 

l. v. 4 

85 9 

GETS 
TP A J 

SET 

A '0 6 

L v- 4 

656 

GE T 0 
SET J 3 

A 4 i A 

L 4 

ft CO 

GE T 0 
SET J 3 

A A 2 A 

L ’* 4 

3 6 4 

GET 0 
G E T U D 

A 4 3 A 

L-; 4 

86.6 

GETO 
S ET J :> 

A A 4 A 

L * 4 

6 7 3 

GETO 
SET J 3 

A A A B 

L * 4 

376 

GETO 
SET JA 

A A 5A 

L 4 

38 0 

GET 0 
SET J ^ 

A A 58 

L ■< 4 

8 8 4 

GETO 

SET'J=> 

A 4 6 A 

L 4 

3B3 

GETO 
SE T J 3 

A A 6 B 

L >■ 4 

696 

GE T 0 
S ET !P 

A 4 A 

L * 4 

3 96 

GETO 
S ET J ^ 

A 0 X 

L -• 4 

9G A 

TP A J 
0 MA S S 
P R I N r 
TPAJI 
PR I NT P 
GSTART 

A. BY 

L - 4 

>180 

GETO 
0 “A S S 
P PINT 
SETJ3 
PP I NTT 
QSTAST 
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COM N'CN LOG 1 C4 ( c ONI I NUE'D) 


V A? 


I ABL E-: 

TYPE 

A OOP. 

SUBROUT I 

A 4 7 8 

L+ 4 

1 60 0 

GET 3 
SET J 3 

A4 8A 

L * 4 

1 6C4 

GET 3 
SETJ 3 

A 4 8 B 

Lv& 

3 6C8 

GST 3 
SETU 3 

A 5 1 A 

L * 4 

3 63 2 

GET 3 
SETJ 3 

A 5 ?. A 

L^4 

1 6 3 6 

GET 3 
SETUP 

A 5 3 A 

L 4 

1 620 

GET 3 
SET J 3 

ALA A 

L •" 4 

1 624 

GET 3 
SET J 3 

ALA B 

L* 4 

3 6 28 

GET 3 
S E T U 3 

A55A 

L* 4 

3 632 

ge: t 3 

SET J 3 

A 5 5 3 

L * 4 

3 636 

GET 3 
SETU 3 

A 5c A 

L* 4 

3 64 C 

GET 3 
SET J 3 

A56B 

L * 4 

1 64 4 

GE T 3 
£ ETU 3 

A57A 

L* 4 

1 64 5 

GET 3 
SETU 3 

A 6 73 

L 4 

1 6 5? 

GE T 3 
SETU 3 

A 58 A 

L* 4 

1 656 

GET 3 
SETU 3 

A 5 88 

L* 4 

3. 66 C 

GET 3 
SETU 3 

A t. 1 A 

L 4 

3 664 

GE T 3 
SETU 3 

A 6 3 A 

L 4 

666 

GET 3 
SETU 3 

A 6 3 A 

L v 4 

1 67?. 

GET 3 
S ET J 3 

A 64 A 

L * 4 

1 676 

GET 3 
SETUP 

A 6 4 B 

L ^ 4 

3 680 

GET 3 
SETU 3 

A 6 5A 

L 4 4 

1 684 

GET 3 
SETJ 3 

A 6 5 B 

LA 4 

3 668 

GETS 

SETU 3 

A66 A 

L »:• 4 

1602 

GET 3 
3 ETU 3 

A 666 

L “ 4 

1 6 96 

GET 3 
SETU 3 

A 6 7 A 

L * 4 

170 0 

GE T 3 
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COMMON 


L nr, I C 4 ( CON 7 I NUO D ) 


am i apl c 

T Y.-E 

ADD-: 

so or-: c:ut l 
SOT J 3 

A 7 B 

Li A 

1 70 4 

GOTO 
SOT J 3 

A A 

L * 

2 7 C A 

GO T 0 
SOT U 3 

Ac U8 

L - r 4 

1712 

GO. T 0 
SETUP 

A 4 -J A 

L A 

1716 

GOTO 
S ET U 3 

A f: A 

L~4 

1 7 2C 

GETO 

SETUP 

A i 7 A 

L" A 

2 724 

GOTO 
SET J 3 

AcO A 

L. * A 

1726 

GOTO 
S E T U 3 

< 

iTv 

< 

L.--< A 

1 7 3 2 

GO T 0 
SETUP 

A 7 0 A 

L*A 

1 736 

GOTO 
SET J 3 

A 3 0 A 

L -v 4 

1 74 C 

GETO 

SETUP 

A 7 0 U 

L-4 

1 74 6 

GETO 

SETUP 

A • 7 A 

L 4 

i 74 e 

GOTO 

SETUP 

A ! ■' *5 

L* a 

17 5? 

GOTO 

SETUP 
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COMMON VARIABLE CHOSE REFERENCE TABL£ 


COM M U * -1 

i> CiL G Y 3 

L EN3 

T H E 2 4 0 

V Af< I A B L. E 

TYPE 

AuDP 

E UBP OUT I ME 

GM 

H ^ c 

c 

FF vi 
TP A J 
BEG I N 
5 V. 1 flj 
QSTAPT 

P A l. I J S 

R* 8 

563 

TRA J 
BEGIN 
S V» 1 N 3 

APL 

f< * 8 

1 120 

BO P E 
a EG I VJ 
S W I M3 
f:f mppt 

GPRI NT 
bp p. i m r 


B. HILTOP INPUT 

1. NAMELIST . Inputs to HILTOP are given through the NAME- 
LIST feature of the IBM Fortran IV programming language. The input NAMELIST 
is named MINPUT, and every input required or used in the program is declared 
by name in the list. The general form for assigning an input value to a quantity is, 
simply, 

NAME = VALUE 

where NAME is the name assigned to the variable and is included in the NAME- 
LIST, and VALUE is a numerical or logical quantity consistent in form (i.e. , 
logical, integer, or real) with NAME. Unless otherwise specified, all MINPUT 
names commencing with one of the letters I through N represent integers, 
whereas all names commencing with one of the letters A through H or O 
through Z are double precision floating point numbers. Each NAMELIST case 
must begin with the characters 

& MINPUT 

commencing in card colume 2 and followed by at least one blank, and end with 
the characters 

&END 

preceded by at least one blank. Card column 1 is ignored on all NAMELIST 
input cards. Multiple data assignments on a single card are permissible if 
separated by commas. Blanks in the variable field, VALUE, are taken as 
zeroes. A comma following the last VALUE on a card is optional on the IBM 
system. The order of the input data assignments is arbitrary; i. e. , they need 
not be in the same order as listed in the NAMELIST. In fact, there is no re- 
quirement that any specific input parameter be represented in the input data 
set. If no value is included in the inputs for a particular parameter, the de- 
fault value is used (see Default Values). For other details regarding the 


135 


NAMELIST feature, the reader is referred to the IBM System 360/Fortran IV 
Language manual, NAMELIST cases may be stacked back-to-back indefinitely. 

2. Definitions of Input Parameters . Specific examples of the 
program inputs are given in the Sample Problems and Results section. Default- 
values of inputs are given in the next section. 

The program inputs, in alphabetical order, are: 


AAI 


AE 

ALPHAA 

ALPHAT 

ALTITU 


AN 


Desired final extra-ecliptic inclination, i. Related to AE, AR, 
and IOUT. [deg] 

Desired final extra-ecliptic eccentricity, e. Related to AAI, AR, 
and IOUT. 

Specific mass of solar arrays, a . [kg/kw] 

3 

Specific mass of power conditioning and thruster subsystem, a 
[kg/kw] 

This input variable is associated with program logic which has 
not been kept up-to-date, specifically, logic pertaining to opti- 
mum departure of a NERVA-type rocket from Earth orbit. This 
variable should be ignored. 

Trajectory-integration exponent n in expression (39). 


AR Desired final extra-ecliptic perihelion distance, r f . Related to 

AAI, AE, and IOUT. [AU] 

ASOL Array of five elements consisting of the solar power law co- 

efficients a^ in expression (18). ASOL(l) > 0 tells the program 
to use the input coefficients rather than the internal coefficients. 
The coefficients are normalized internally, and the program exe- 
cutes the iterations to produce the required remarkable points of 
the power curve (which are printed), 

BI Efficiency coefficient b in expression (16). Related to DI and El. 
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B1 

B2 

B3 

CNI 


CNIX 


CSTR 

CTANK 

CTRET 

DI 

DMRETR 

DPOW 


ECI 


ECIX 


El 


Launch vehicle coefficients b , b , and b in expression (2). 

X £ O 

Used only if MBOOST is negative, 

[kg, m/sec, kg] 

Inclination to ecliptic of primary-target orbit. Input only when 
MOPT3 = 11. Related to ECI, OMI, SAI, SOI, TPI, EMUODD, 
and RADODD. [deg] 

Array of five elements, the first three of which may be currently 
used. Inclinations to ecliptic of intermediate -target orbits. In- 
put CNEX(i) only when MOPTX(i) = 11. Related to ECIX, OMIX, 
SAIX, SOIX, TPIX, EMUODX, and RADODX. [deg] 

Structural factor, k , in expression (8), 

s 

Propellant tankage factor, k, in expression (7). 

Retro tankage factor, k ,, in expression (11). 

rt 

Efficiency coefficient d in expression (16). Related to BI and El. 
[km/sec] 

Retro engine mass, m , in expression (11). [kg] 

rs 

Ratio of housekeeping power p^ to reference power p^. The 

power transmitted to the propulsion system is that generated by 

the arrays less housekeeping power which is constant along the 

trajectory. The power output of the arrays normal to the sun at 

1 AU is p „ + p, . This option should not be invoked on missions 
ref h 

during which large solar distances are encountered where the power 

developed is less than p, . Erroneous results will be obtained. 

h 

Eccentricity of primary-target orbit. Must be less than unity. 

Input only when MOPT3 = 11. Related to CNI, OMI, SAI, SOI, 

TPI, EMUODD, and RADODD. 

Array of five elements, the first three of which may be currently 
used. Eccentricities of intermediate-target orbits. Input ECIX(i) 
only when MOPTX(i) = 11. Related to CNIX, OMIX, SAIX, SOIX, 
TPEX, EMUODX, and RADODX. 

Efficiency coefficient e in expression (16). Related to BI and DI. 
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EMUODD 

EMUODX 

GAMMAX 

GAP 

HOUR 
EBAL 
ENT PR 

TOUT 


Gravitational constant of primary-target. Input only when 
MOPT3 = 11. Related to ECI, CNI, OMI, SAI, SOI, TPI, and 
RADODD. [m 3 /sec 2 ] 

Array of five elements pertaining to the gravitational constants 
of intermediate -targets. These inputs must be ignored at present. 

Maximum permissible value of the power function y when MODE = 5. 

At solar distances less than the value for which y = GAMMAX, the 
solar arrays are assumed to be tilted such that y is maintained at 
the limiting value. 

Propulsion-corner proximity tolerance-interval, Acr. See discussion 
in the section Avoiding Corners in the Propulsion-time Function. 
Whenever the thrust switch function a grazes the zero-axis within 
the tolerance |Ao"| on any trajectory, an internal counter is in- 
cremented, and the trajectory is considered to be in the neighborhood 
of a propulsion-time corner. Positive value of GAP causes forced- 
thrusting case to be inserted, negative value causes bypass to next 
case, whenever the internal counter reaches the related input 
variable NHUNG. 

Hour-of-day of reference date (e.g. , 17.352D0). Related to 
MYEAR, MONTH, and MDAY. 

Ballistic option indicator. Setting IBAL ^ 0 invokes option 1 dis- 
cussed in the section Ballistic Trajectory Option. 

Indicator which specifies print -length when the iteration in sub- 
routine INTERP fails. Value of 0 causes shortprint and 1 causes 
detailed-print. 

Extra-ecliptic mission indicator. IOUT = 1 or 2 indicates that 
extra-ecliptic target conditions are desired, in which the iterator 
dependent variable triggers Y1 (2) through Y6(2) are set equal 
to 1, and for which the input LAUNCH (which see) should probably 
be set to 1, and parameters related to LAUNCH also set appropriately. 
Ordinarily MOPT2 = 3. No retro stage may be employed. 

= 1 i, e, r specified; f = 0. 

p n 

= 2 i, e, a specified; f optimized. 

n 

In the above, i = final extra-ecliptic inclination, e = final 

eccentricity, r = final perihelion distance, a = final semi- 
P 

(continued on next page). 
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I 

I 

I 



IOUT 
(conto ) 

IRK 

IRL 

IROT 


ISPIN 

ITF 


ITPRNT 

JPP 


JPRINT 


major axis, and f = true anomaly at the final time. Final Q 
and co are optimized in both cases. Related to AE, AR, and 
AAI. 

Numerical integration option (currently not used). 

Primer-origin-proximity step-size-control indicator. Value of 
zero causes the bypass of control, leaving the step-size Au 
constant. See discussion in the section, Integration (Thrust). 

A non-zero value of EROT causes the input ecliptic projection of 
the primer vector and its time derivative to be rotated about the 
z-axis through an angle equal to the difference in longitudes of 
the spacecraft between the last trajectory of the previous case 
(or zero if no previous case) and the first trajectory of the current 
case. This feature permits one to use the initial adjoint variables 
from a 2-dimensional trajectory as the initial-guess inputs for a 
3-dimensional trajectory using the ephemeris option. 

Spinner indicator. Not used at present. 

Provides normal termination conditions for runs which require 
more machine time than is estimated. The value specifies the 
number of machine-time seconds (CPU and I/O) required to exe- 
cute the summary trajectory after halting the iteration-sequence, 
[sec] Does not apply if subroutine REMTIM is dummied. 

Indicator for special print from MINMX3 iterator. Non-zero value 
invokes print. 

Jettison indicator j for electric propulsion sj^stem prior to 
primary-target retro-maneuver, as used in expression (9). 

= 0 Propulsion system not jettisoned 

= 1 Propulsion system jettisoned prior to retro maneuver. 

Unit 11 printout-length indicator. A value of zero causes the 
iterator independent and dependent variables to be output only for 
each summary-trajectory; a value of one causes the same output 
additionally at each iteration of an iteration sequence. 
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JT Jettison indicator for electric propulsion tankage prior to 

primary-target retro-maneuver, as used in expression (9). 

= 0 Tankage not jettisoned 

= 1 Tankage jettisoned prior to retro-maneuver, 

KPART Option for automatically selecting improved independent parameter 

perturbations for generating the iterator's partial derivative matrix. 
The option is invoked by setting KPART = N (N > 0), where N is the 
maximum number of allowed steps, as discussed in the section. 
Perturbation Step Size Selector. KPART must be set back to zero 
if not desired on subsequent cases. 

LAUNCH Launch mode selector, pertaining to the optimization of the de- 
parture asymptote declination, invoiced by LAUNCH = 1. Related 
to X10, Y10, X17, and Y17. 

LOADX Intermediate -target initial-guess feature. Should be used with 

NSET(5) = 1, and then set to zero on the subsequent case. A non- 
zero value of LOADX will invoke this feature, whereby the primer 
A and its derivative A will be loaded into the iterator independent- 
variable arrays at each intermediate-target provided that the trigger 
of the independent variable is on. The sole purpose of this capa- 
bility is merely to generate an initial-guess for a multiple-target 
mission, where the values loaded into the iterator arrays represent 

c P 

continuous A and A at each target. 

MAXHAM Maximum number of times that the program will print the warning 
message BAD HAMILTONIAN on any given computer run. 

MBOOST Launch vehicle selector. 

= 0 ATLAS (SLV3X) /CENTAUR 

1 TITAN III C 

2 TITAN III C (1207) 

3 TITAN III X/CENTAUR 

4 TITAN III X (1207) 

5 TITAN III X (1207)/ CENTAUR 

6 SATURN IB/LM 

7 SATURN IB/CENTAUR 

8 SATURN IC/SIVB/ CENTAUR 

9 TITAN III X (1205)/ CENTAUR 

10 TITAN III B (CORE)/CENTAUR 

11 TITAN in D (1205)/CENTAUR 
(continued on next page). 
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MBOOST 

(cont) 

MDAY 

MODE 

MONTH 

MOPT 

MOPTX 


~12 DELTA 

13 TITAN III D 

14 TITAN III D (1205)/CENTAUR/TE364 (2250) 

15 TITAN m E/CENTAUR 

16 SHUTTLE /TRANSTAGE 

17 SHUTTLE/DELTA 

18 SHUTTLE/AGENA 

19 SHUTTLE/CENTAUR 

20 SHUTTLE/CENTAUR/BURNER II (2300) 

NEG Use input booster coefficients Bl, B2, and B3. 

Day-of-month of reference date (e.g., 26). Related to MYEAR, 
MONTH and HOUR. 

Power variation option selector. The value of MODE is equal to the 
option-number of the power-curve, discussed in the section, Electric 
Propulsion System (which see). Possibly related to ASOL and 
GAMMAX. MODE = 1 has been eliminated. 

Month-of-year of reference date (e.g. , 8). Related to MYEAR, 

MDAY, and HOUR. 

Ballistic option indicator. Using MOPT invokes option 2, discussed 
in the section. Ballistic Trajectory Option, as follows: 

o 

= 0 No action (use input A , A , and v ). 

o o 00 o 

= 1 Generate ballistic solution with flyby end conditions. 

= 2 Generate ballistic solution with orb iter end conditions. 

Related to REVS. 

Array of five elements, the first three of which may be currently 
used. This array specifies the target-number, or planet-number, 
of the successive intermediate -targets, and a value of zero indicates 
absence of the intermediate -target. A zero-entry must not precede 
a non-zero entry. Planet selection is the same as for MOPT2. 
MOPTX(l) pertains to iterator parameters X41-X50 and Y41-Y50; 
MOPTX (2) pertains to X51-X60 and Y51-Y60; and MOPTX (3) pertains 
to X61-X70 and Y61-Y70, Times at the targets are X48, X58, and X68. 
Not to be used unless MOPT2 ^ 0. 
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M0PT2 Launch planet number and ephemeris-option indicator. 

= 0 Analytical planetary ephemeris is not used. 


^ 0 Analytical planetary ephemeris is used and the specific 



launch planet is selected as follows: 

= 1 

Mercury 

= 27 

Hebe 

2 

Venus 

28 

Iris 

3 

Earth 

29 

Flora 

4 

Mars 

30 

Achilles 

5 

Jupiter 

31 

Amor 

6 

Saturn 

32 

Hidalgo 

7 

Uranus 

33 

Alinda 

8 

Neptune 

34 

Grigg-Skjellerup (1977)* 

9 

Pluto 

35 

Kopff 

10 

Ceres 

36 

Grigg-Skjellerup (1982)* 

11 

Input Target ** 

37 

Ganymed 

12 

D'Arrest (1982)* 

38 

Ivar 

13 

Encke (1980)* 

39 

Beira 

14 

Icarus (1987)* 

40 

Kepler 

15 

Eros 

41 

Giacobini-Z inner (1985)* 

16 

Geographos (1983)* 

42 

Borrelly (1987)* 

17 

Encke (1977)* 

43 

Tempel II (1988)* 

18 

Encke (1984)* 

44 

Tempel II (1983)* 

19 

Encke (1987)* 

.45 

Tuttle-Giacobini-Kresak 

20 

Halley 

46 

Schaumasse 

21 

Betulia 

47 

Honda-Mrkos-Pajdusakova 

22 

Toro (1983)* 

48 

Giacobini-Z inner (1979)* 

23 

Pallas 

49 

Icarus (1987)* 

24 

Juno 

50 

Toro (1987)* 

25 

Vesta 

51 

Geographos (1987)* 

26 

Astraea 




MOPT3 Planet number of primary target. Planet selection is the same as 

for MOPT2. If ephemeris is not used, MOPT3 is used only for 
retro-stage mass computations. 

MOPT4 Array of ten elements, specifying up to ten post-swingby targets. 

Planet selection is the same as for MOPT2, and a value of zero 
indicates the absence of a post-swingby target. A negative value 
in MOPT4(l) selects multiple ballistic swingbys, rather than a set 
of single swingbys in which case also set MAXHAM = 0. Negative 
values (in absolute value) produce planet selection the same as for 

MOPT2. When MOPT4(l) < 0, the remaining elements of MOPT4(i) 

(continued on next page) 

*Year -value indicates apparition for which internal orbital elements are most accurate. 

**Input corresponding orbit elements (see CNI, CNEX). None are available for the 
launch planet. 
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M0PT4 

(cont) 


MPOW 


MPRINT 


MPUNCH 


< 0 and > 


MREAD 


maybe positive or negative. See the. section, Swingby Continuation 
Anatysis for details and Sample Case H for an example-case. 

Should be used only for primary-target flyby missions. Related 
to T2, MSWING, NS WING and XSWING. 

Flag used in conjunction with the solar array degi'adation option. 
Value of zero results in the optimum orientation of the arrays 
relative to the sun line. A non-zero value forces the arrays to 
an orientation yielding the maximum power achievable at that 
instant. Related to TPOWER. 

Indicator for printing the summary-trajectory (final trajectory of 
a case) as a function of time or for invoking extra printout. 

= 0 Small-size block print at thrust switch points only 
(SWITCH POINT SUMMARY page). 

= 1 Same as = 0, except expands to become a standard print- 
block of parameters for each computed point along the tra- 
jectory, including the trajectory extension controlled by the 
input variable TGO. 

= 2 Same as = 0, except each block contains extra lines consisting 
of target-relative coordinates and target magnitudes. 

= 3 Combination of = 1 and = 2. 

Punched-card and trajectory-tape generation control. 

= 0 No special output. 

1 Punch final values of independent parameters. 

2 In addition, punch selected mission analysis parameters used 
for graphic documentation or other purposes. 

100 Punch trajectory output used with the ASTEA program. The 
absolute value of MPUNCH determines the frequency of tra- 
jectory points output, e.g. , -3 would result in the punching 
of every third integration point. 

101 Trajectory tape output used with the ASTEA program. The 
absolute value less 100 determines the frequency of trajectory 
points output. Related to NTAPE. 

Card input option (iterator independent variables) 

= 0 No special cards input. 

(Continued on next page). 
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MREAD 

(eont) 

MSWING 


MTMASS 


= 1 The independent variables generated by a previous run bj' the 
MPUNCH = 1 or 2 option are input following the NAMELIST 
case, as discussed in the section, Program Output, 

Array of ten elements, used only when running multiple -target 
ballistic swingbys, such that MSWING (i) corresponds to MOPT4(i) 
and selects the type of swingby maneuver, desired at the respective 
swingby target. Used only if MOPT4(l) < 0, The shooting method 
(MINMX3 iterator) is used, and values of -1, -2, or -3 correspond 
to a swingby passage distance initial guess of rp = ® (i„e. , con- 
tinuous heliocentric velocity). Each element MSWING (i) may have 
any of the following values: 

= - 1 Go* directly for unpowered swingby; if and only if it fails, 
go for powered swingby having flight time T2(i) = initial 
guess. 

= - 2 Go directly for powered swingby only, having T2(i) = flight 
time of post-swingby leg. 

= - 3 Go directly for unpowered swingby; then, whether it succeeds 
or not, go for powered swingby having T2(i) = flight time. 

= - 4 Go directly for unpowered swingby, but using initial velocity 
guess loaded into XSWING(j, i), j = 1,2,3, similar to 
MSWING (i) = - 1. 

= - 5 Same as = - 2, except use initial guess as in = - 4. 

♦"Go for" means "attempt to obtain (solution)". 

Related to MOPT4, T2, XSWING, and NSWING. 

Mission-type selector pertaining to the primary target. 

= 0 Flyby mission. 

1 Orb iter (high-thrust retro-maneuver without velocity loss). 

2 Orbiter (high-thrust retro-maneuver with velocity loss). 

3 Specified arrival excess speed v . 

If v =0, rendezvous mission 
»n 

If v >0, controlled flyby mission 
=°n 

No retro-maneuver in either case. 

4 Orbiter (Electric propulsion system performs spiral maneuver. 

Arrival excess speed v must be specified as zero). 

»n 

(continued on next page) 
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MTMASS 

(cont) 

MUPDAT 

MYEAR 

NDIST 

NHUNG 

NORMAL 

NPERF 

NPRINT 


Other parameters which may be related to MTMASS are DMRETR, 
CTRET, RPER, RAP, THRET, SPtRET, JPP, and JT„ 

Flag indicating whether iterator independent variables at end of one 
case are to be updated for use as first guesses of next case. 

= 0 Do not update independent parameters. 

1 Update independent parameters for next case to be those obtained 
at end of iteration on the current case. 

Year of reference date (e.g. , 1982). Related to MONTH, MDAY, and 
HOUR. 

Identification number of celestial body to be used as the reference 
for the communication distance and angle measurement printed in 
the Extremum Point Summary Table . Identification code is the 
same as for MOPT2. 

Maximum number of propulsion-corner-proximity occurrences 
allowed in a given iteration-sequence. Related to GAP. 

Automatic adjoint-variable scaling. 

= 0 No action. 

1 All A and A are scaled such that becomes unity. 

Identification number of end condition that is to be used as the per- 
formance index when employing the direct parameter optimization 
feature (Improve Mode). The identification code is the same as 
the i in the Yi end condition array. 

Print selection flag. Permits selection of amount of printout 
desired on each case. 

= 0 Print only the case summary. 

1 Print switching point summary of final trajectory. 

2 Print MINPUT and case setup. 

4 Print trajectory summary on each iteration. 

8 Print partial derivative matrix each iteration. 

Combinations of options obtained by summing options desired. 

If NPRINT > 15, printout consistent with NPRINT = 0 is obtained, 
(continued on next page). 
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NPRINT If the sign of NPRINT is reversed to negative, the iterator inde- 

(cont) pendent and dependent variables additionally are printed for every 

trajectory which HILTOP generates (including neighboring trajec- 
tories). 

NSET Iteration-sequence control array. 

NSET(1) Not used for input. 

NSET (2) Not used for input. 

NSET (3) Maximum number of iterations permitted in attempting 
to satisfy constraints in satisfy mode. If zero, no upper 
limit imposed. 

NSET (4) Flag indicating whether constraints are to be satisfied 
prior to entering improve mode. 

= 0 Satisfy constraints first. 

1 Proceed immediately to improve mode. 

NSET (5) Maximum number of iterations permitted after entering 
improve mode. Setting NSET (5) = 1 causes iterator to 
be bypassed and computes single trajectory to obtain 
printout. 

NSWING Swingby continuation analysis option indicator. NSWTNG must be 

negative and has the same definition as MSWING (which see); 
NSWING must be used when MOPT4(l) > 0, and may be used when 
MOPT4(l) < 0. If MOPT4(l) < 0 and MSWING(i) - 0, then MSWING(i) 
will be set to the value of NSWING. Related to MSWING, MOPT4 
and T2. 

NSWPAR Iterator independent -variable perturbation-increment control. 

= 0 No action. 

1 Allows the iterator to vary a given independent-variable 
perturbation Ax whenever a neighboring trajectory is 
detected which has a different number of thrust switch points 
than the associated nominal trajectory. Ax is varied until 
the same number of switch points is achieved. 

NTAPE Specifies the unit-number for the ASTEA trajectory tape. Pertains 

to when MPUNCH * - 101. 
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OMI 


OMIX 

POWFIX 

PSIGN 

RADODD 

RADODX 

RAP 

REVS 

RPER 

SAI 

SAIX 


Ascending node angle (with respect’ to vernal equinox) of primary- 
target orbit. Input only when MOPT3 = 11. Related to CNI, ECI, 
SAI, SOI, TPI, EMU ODD, and RADODD. [deg] 

Array of five elements, the first three of which may be currently 
used. Ascending node angles of intermediate -target orbits. 

Input OMIX(i) only when MOPTX(i) = 11. Related to CNIX, ECIX, 
SAIX, SOIX, TPIX, EMUODX, and RADODX. [deg] 

Launch-vehicle -independent (i.e. , no launch vehicle) trajectory 
option in which the value of POWFIX is the spacecraft’s reference 
power, [lew] 

Flag defining the sense of the launch hyperbolic excess velocity 
relative to the initial primer vector. A value of +1. results in 
the assignment of the geocentric right ascension of the excess 
velocity equal to that of the initial primer vector. A value of -1. 
causes the geocentric right ascension of the excess velocity to be 
180 degrees from that of the initial primer. 

Radius of primary target. Input only when MOPT3 = 11. Related 
to CNI, ECI, OMI, SAI, SOI, TPI, and EMUODD. [meters] 

Array of five elements pertaining to the radii of intermediate 
targets. These inputs are not used at present. 

Apoapse distance of capture orbit about primai'y target. 

[planet radii] 

Number of complete revolutions of the ballistic trajectory gene rated 
when the associated input MOPT is used. Must be a positive whole 
number. 

Periapse distance of capture orbit about primary target. 

[planet radii] 

Semi-major axis of primary-target orbit (must be positive). Input 
only when MOPT3 = 11. Related to CNI, ECI, OMI, SOI, TPI, 
EMUODD, and RADODD. [AU] 

Array of five elements, the first three of which may be currently 
used. Semi-major axes of intermediate-target orbits (must be 
positive). Input SAIX(i) only when MOPTX(i) = 11. Related to 
CNIX, ECIX, OMIX, SOIX, TPIX, EMUODX, and RADODX. 

[AU] 
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SOI 


SOIX 


SPIRET 

STATE 


STEP1 

STEP2 

TCOAST 

TDV 


TGO 


Argument of perihelion of primary-target orbit. Input only when 
MOPT3 = 11. Related to CNI, ECI, OMI, SAI, TPI, EMUODD, 
and RADODD. [deg] 

Array of five elements, the first three of which may be currently 
used. Arguments of perihelion of intermediate-target orbits. 

Input SOIX(i) only when MOPTX(i) = 11. Related to CNEX, ECEX, 

OMIX, SAIX, TPIX, EMUODX, and RADODX. [deg] 

Retro-stage specific impulse (pertaining to the retro-maneuver 
at the primary target), [sec] 

Array of six elements containing the Cartesian position and velocity 
components of the primary target. Use only when MOPT2 = 0 and 
the trigger settings of Y1 (2) through Y6 (2) are 0 or 1. [AU, AU/tau] 
(tau = 58.132440991 days) 

Thrust-phase computation step size, Au. Related to AN. 

Coast-phase computation step size, AjS. 

Array of twenty elements, consisting of the durations of the coast 
phases corresponding to the coast-phase start-times input in the 
associated array TOFF, [days] 

Time of occurrence of an impulsive deep space burn, in days from 
the start of the trajectory, which may be used only if the entire 
trajectory is ballistic (i.e. , electric propulsion is not permitted 
with this option, nor is a third intermediate target). Iterator in- 
dependent variables X64, X65, and X66 must be turned on, as 
these are used as the AV vector components of the deep space 
burn in EMOS. Also, set MAXHAM = 0. The following special 
feature is available regarding a first intermediate-target. If 
1.D5 < TDV < 2.D5, then the burn occurs (TDV - 1.D5) days 
after passage of that target; if TDV > 2.D5, the burn occurs 
(TDV - 2.D5) days before passage of that target, [days] 

Ballistic trajectory-extension print option. When zero, no action. 
When positive, TGO = the number of days that the trajectory is to 
extend ballistically beyond the primary-target when no swingby- 
continuation is requested, and ballistically beyond the (last) post- 
swingby target when swingby-continuation is requested (in addition 
to the post-swingby trajectory segment itself). Any negative value 
will invoke printout of only the post-swingby trajectory segment or 
segments when swingby-continuation is requested. Applies also to 
trajectories with multiple swingbys. [days] 
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THRET 

TOFF 

TP I 
TP'IX 

TPOWER 

TSCALE 

T2 

XANG1 

XANG2 

XSWING 


Retro-stage thrust, f , used only when MTMASS = 2. [lbs] 

Array of twenty elements, consisting of the times, in days from 
the start of the trajectory, at which imposed coast phases are to 
begin. Times must be in ascending order. Related to TCOAST. 
[days ] 

Time from reference date (MYEAR, etc.) to perihelion passage, 
for the primary target. Input only when MOPT3 = 11. Related 
to CNI, ECI, OMI, SAI, SOI, EMUODD, and RAD ODD. [days] 

Array of five elements, the first three of which may be currently 
used. Times from reference date (MYEAR, etc.) to perihelion 
passages, for the intermediate targets. Input TPIX(i) only when 
MOPTX(i) = 11. Related to CNIX, ECIX, OMK, SALX, SOIX, 
EMUODX, and RADODX. [days] 

Solar-cell degradation characteristic-time; nuclear electric pro- 
pulsion radioactive-decay chai'acteristic-time. Related to MPOW. 
[days ] 

Iterator dependent -variable tolerance-interval scaling factor; 
scales all tolerances multiplicatively by the amount TSCALE. 

Array of ten elements consisting of initial estimates of swingby- 
continuation trajectory-segment flight-times, i.e., T2(i) corres- 
ponds to MOPT4(i). [days] 

Latitude of the launch site. Used only if LAUNCH is non-zero. 
Related to XANG2. [deg] 

Maximum parking orbit inclination permitted by range safety 
considerations. Used only if LAUNCH is non-zero. Related 
to XANG1. [deg] 

Array of velocity vectors consisting of initial velocity'- guesses of a 
given post-swingby trajectory segment. Used only when either 
NSWING or MSWING has a value of -4 or -5. See especially the 
description of MSWING = -4. Velocity consists of exactly the same 
values as found in the VI, V2, V3 locations of the trajectory block 
print (first block). Related to MSWING, NSWING, MOPT4, and T2. 
[AU/tau] 



xo 


Array of seven elements, the first six of which contain the Cartesian 
position and velocitjr components of the launch planet. The seventh 
element is not used for input. Used only when MOPT2 = 0. 

[AU, AU/tau] 


The following describes the iterator independent and dependent variable 
arrays of the boundary value problem. Input pertaining to the individual inde- 
pendent parameters is contained in the arrays XI through X70. The independent- 
parameter arrays have five elements for each variable, as follows (where i = 

1, 2, 3, .... 70): 

Xi(l) Input value of parameter. Must be input regardless of trigger 

setting. If trigger is on (i.e. , Xi(2) = 1), input value is used as 
initial guess of independent parameter and is varied at each subse- 
quent iteration. If trigger is off, the parameter is not used as an 
independent parameter and is not changed. 

Xi(2) Trigger indicating whether parameter is to be an independent 

parameter in boundary value problem. 

Xi(2) = 0 Not an independent parameter. 

(Trigger is ’'off"). 

1 Use as independent parameter. 

(Trigger is "on"). 

Xi(3) Maximum change to parameter permitted in a single iteration. 

Should be a positive quantity. Used only if trigger is on. Units 
are same as that of the parameter. 

Xi(4) Perturbation increment used to compute partial derivatives by 

finite differences. Used only if trigger is on. Units are same 
as that of the parameter. 

Xi(5) Weighting factor. Should be a positive quantity. A value of 1, is 

generally recommended. The larger the weighting factor, the more 
the parameter is inhibited from varying. Used only if trigger is on. 

The independent variables are as follows: 
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XI 

A o< 1) 



X2 

A o (2) 


> Initial primer vector. 

X3 

A 0 <3> 

J 


X4 

A 0 d) 



X5 

A o (2) 


^ Initial primer derivative. 

X6 

A o (3) 

J 


X7 

X 

VO 


Initial mass-ratio adjoint-variable. 

X8 

X 

T 


Propulsion-time adjoint-variable. 

X9 



Not used. 

X10 

6 


Geocentric declination of launch hyperbolic 


excess velocity, [deg] 

There is no conversion from input to internal units for any of the adjoint variables. 

2 

Xll Reference thrust acceleration, g„ [m/sec ] 

X12 Electric propulsion system jet exhaust speed, c. 

[m/sec] 

X13 Launch hyperbolic excess speed, v 

[m/sec] 0 

X14 Hyperbolic excess speed at primary target, v 

[m/ sec ] n 

X15 Initial time, t , measured from the reference 

date (MYEAR, etc.), [days] 

X16 Time at the primary target, t measured from 

the reference date (MYEAR, etc.). 

[days ] 

X17 Launch parking orbit inclination, i. Used only 

if LAUNCH = 1. Optimized internally by the 
program if both X17 and Y17 triggers are 
off. [deg] 
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XI s 

X 



o 

Initial spacecraft heliocentric velocit}'. Not 

XI 9 

y 0 

required unless one of the three triggers is on. 


o 

[AU/tau] (tau = 58.132440991 days) 

X20 

z 

0 


X21 


Constant thrust cone -angle, 0. Non-zero value 
invokes the constant-0 constraint. 0 < 0 5 180°. 
Zero-value implies that 0 is optimized along the 
trajectory (variable 0). [deg] 


X22 through X29 are currently not used (although some locations following X21 
are reserved for additional constant thrust cone-angles). 

X30 A Degradation-time adjoint -variable. 

S 

X31 through X40 are currently not used. X41 through X50 pertain to the first 
intermediate target, X51 through X60 pertain to the second intermediate target, and 
X61 through X70 pertain to the third intermediate target. The corresponding 
intermediate -target parameters are ignored if the intermediate target is absent. 
Subscripts 1, 2, and 3 pertain to the first, second, and third intermediate 
targets, respectively. 

X41 
X42 
X43 
X44 
X45 
X46 
X47 
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X48 Time at the first intermediate tax-get, t^, 

measured from the reference date (MY EAR, 
etc.). [days] 

X49 Sample-mass factor, k , , for sample- 

^ samp 1 K 

retrieval at first intermediate target* 

X50 Drop-mass factor, k . , , for instrument - 

1 drop 1 

package dropoff at first intermediate target. 

The independent variables X51 through X60 and X61 through X70 are 
identical to X41 through X50 except that they pei’tain to the second and third 
intei'mediate targets, respectively. A third intermediate target ma 3 ' not be 
present when simulating ballistic missions having a deep space burn (See TDV), 
in which case XG4, X65, and X66 are used as follows: 


X64 

Ax 

X65 

o 

Ay 

X6G 

Az 


> Deep-space velocity-increment. [AU/tau] 


Inputs pertaining to the individual dependent parameters are contained in 
the arrays Y1 through Y70. The dependent -parameter arrays have three ele- 
ments for each variable, as follows (where i = 1, 2, 3, . . . , 70): 

Yi(l) Desired value of the dependent parameter. 

Yi(2) - Trigger. If off (i.e., equal to zero), the para- 

meter is ignored and is not considered a dependent 
parameter. Then the other two inputs pertaining 
to the parameter need not be input. If trigger is 
on, (i.e. , not equal to zero), the parameter is 
considered to be a dependent parameter or con- 
straint. Certain of the parameters may have 
up to three non-zero trigger settings. These will 
be discussed individually below. 

Yi(3) Tolerance of desired value (full interval width). 
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It should be noted that the transversality conditions, which comprise some of the 
parameters, are developed under the assumption that all constraints are of the 
point constraint type. Therefore, the satisfy-mode is sufficient in solving any 
optimization problems for which a complete set of transversality conditions is 
available. 


The dependent-parameter arrays are as given below. T(x) represents 


"the, transversality condition associated with x" and the function T(x) will have 
different values depending upon the constraints imposed on the problem. See 
NOMENCLATURE for definition of symbols and subscripts. 



Y1 

i — i 

< 
i i 

<1 

a[AU] 

Solar distance* [AU] 

T(o> 

Y2 

Ay [AU] , 

n ' 

e 

T(9 t )* 


T < 6 t> 

Y3 

Az [AU] ! 

n 

i[deg] 



T(t„) 

Y4 

Ax [AU/tau] 
n 1 

T(O) 

T(x n ) 


T<4> 

Y5 

Ay [AU/tau] , 
n 1 

T(tO) 

T( y n ) 

( optimal 
flyby 

V 

“O 

Y6 

Az [AU/tau] , 
n 1 

T(f) 

T(z ) 
n 


T(X) 



♦Applicable only for two-dimensional motion in the xy plane. Also requires that 
MOPT2 = 0. 


Under Trigger 1 above, the first set of conditions applies to ordinary targeting 
conditions for position and velocity, and also to extra-ecliptic conditions to be 
satisfied when IOUT = 1; the second set of conditions applies to extra-ecliptic 
missions when IOUT = 2. T{£2), T(to), and T(f) are symbols for the trans- 
versality conditions yielding optimum final node angle, argument of perihelion, 
and true anomaly, respectively. 
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NERVA 



• Trigger 1 

Trigger 2 

Trigger 3 

Y7 

V 

n 

X 

uxx ' 

m x [kgJ 
net ° 

Y8 

T (T) 

T [days] 


Y9 

Current^ not used. 


Y10 

T(6) 

6 [deg] 

Used only if LAUNCH ^ 0. 

Yll 

T(g) 

g [m/sec 2 ] 

p ,[kw] 
ref 

Y12 

T(c) 

c [m/sec] 


Y13 

T (v ) 

v [m/ sec ] 

OOQ 


Y14 

Tfv ) 
°°n 

v [m/ sec ] 

oo n 

extra-ecliptic inclination [deg] 

Y15 

T (t 0 ) 

t [days ] 
o 


Y16 

T (t n ) 

t [days] 
n 

t - t [da3>s ] * 
n o 


*Time transversality with flight time fixed is assigned to Y15 under trigger 1. 


Y17 

T(i) 

i[deg] , 

where i = parking 
orbit inclination. 

Used only if LAUNCH ^ 0 

Y18 

T(x o ) 

x [AU/tau] 
0 


Y19 

T(y o ) 

y ^ [AU/tau] 


Y20 

T ( Z ' 0 ) 

z [AU/tau] 
0 


Y21 

T(0) 

0 [deg] for 0 = 
constant with time. 

Y22 through 

Y29 are currently not used. 


Y30 

T(s) 

s [days ] 

(Degradation time) 


Y31 through Y40 are currently not used, Y41 through Y50 pertain 
to the first intermediate target. 

Y4X Ax^AU] 
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Trigger 1 


Trigger 2 


Y42 

Ay^ [AU] 


Y43 

Az l CauJ 


Y44 

Ax^ [AU/tau] 1 

T(x x ) 

Y45 

Ay^ [AU/tau] ’ 

T(y x ) 

Y46 

Az^ [AU/tau] ' 

T(z x ) 

Y47 


v „ [m/sec] 
“1 

Y48 

T(t x ) 

t [days] 

Y49 

m „ [kg] 

samp 1 


Y50 

m Dkg] 

drop 1 



optimal 

flyby 


Y51 through Y60 and Y61 through Y70 are identical to Y41 through 
Y50 except that they pertain to the second and third intermediate targets, 
respectively. 

3. Default Values of Input Parameters. The following is a complete, 
alphabetical list of the default values of program input quantities having non-zero de- 
fault values, except for the iterator arrays. All other inputs are zeroed. The de- 
fault values of the iterator arrays Xi(l), Xi(2), Yi(l), and Yi(2), for i = 1,2, 3, * . . , 70, 
are zero, and the default values of Xi(3) through Xi(5) and Yi(3) for the same 
range of i are listed in the listing of program inputs of Sample Case H. 


ALPIiAA 

15. 

IRK 

1 

ALPHAT 

15. 

ERL 

1 

AN 

1.5 

ITF 

3 

AR 

1. 

MAXHAM 

5 

BI 

.76 

MDAY 

1 

CTANK 

.03 

MODE 

4 

CTRET 

1/9 

MONTH 

1 

DI 

13. 

MOPT3 

10 

GAMMAX 

1. 

MUPDAT 

1 

GAP 

.0001 

MYEAR 

1975 

HOUR 

12. 

NDIST 

3 
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NHUNG 

25 

STATE (1) 

1 . 

N PR I NT 

7 

STATE (5) 

1 . 

NSET (3) 

300 

STEP1 

.03125 

NSET{5) 

300 

STEP2 

.125 

NSWPAR 

1 

TDV 

- 1 . 

NT APE 

17 

TGO 

- 1 . 

POWFIX 

- 1 . 

THRET 

400. 

PSIGN 

1 . 

TOFF 

20*-l. 

RADODD 

1 . 

T POWER 

10. **30 

RAP 

38. 

TSCALE 

1 . 

RPER 

2. 

T2(i) 

50*i 

SAI 

1 . 

X0(L) 

1 . 

SPIRET 

300. 

X0{5) 

1.0000015 
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C . PR OGR A M OTJ TPUT 

Program output under normal termination conditions provides a listing 
of the program inputs, a. description of the iterator independent and dependent 
variables, the iteration history, the thrust switching history (which optionally 
expands to become the trajectory block print), an iterator summary-page des- 
cribing the results of the iteration, the table of extrema of selected trajectory 
functions, the mission schedule giving calendar dates and target positions and 
velocities, and the performance summary page describing the spacecraft masses. 
These output components are best understood by perusing the various sample 
cases which are included in this document. 

The listing of the program inputs consists of a listing of the iterator in- 
dependent and dependent variable arrays, followed by a presentation of the re- 
maining program inputs, which is alphabetical except that floating-point quantities 
precede integer quantities. 

The page describing the iterator independent and dependent variables con- 
tains only those X and Y parameters (from the listing of program inputs) having 
non -zero trigger values, as these are the parameters relevant to the boundary 
value problem to be solved., The index of each independent and dependent variable 
is given, together with the value of each independent variable, the desired value 
of each dependent variable, the independent -variable iteration-step limits and 
neighboring-trajectory perturbation-increments, the dependent -variable acceptabi- 
lity-tolerances, and the iterator weights. 

The iteration history consists of a summarj^-print of the iterator independent 
and dependent variables at each iteration-step, and also the thrust switching times, 
spacecraft masses, and propulsion system parameters. Depending on the value of 
the input variable NPRINT, the iterator partial derivative matrix -dy./Bx. may 
be printed, and also a summary-print of every trajectory which the program gener- 
ates, including the neighboring trajectories, may be obtained. In the partial 
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derivative matrix, are the relevant dependent parameters, x. are the 
relevant independent parameters, and i is the row index,, 

The thrust switching history expands to become the trajectory block- 
print according to the value assigned to the input variable MPRINT. The indi- 
vidual components of the block-print are described in the section, Auxiliary 
Computations. The block-print occurs at each compute-step, which generally 
corresponds to fixed increments in the trajectory independent variable, u on 
thrust arcs and & on coast arcs. This is not available as the iteration proceeds 
but only on the final, summary-trajectory of each case. 

The iterator summary page displays all of the iterator independent 
variables (at the end of the current iteration), since many of them affect the 
trajectory computations even when they are held fixed. Also displayed Is the 
"switch-count history", which is a listing of the total number of thrust switching 
points on the nominal trajectory of each iteration-step, and the number of both 
thrust -phase and coast -phase computation steps for the current summary-tra- 
jectory, which is the final trajectory at the end of the iteration-sequence for the 
given case. 

The entries in the table of extrema of selected trajectory functions are 
described in the section. Auxiliary Computations. The mission schedule gives 
the position, velocity, solar distance, and ecliptic latitude and longitude of 
the spacecraft and each target at the target intercept or rendezvous time, and 
also the two-body transfer angle between the launch planet and each target. 

The performance summary page includes a mission-type and launch 
vehicle description, electric propulsion system structural parameters and 
mass summary, extreme trajectory and performance conditions, launch and 
primary-target hyperbolic excess speeds and, if applicable, the high thrust 
retro maneuver and capture orbit summary pertaining to the primary target. 
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I 

A print option is available which allows the extension of the summary- 
trajectory ballistically beyond its normal endpoint, which is useful for determining 
where the spacecraft goes after the primary mission objectives have been accom- 
plished. This print option is controlled by the program input quantity. TGO, and 
consists of the trajectory block-print, extremum table of selected functions, and 
an additional mission schedule entry. This print option is displayed automatically 
whenever the primary-target swingby-continuation-analysis is requested. 

HILTOP optionally provides punched-card output for each case under 
normal program termination conditions. When the program input variable 
MPUNCH = 1, the punched-card output consists of seven cards for each case, 
containing all seventy of the iterator independent variables consecutively, each 
card containing ten independent variables in A8 format, 10A8. It is these 
seven cards which may be used to initiate a trajectory on a subsequent run, by 
using the program input variable MREAD = 1. When the program input variable 
MPUNCH = 2, five additional cards are punched for each case, and these cards 
contain the pertinent information which summarize a converged, optimal tra- 
jectory. No cards are punched for a given case if a trajectory or iterator error 
condition exists for that case, including convergence failure or obtaining maxi- 
mum allowable iterations. The contents of the five summary cards (cards 8 
through 12), having units identical to the program input units, are as follows: 


Card 8 Format 

power function indicator, MODE 12 

numerical integration indicator, IRK EL 

launch vehicle indicator, MBOOST 12 

ballistic solution indicator, MOPT 11 

day of month, MDAY 12 

month of year, MONTH reference date 12 

years since 1900, MYEAR - 1900 12 


160 



Card 8 (cont) 

Format 

o 

A, A rotation indicator, IROT 

EL 

retro maneuver propulsion system jettison 
indicator, JPP 

EL 

retro maneuver tankage jettison indicator, JT 

El 

launch planet indicator, MOPT2 

12 

primary target indicator, MOPT3 

12 

mission-type indicator, MTMASS 

12 

ballistic solution indicator, EBAL 

EL 

blank 

EL 

launch mode indicator, LAUNCH 

n 

star-sighting indicator, ISTAR 

12 

spacer 

50x 

case number 

12 

card identifier 

12 

Card 9 

Format 

specific mass, a = a + (L + Ap) a 

t 2. 

A4 

tankage factor, k^ 

A4 

structure factor, k 

s 

A4 

retro engine mass, m^ 

A4 

retro tankage factor, k^ 

A4 

primary-target orbit periapse distance, 

A4 

primary-target orbit apoapse distance, r^ 

A4 

retro stage thrust magnitude, f 

r 

A4 

retro stage specific impulse 

A4 

hour of the day (reference date) 

A4 

extra-ecliptic final inclination 

A4 

efficiency coefficient, b 

A4 

efficiency coefficient, d 

A4 
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Card 9 fcont) 

Format 

efficiency coefficient, e 

A4 

launch vehicle coefficient, 

A4 

launch vehicle coefficient, b 

2 

A4 

launch vehicle coefficient, b^ 

A4 

extra-ecliptic perihelion distance 

A4 

extra-ecliptic final eccentricity 

A4 

case number 

12 

card identifier 

12 


Card 10 Format 

input primary-target position and velocity 6A4 

input initial state and mass ratio 7A4 

reference power in no-launch -vehicle mode, p „ A4 

ref 

thrust-phase computation step size, Au A4 

coast-phase computation step size, A/3 A4 

unused location A4 

input primary-target gravity constant, (i, A4 

v 

input primary-target radius A4 

case number 12 

card identifier 12 


Card 11 Format 

input primary-target semimajor axis A4 

input primary-target eccentricity A4 

input primary-target inclination A4 

input primary-target ascending node A4 

input primary-target argument of perihelion A4 

input primary-target perihelion time A4 

exponent in step-size law, n A4 
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Card 11 (eont) 

Format 

launch- site latitude, L 

A4 

range safety limit, i 

max 

A4 

launch asymptote declination, 6 

A4 

launch parking orbit inclination, i 

A4 

primary-target communication distance, r 

c 

A4 

primary-target communication angle, a 

0 

A4 

retro-stage structure and tankage mass, m 

rsi 

A4 

retro-stage propellant mass, m^ 

A4 

retro-maneuver incremental velocity, Av 

A4 

retro-maneuver velocity loss 

A4 

final solar distance 

A4 

power at primary target 

A4 

case number 

12 

card identifier 

12 

Card 12 

Format 

launch Julian date 

A4 

flight time to primary target 

A4 

initial spacecraft mass, m Q 

A4 

power plant mass, m 

ps 

A4 

propellant mass, m^ 

A4 

tankage mass, m 

V 

A4 

structure mass, m 

s 

A4 

reference power, p , 
ret 

A4 

efficiency, r] 

A4 

net spacecraft mass, m 

A 4 

reference thrust (newtons) 

A4 

maximum solar distance 

A4 

minimum solar distance 

A4 
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Card 12 (cont) 


Format 


A4 
A4 
A4 
A4 
A4 
A4 
12 
12 

No information regarding intermediate targets is punched. It is not 
feasible to design a card-punch routine comprehensive enough to meet every- 
body's possible requirements, and therefore the user seriously contemplating 
using punched cards for summary purposes should design his or her own 
punched quantities and formats. 

On option, the program will write the final trajectory on a magnetic tape, 
or punched cards, for input to the ASTEA program (see definition of input para- 
meter MPUNCH). When punched, each trajectory point requires five cards. The 
content and format of these cards are as follows : 


maximum power 
maximum thrust 
electric-propulsion burn time, T 
travel angle, 0 

V 

retro -maneuver burn time, t^ 
primary-target capture-orbit periapse-speed 
case number 
card identifier 


Card 1 TI, WI, PRATIO, LEND, IT HR, where 


(3D15.0, 212) 


TI = the time since the start of the trajectory (days). 

WI = the mass ratio, current mass over initial mass 

(unitless). 


PRATIO = the power ratio, current power available over 
power available at one AU from the sun. For 
nuclear electric spacecraft, PRATIO = 1 
always (unitless). 

IEND = Flag which indicates last trajectory point. 


= 0 not last trajectory point. 

= 1 last trajectory point. 
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Card 2 
(6D12.0) 


Card 3 
(4D15.0) 


Card 4 
(6D12.0) 


Card 5 
(6D12.0) 


ITHR - Flag which indicates thrusting or coasting. 

= 1 Indicates thrusting (and thrust-on switching points). 

= 0 Indicates coasting (and thrust-off switching points). 

RI(I), 1 = 1, 2, 3; RDI(I), 1 = 1, 2, 3, where 

RI = spacecraft position vector at current time-point 

(AU). 

RDI = spacecraft velocity vector at current time-point 
(AU/tau). 

ELI(I), 1 = 1, 2, 3; CSI, where 

ELI = solar radiation pressure unit vector. 

3 2 

CSI = solar radiation pressure coefficient c g (AU /tau ). 

PI(I), 1 = 1, 2, 3; PDI(I), 1 = 1, 2, 3, where 

PI = target planet position vector at current time-point 

(AU). 

PDI = target planet velocity vector at current time-point 
(AU/tau). 

ELAM(I), 1 = 1, 2, 3; ELAMD(I), 1 = 1, 2, 3, where 

ELAM = Primer vector (Lagrange multipliers adjoint to 

spacecraft velocity), in any units. This is merely 
a vector which lies along the thrust vector. 

ELAMD = time derivative of ELAM (any units/tau). 


If the trajectory is output on tape, the contents of the tape are nothing more 
than the contents of the five cards described directly above. The format of each 
tape record is (30 A8) , so that each record consists of 30 words, each of length 
8 in the A-format. The organization of the tape record is as follows: 

TI, WI, PRATIO, ATHR, (RI(I), 1 = 1, 3), 

(RDI(I), 1 = 1, 3), (ELI (I), 1 = 1, 3), CSI, (PI(I), 1 = 1,3), 
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(PDI(I), 1 = 1, 3), (ELAM(I), 1 = 1, 3), (ELAMD(I), 1 = 1, 3), 
AEND, UNIDUM, UNIDUM, UNIDUM 

where the names and their units are given directly above in the description of the 
cards. UNIDUM (universal dummy variable) represents unused words. Finally, 
ATHR and AEND are real variables which replace the integer quantities IT HE. 
and IEND, respectively, which are punched on the cards. 
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D. JOB CONTROL FOR HILTOP 

1. Program Execution . The HILTOP program is stored at the 
Goddard Space Flight Center IBM 360 model 91 computing facility in object module 
form as a member of a partitioned data set on a user disc pack. The job control 
cards sufficient to access and execute the program are as follows : 


// EXEC 
//GO.GVE 


//GO. FT(J7l : 0iil W 
/ /GO. (-'Tlll-'OOl !)i) 
// GU . F T12F 001 1)0 
//GU.DATAli O') * 


LU Ai) EiL !i EG 1 0 II . GO = 2 /) I; K , PAX 1=’ EP=‘tAf !I,C, I LE~ 
L I fl U!) Do N=; i: ; LOADi ZOF I M Xi ! L ) , 0 1 5P=SMiL nC:J = 


^ J GOU T=G , l> C 0 = ( il E CF: 1 =F n , L A E C L = 3 0 , 
s ysuu t , ;j on = ( n f. cn i =f n , i. n eg l = n , 

0 VSOt) T~!{ , D 0 ! ; =( W E CF! i =F :) / LP* EC L= P. 0, 


3 4 0 <g ' 
r.EOFM = 
r.LK3IZE = 30) 
n Lies f L E = 00 ) 
nLK5izr ; = :jo) 



The NAMELIST input data cards follow the last job control card listed above. 


2. HILTQP Program Machine Requirements . When compiled on the 
IBM 360 Model 91 computer at the Goddard Space Flight Center under the Fortran 
H Level 21.7 compiler with the computer optimization level equal to two, the 
HILTOP program occupies approximately 312, 000 decimal bytes in core. This 
includes the core requirements for the following IBM library subroutines which 
the program uses: 


EHCCOMH2 

IHCECOMH 

IHCEFIOS 

IHCEFNTH 

IHCERRM 

IHCETRCH 

IHCFCVTH 

IHCFDXPD 

IHCFDXPI 

IHCFIOS2 

IHCLASCN 


EHCLATN2 

IHCLEXP 

IHCLLOG 

IHCLSCN 

IHCLSQRT 

IHCLTANH 

IHCNAMEL 

IHCUATBL 

IHCUOPT 

REMTIM 
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The program is written almost entirely in double -precision Fortran IV 
using the non-standard Fortran statement IMPLICIT REAL*8 (A-H, O-Z). This 
results in the assignment of an 8-byte word location to each real variable name 
which begins with one of the letters A through H or O through Z, unless the 
name is specifically declared to be of another type. An 8-byte word contains 15 
decimal digits on the machine cited. As in standard Fortran IV, names commenc- 
ing with one of the letters I through N represent integer variables of 4 -byte word 
length. 


The peripheral equipment referenced by the HILTOP program is the card 
reader, assigned to unit 5, the high speed printer, assigned to unit 6, the card 
output, assigned to unit 7, and two arbitrary output devices, assigned to units 
11 and 12. The execution step requirements of the program are a little less 
than 350,000 bytes of Main Core Storage. 
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IV. 


SAMPLE PROBLEMS AND RESULTS 


A list of the necessary program inputs and a copy of the resulting 
printed output are presented in this section for each of eight sample problems. 

The sample problems were selected to display the use of most of the important 
features of the program. The first problem is an orb iter mission exhibiting 
the use of a high thrust retro engine and a launch vehicle with input reference 
characteristics; the second problem is an asteroid rendezvous mission and em- 
ploys two cases to exhibit the indirect and the direct optimization techniques; 
the third problem is a deep space probe that makes use of the two-dimensional, 
open-angle formulation and yields the minimum flight time solution for specified 
net spacecraft mass; the fourth problem is a Jupiter flyby mission that displays 
the ballistic swingby continuation feature, the fixed cone angle feature, the con- 
strained propulsion time feature, and the launch vehicle independent formulation; 
the fifth problem is an Encke rendezvous mission encountering two asteroids en- 
route; the sixth case is an extra-ecliptic mission that displays the housekeeping 
power option and exhibits the effects of high launch asymptote declinations on launch 
vehicle performance and the low thrust trajectory; the seventh case is a comet 
rendezvous mission which includes the effects of solar array degradation due to 
radiation effects; and the eighth and final case displays the HILTOP's powerful 
capability for ballistic mission design and optimization. The specific mission 
chosen is a cometary flyby past G iac ob ini -Z inner followed by a deep space burn 
10 days after passage, a return to and swingby of Earth (unpowered), a second 
swingby of Earth (powered) , and finally encountering the comet Borrelly nearly 
1023 days from launch. The tremendous flexibility for creating imaginative, 
multi-target mission profiles is demonstrated in this example. 

Machine -time quotes are not given for the sample problems because the 
machine-time values are not necessarily representative of the machine-time 
requirements for running the HILTOP program to perform the task of electric 
propulsion mission analysis. A single trajectory on the IBM 360 Model 91 
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computer takes from about .003 minutes of CPU time for "fast" missions sucb 
as direct outer-planet flybys to about .012 minutes for "slow" missions such as 
Mercury orbiters. However, it is not the time-per-trajectory that is important, 
but rather it is the total number of trajectories required in an iteration-sequence 
that counts, and quite often that number is not as small as in the sample problems 
displayed. 
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A. 


MERCURY ORBITER 


The objective is to place maximum net spacecraft mass into a 1.2 by 
22 o 8 radii orbit about the planet Mercury for specified launch date and flight 
time of 510 days. The analytic planetary ephemeris is employed to locate 
Earth and Mercury on the appropriate dates and to evaluate the initial and 
final state vectors. The long flight time for this mission (compared to ballistic 
transfers) is common to the class of optimum electric propulsion trajectories 
characterized by a 2\ revolution spiral about the sun. Solar electric propulsion 
with power ratio Option 4 (i.e. , peak power is maintained near the sun) is 
assumed. The booster parameters are input and represent the TAT (3C)/Delta/ 
TE 364-3. The final capture orbit insertion maneuver is performed with a 
chemical retro stage, and the retro velocity increment computation includes 
the finite thrust velocity penalty. 


uM I 
Xb( 
Yl( 
Yl'l 
itoP 
Ul = 
Xi = 
X£> = 
XI (5 


h’PU 
2 ) = 
2 ) = 
( 2 ) 
T2 = 
5 U i 
~ 4 < 
- 7 . 
= 5. 


T XI ( 2 ) = 
1 • I) ij / X 7 - 
1 . i> !// Y 2 ( 
= 1 . in)/ Y 1 
3 / MOPT3= 
15.33200 
07101)-!, 
0777U-1, 
11)2 .iEk'i) 



The entire electric propulsion system and the electric propulsion propellant 
tankage are jettisoned prior to the retro maneuver. The parameters optimized 
in this case include the reference thrust acceleration (and therefore the reference 
power ) , the jet exhaust speed, and the launch and arrival excess speeds. The 
NAMELIST input data set used to generate this case is reproduced above and 
the complete output obtained with NPRINT = 15 and MPRINT = 0 is displayed 
on the following pages. 
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CERES RENDEZVOUS 


B„ 

The objective of this mission is to rendezvous with the asteroid Ceres 
with maximum net spacecraft mass. The rendezvous is accomplished by simply 
matching the heliocentric position and velocity of Ceres on the arrival date. Ex- 
cept for the launch excess speed provided by the Titan III B (Core)/ Centaur, the 
entire mission is performed b 3 r a solar electric propulsion system, i.e. , no re- 
tro stage is employed. To show how one inputs an arbitrary planet ephemeris 
for a primary target, the orbital elements of Ceres are input directly; however, 
the same results would be obtained for this case by setting MOPT3 = 10. Although 
the launch and arrival dates are left open, the flight time is constrained to 655 
days. The program optimizes the reference thrust acceleration, the jet exhaust 
speed, the launch excess speed, and the launch and arrival dates. Rendezvous is 
achieved by driving the arrival excess speed to zero, and retro-stage mass com- 
putations are bypassed by setting MTMASS = 3. 

The output for two separate cases is shown for this mission. The first 
case depicts the solution as obtained using the indirect method involving the 
satisfaction of the appropriate transversality conditions. Because the inputs 
were very close to the desired solution and because the indirect method generally 
exhibits quadratic convergence characteristics in the vicinity of the solution, only 
four iterations were required for convergence. The second case displays the use 
of the direct optimization capability of the program. This method ignores the 
transversality conditions and attempts to improve the performance index directly 
while maintaining approximate satisfaction of the specified end conditions. Note 
that this case is simply added behind the first case on a single run submittal. 

Using the value of MUPDAT = 0 assures that the independent variables of the 
two cases are identical. Starting from the same initial conditions as case 1, 
the direct optimization procedure requires 33- iterations to converge. 
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This example is somewhat unfair to the direct method in that it only 
shows the input setup; it shows neither the strong convergence property of the 
direct method when the inputs are far from the solution nor the flexibility of the 
method in choosing the performance index. It does show, however, that con- 
siderable machine time is required once the iteration approaches the desired 
solution. The inputs for the two cases are listed below, followed by the program 
output for the two cases resulting from the settings of NPRINT = 3 and 
MPRINT = 0. 


Ai l I MPUT Xl( 2 ) =1. DO, X2( 2 )=1. DO, X3( 2 )=1. DO, X4( 2 )=1. DO, X5( 2) =1.0.0 

Xu (2) =1.1)0, X11C 2)=1.D0, X12( 2) =1.00, XI 3 ( 2 ) =1. U 0, X 1 ii ( 2) =1. !)•:}, XHj( 2)=1.D0 

Y 1 ( 2 ) = 1 . D 0 , V 2 ( 2 ) = 1. DO, Y3 ( 2 ) =1 . DO, Y4 ( 2 ) =1. i»0, Y 5 ( 2 ) =1 . 00, YG ( 2 ) =1. D C 

Yll( 2 ) =1 . DiJ, Y12 ( 2 ) =1 . DO, Y 13 ( 2 ) =1 . D 0, Y15( 2 ) =1 . DO, Y1-; =fi 55 . DO, 3. H0 

N Pit I NT = 3,iUi PDAT=0, LAUf-JC!l«0, NO ill 1 A L = 1 , M D 0 0 S T = 1 0 , A L P ! I A A = 1 2 . 5 D 0 , M 'C- 1 A G C = 3 


ALPMAT=12 . 5D0 


MO PT2 =3 , : It) PT3 = 1 1 , i 1 Y EAR = 10 7ii , i iO NT! I = 1 1 , ! ID AY =10, i 10 ij ii=C . 5 2 3D 0 , X 7 = 
SAI = 2. 7075D0, EOt =. 0750!) 0, ON I =10 . G 071)0, Oil I =00. 5 1374 SP0, GO 1=71. 
XI =d . 77 in -1, X2 =G . 5S 3D-4 , X3 = - 2 . 5330-1, X4=-l. 14 2D-1, X5=4 . 0030-1 
XG = 3. 3 740-2, XI 1=4 . 4 37D-4 , X12 = 3 . 13404, X13 = l . 54 GO 3, X15 — 5. ID 2 
XI (3 = 1. 440!) 2 Gt END 

Ail I l.'PiJT Yll(2) =!). DO, Y12(2)=0.00,Y 13(2) =0.0 0,Y15(2) =0.0 0 
Y7 = l. 03, 3.D0,0. l)0.,IJPEi(F = 7 ft END 


. 20 GOO 
3520 1320 0 
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PROGRAM INPUTS 
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i 17 TRAJECTORIES WITHOUT PARTIAL DERIVATIVES AND 33 TRAJECTORIES WITH PARTIAL DERIVATIVES REQUIRED FOR THIS CASE 
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CASE 2 ITERATOR PARAMETERS 
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EX1REHUH POINTS OF SELECTED FUNCTIONS 
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OtP OECU iOEGI P*RK INC (OEGJ OEP VINF (M/SEC) C3 <KP**2/SEC*+2) APR VINF <M/SEC> C4 { KH* ♦ 2/SEC**2 ) 

-18,7452 26*5000 1646*95302 2*391626 0*95606 0*000001 


C. DEEP SPACE PROBE 

The objective of this mission is to transfer a specified net spacecraft 
mass from the Earth to a heliocentric distance of 10 AU in minimum time. 

The case demonstrates the use of the two-dimensional, open-angle transfer 
end conditions. A constant power profile consistent with nuclear electric pro- 
pulsion is assumed and the Titan III B (Core)/Centaur launch vehicle is employed. 
The velocity at the final distance is left unconstrained, and the travel angle, 
reference thrust acceleration, jet exhaust speed, and launch excess speed are 
optimized. Convergence is obtained after two iterations. The program inputs 
required for the case are listed below and are followed by the output resulting 
from using the default values of NPRINT and MPRINT. 


&-I INPUT X1(2)=1.!)i) / X2(2)=1.OC/X1»(2)=1.DO/'{5(2)=1.DO / X11(2) e 1.D0 / /{7 = 1..9O 
X12(2)=l.»)0 / X13(2)=l.r)O / Xifj(2)=l.P0 / Yl = l.[)l,2.DQ/Y2(2)=2.nQ / Y!‘(2)=2.DO 
Y5{ 2) =2 . 00, Y7=4 . 02, 3. DO, Yll( 2) =1 .00, Y12( 2) =1 . 00, Y13( 2 ) =1 . DO, LAOUCM=0 
i-iOGOC T=10/i IU PT3 =11,(100 E= 3, XI . 4 3290" 1, X2=2 . 0500/ X4 =- 2 . 5-35500 

X5=4 .4 3190-1, XI 1 = 3 . 2 5100-4, X 12=4 . 35 5:39904 , X13=l. 5239 703, X 10=1. 103 3EN0 
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CASe 1 (CONVERGED) PERFORMANCE SUMMARY 
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Do 


JUPITER FLYBY WITH BALLISTIC SWINGBY CONTINUATION 


The objective of this mission is to place maximum payload on a flyby tra- 
jectory past Jupiter and to calculate for that optimum flyby trajectory the four 
appropriate sets of Jupiter flyby conditions which would result in ballistic con- 
tinuation trajectories to Saturn, to Uranus, to Neptune and to Pluto. In addi- 
tion to the ballistic-swingby-continuation feature, this case demonstrates the use 
of the launch -vehicle -independent option, which is invoked by setting POWFIX to 
a positive number, the value of which is the reference power, p of the space- 
craft, in kilowatts. The specific case presented here is a 550 day Earth-to- 
Jupiter transfer with a launch excess speed of 4 km/sec and with a reference 
power of 15 kilowatts. In addition, two constraints are imposed on the solution: 
(L) the thrust cone-angle 0 is constrained to be constant over the thrust period 
and (2) the propulsion time is constrained to be 200 days. The parameters 
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optimized for this case include the reference thrust acceleration, the thrust cone- 
angle, and the launch date. The velocity at Jupiter arrival is left open. This 
particular case required 6 iterations to converge. The inputs for the case are 
listed above and the program output resulting from setting NPRINT = 3 and 
MPRINT = 0 is displayed on the following pages. The pertinent ballistic swingby- 
continuation trajectory data are displayed on the pages following the PERFORMANCE 
SUMMARY page. 
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If TRAJECTORIES WITHOUT PARTIAL DERIVATIVES ANO 6 TRAJECTORIES WITH PARTIAL DERIVATIVES REQUIRED FDR THIS CASE 
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FOR SOLUTION HAVING 11.12 PASSAGE DISTANCE 
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T WO-eCIOY TRANSFER ANGLE BETWEEN EARTH and NEPTUNE IS 240.3072 degrees 
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E. ENCKE RENDEZVOUS WITH DOUBLE ASTEROID FLYBY 

This example exhibits the multiple target mission capability of HILTOP. The 
specific example. consists of an 830 day rendezvous mission to Encke, arriv- 
ing 30 days prior to perihelion passage in the 1987 apparition, and encountering the 
asteroids Metis and Amherstia at approximately 255 and 551 days, respectively, 
into the mission. The launch vehicle is the Titan III E/Centaur while the solar 
electric propulsion system is characterized by a reference power of 10 kilowatts, 
a specific impulse of 2900 seconds and an efficiency factor of 0.61. Because this 
particular mission results in a nominally high launch asymptote declination, the non- 
coplanar launch maneuver logic is implemented by setting LAUNCH = 1. The asymp- 
tote declination is optimized subject to the constraint that the parking orbit inclination 
be no greater than 32.5 degrees. The constrained declination case was obtained in 
a sequence of fixed declination trajectories, commencing with the unconstrained case. 
In this sequence, it was observed that the initial primer vector was attempting to 
sweep through the North Pole. To permit this to happen, it is necessary to define 
the right ascension of the asymptote to be 180 degrees from the initial primer. 

This is achieved by setting PSIGN equal to -1. Many helpful suggestions for the 
use of HILTOP in studying multiple target missions are given in Reference [l2 3 . 

The inputs to the sample case shown represent converged values for a 
reference power of about 9 kilowatts. From these starting conditions a converged 
solution for 10 kilowatts was achieved after 5 iterations. The complete input 
data set is reproduced below and is followed by the program printout associated 
with NPRINT = 3 and M PRINT = 0. 
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Uii I JJPUT 

Xl<2)=l.DO,X2(2)=l.nO,X3(2)=l.nn,X4(2)=l.n!),X5(2)=l.D!; 
X7=l.:>0 / Xi0( 2)=l. r J0 / Xll( 2)=l.')0 y X12=2.C43n235 ! >4,X13(2) 
X15*-3GO.riO, X1G=- 30.D0 

X4 1( 2) =i . DO, X4 2 ( 2 ) =1.00, X43( 2) =1 .no, X44 ( 2 ) =1 . DP, X4 5( 2 ) 
X4 3=-i> 03 . 00, 1.0 0 

X5i(2)=l.D0,X52( 2)=1.00,X53(2)=l.nn,X54(2)=l.Di;,Xr>5(2) 
X 5 3 =- 3 Oii .00,1.00 

Y1(2)=1.D;),1.L)-5,Y2( 2) =1.0:1, l.D-5, Y3C 2) =1.00,1.0-5, Y4( 
Y 5 ( 2 ) = 1 . 0 0 , 1 . 0 - 5 , Y 0 ( 2 ) = 1 . 0 0 , 1 . 0 - 5 , Y 1 D ( 2 ) =1. DO, l.D-5, Y1 
Y13( 2 ) =1 . DO, 1. D-5, Y41(2)=l. D 0, 1. 0-5, Y42(2)=1.D", l.D-5, 
Y44 ( 2) =2 . 00, i.D-5, Y4 5( 2 ) =2. DO, 1 .D-5, Y4G ( 2) =2. 00, 1 . D- 5, 
Y51( 2) =1,00, l.D-5, Y52( 2) =1. 0 0,1.0-5, Y53( 2) =1.00, l.n- 5 , 
Y5 5( 2) =2 . 00, 1 . D-5, Y5o( 2) =2. DO, 1. 0-5, Y:»3( 2) =1. 00, i. 0-5 
CNI X( I) =5. 53100, cCI :'(!) = . 12 2 O 0 ,.> ! I ;<(1) =0 0 . 7JD0,SD I N(l) 
3 ;M X ( 1 ) = 2 . 3 0 0 3 : ) 0 , T P 1 X ( 1 ) = - 1 :) 0 . D 0 

cm X( 2) =13. 035DO, PCI X( 2) =. 2732DO,OMI X( 2) =330. 11700 , 3J I 
3AI X( 2) =2 . C 700.D0, TPI X( 2 ) =- 723. DO 

M Ti l A3 3 = 3 , ; ID AY =17, MD NT! 1 = 07, M Y EAN= l’J 3 7 , : 10 D E = 5 , : 11 PT2 =3 , 10 
|.;S'»JU3T = 15,ALPifAA=15.!iO,ALP!IAT=15. 235D0, PTAN:C = . 100, •} I =. 
iJP.i 1 NT = 3, LAD NCil = l, XAMfi 1 = 23 . 5D0, XANG2 = 3 2 . 50 0, AAP=- 1. 0-4 
Xl=-3 . 7234C 3 7000 UOD-ui, X2=-2 . 9597-33200000O 00, X3=-3 
X4 =-1.0251733000000 00, X5= 5.0174474000000-01, XP = 3 
XI 0 =“4 . C 70510 7000000 .01, XI 1= 4 . 1 2 5 1 G 4 3 0 0 3 0 3 D - 0 4 , X 1 3 = 


X41=-l. 10 574 J 30000OD 
X44 = 1. 753025300000:) 
X4 o=-i> . 0^2032 bO 00. )C!) 
X53=- 3. 5400020000000 
X5u = I. 3 14 0 7 2 200 COO 1 ) 

< i Eli.) 


,xr.( 2)=1.D0 
=l.nn 

= l.i)0,X40(?)=i.no 
=i.no,X5r.(2)=i.r>n 

2 ) = l . n , i .n-s 
i=i.LU / 5.n:.,i..n-(i 
Y 4 3 ( 2 ) = 1 . D 0 , l.D-5 
Y4 0( 2) =1. D *■, l.D-5 
Y 54(2) =2 . DO, l.D-5 

=4 . 7 onno 

.”( 2) = 255. 331D0 

PT3 = lD / ; !;) PTX = 2* 11 
fili')f),ni =0. DO 
, PG I Hff = - 1 . DO 
. 0227002000000 00 
. 7 0 04 2130 0 0 '■ 0!) - o 1 
.1252153OO0CCD 03 

0 1 , X4 2 = - 7 . 34 G 2 75 7 0 0 00 0[) O0,X4 5 = 1. 7253003000000 51 
00, X4 5=-l . 2 3 5'.) 3 (>500 00 CD O0,X4-r.=-3. ''007527000000 00 

02, X51=-3.3ir.2113!’.C000D OC, X52=-4 . 0113003000000 OP 
00, X54* 1, 0232030000000 30,X55 = 3 . 137194 4 0 0 0 00 - 0 2 
00,XS3=-3. 004 304 *.) 0 0 0 u On 02 
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OEP OECL t DEO | PARK INC (OES) OEP VI NF (M/SEC) C3 (KM**2/SEC*«2 > APR VINF IM/SECI C4 ( KM*«2/SEC**2 » 

•46*1900 32*5003 8002*35326 64*037658 0*08770 0*000000 

AT LIMIT 



F. EXTRA ECLIPTIC MISSION 

The objective of this mission is to place maximum payload in a circular 
orbit of 1.001 AU radius inclined 45 degrees to the ecliptic. The mission is 
2^ years (912.5 days) in duration and departs from Earth parking orbit on % 

April 21, 1979. The solution given in this example is contained in the class of 
six-burn solutions which, for the mission duration assumed, tends to restrain 
the trajectory from deviating far from the nominal 1 AU solar distance. The 
specific case chosen uses the Titan III E/Centaur launch vehicle. The reference 
power (the power delivered to the power conditioners at 1 AU) is 20.35 kw 
and the specific impulse of the thrusters is 3000 seconds. The launch excess 
speed is optimized. The extra-ecliptic end conditions are invoked by setting 
IOUT = 1 and defining the desired values for AR, AE and AAI, the radius, 
eccentricity and inclination, respectively. 

* 

This case exhibits the use of several optional features of the program. 

A total of 0.65 kw of power developed by the solar arrays is reserved for 
housekeeping (non-propuls ive) uses. This option is triggered with the input 
DPOW which is the ratio of housekeeping power to reference power. The power 
delivered to the power conditioners at distances below 1 AU is not permitted 
to exceed that delivered at 1 AU. This constraint is invoked by setting MODE 
equal to 5 and GAMMAX (the maximum permissible value of the power factor 
y) equal to 1. The effects of launch asymptote declination are included in the 
launch vehicle performance model by setting LAUNCH equal to 1. The equatorial 
inclination of the launch parking orbit is limited to a maximum of 36 degrees 
through the input parameter XANG2. Since the geocentric declination of the launch 
asymptote for extra ecliptic missions is usually much greater than this inclination 
limit, the solution will include a non-coplanar injection maneuver from the launch 
parking orbit. The declination of the launch asymptote is optimized. Finally, the 
option of inputting the coefficients of the power profile is illustrated. The inputs 
for this case are listed on the next page. 
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It should be noted that the choice of final orbit radius of 1,001 AU rather 
than 1.0 AU was made to alleviate numerical difficulties arising as a result of 
the corner in the power curve at 1 AU. Neighboring trajectories terminating 
on opposite sides of the corner point tend to possess different partial derivatives 
(i.e. , they will behave differently when subjected to the same perturbation). Con- 
sequently, if the final desired distance were exactly the point of the discontinuity, 
one might expect convergence retardation when the end conditions are nearly 
satisfied. 


AM INPUT Xl(2)*l.DO / X2(2)=l.nO / X3(2)=l.DO,X4(2)=l.nO,X5C2) =1 . DO 

Xn(2)=l.D0,X7=l.D0,X10(2)=l.D0,Xll(2)=l.D0,X12=2.941995D4 

X13(2)=l.D0,X15=3.10370Dl,Xlfi=9.43537D2 

Yl(2)=l.D0,Y2(2)=l.n0 / Y3(2)=l.D0 / Y4(2)=l.D0,Y5(2)=l.D0 / Y6(2)=l.D0 

Y10< 2 ) -1. DO, Yll-20. 35D0y3.no, Y13( 2 ) bI.DO 

LAUNCH = 1 / MB00ST=15 / MTMASS =3, MOD E = 5 

MOPT2 = 3 / MOPT3 = 0 / MYF.AR=19 79,MONTH = 3,HnAY = 21 

BI=.63D0yDI =O.DO / CTAN!C=. 035DOyGAMMAX = l.nO / DPOV/ = 3. 194103194 103D-2 
I 0 UT= 1 , AAI -4 5.DOyARBl.OOinOyXANGlB28.5DOyXANG2B36.no 
ASOL-l.4382DOyO.nfl/-. 223 5D0, 0. DO, -. 214 7D0 
XI — 1. 8492G9016S36D-02, X2 — 3 . 4 44 54 586 96 8 7D- 0 1, X3 = - 5. 292S94680826D 00 
X4 = 7. 2838600516620-01, X5 = 6 . 4 8995 82 7726 3D- 01, Xfi =- 3. 2 7 764 1 8049 73D- 01 
X10—3. S09150729556D 01,X11= 2 . 790574 06 12 56D-04, X13= 5 . 32 8606 53464 5D 03 
A END 
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G. 


COMET RENDEZVOUS MISSION 


The objective of this mission is to deliver maximum payload to the comet 
Tempel n, rendezvousing at perihelion of the comet's path in the 1988 apparition. 
For the case shown, the arrival date is fixed at September 16, 1988; the launch 
date and launch excess speed are optimized to yield maximum payload for the 
nominally 4-year class of solutions. The launch vehicle assumed is the Titan 
III E/Centaur. The electric propulsion system parameters are representative 
of projections of the SERT m spacecraft. Specifically, the reference power is 
8.671 lew, the specific impulse is 2900 seconds, the efficiency is 0.63376 
and the specific propulsion system mass (arrays plus power conditioning and 
thruster subsystem) is 30.285 kg/kw. 

This example illustrates the use of the solar array degradation option. 

This option is invoked by setting the characteristic degradation time TPOWER 

10 

to a value less than 10 . In this case, it is set to 7121* days which means 

that if the arrays were situated at 1 AU and oriented normal to the sun line for 
this amount of time, the power developed by the array would degrade to l/e of 
its initial power output. The input MPOW = 1 forces the arrays to be oriented 
normal to the sun line throughout the mission. Setting the triggers X30(2) and 
Y30(2) to 1 results in the optimal adjustment of the trajectory to accomodate 
the degradation. The complete set of inputs for this case follows and the resultant 
program output begins on the next page. 


&M INPUT Xl(2)=l.O0,X2(2)=l.D0,X3(2)=l.D0,X4(2)=l.DO,X5(2)=l.DO 
X6( 2) =1 .00, X11(2)=1.DO/X13(2)=1.DO/X15(2)=1.DO / X7=1.DO 
Y1(2)=1.00,Y2(2)=1.00,Y3(2)=1.00,Y4(2)=1.00,Y5(2)=1.00,Y6(2)=1.00 
Y13(2)=l.D0,Y15(2)=l.O0,ALPHAT=15. 28 5DO,MBOOST=15,M7?1ASS = 3 
MO PT2 = 3,M0 PT3=4 3,MYEAR = 19 88,MONTH = 9 ,M0AY = 16, H0UR = 1 7. 2200 
X12=2. 84 39 28 504/BI®. 633 7600/01 =0.00, CTANK = . 1D0 / Y11=8.671D0, 3. DO, 1.0-3 
TPOWER® 7121. D0,MPOW=l,X30( 2) =1.00, Y30( 2) =1.00 
Xl= 5. 0729985882240 00, X2= 5. 350 1 836 72 1260 00, X3=-2. 159808142389D 00 
X4=-2. 0693149333120 00, X5= 4.5652477413050 00, X6= 3 . 31 73220885030-01 
Xll= 1. 6479796541510-04, X13= 6 . 718039 798493D 03, XI 5=- 1. 51762 52844 70D 03 
X30=-4. 6574062226460-03 &F.N0 


♦Corresponds to 5% degradation per year. 
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Boor TRANSFER ANGLE BETWEEN EARTH ANO TEMPEL 11(1988) (S 10.9158 DEGREES 



LAUNCH VEHICLE IS TITAN III E/CENTAUR 
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H. MULTIPLE BALLISTIC SWINGBY MISSION 

This computer run demonstrates how the program may be effectively used 
in the investigation of all-ballistic missions (with electric propulsion absent). 

In general, trajectories are forced to be all-ballistic by setting the pro- 
pulsion-time adjoint variable X8 to a ’’large” negative number (X8 = - 1.D3) 
with respect to the mass ratio adjoint variable X7 - 1.D0. These settings (to- 
gether with the reference thrust acceleration Xll = l.D-4, jet exhaust speed 
X12 = 1.D0, and non-zero primer vector components X2 = .IDO at the launch 
planet and X42 = .IDO at the first intermediate-target) guarantee that the 
electric propulsion thrust switch -function will be maintained negative, yielding 
ballistic flight at all times. The electric propulsion specific masses ALPHAA 
and ALPIIAT are also set to zero, so that the net spacecraft mass equals the 
launch vehicle payload. 

The particular ballistic mission simulation demonstrated here involves 
a March 6, 1985 launch of a 1, 635 kg payload by an Atlas/Centaur launch 
vehicle (using the launch vehicle selection default-value of MBOOST = 0; actually, 
changing MBOOST alters the payload mass computation but not the C or the 

O 

trajectory). 

The primary target is the Earth, specified by MOPT3 = 3, and there is 
one intermediate target, the comet Giacobini-Z inner, specified by MOPTX = 41. 
The launch occurs fifteen days before (X15 = - 15. DO) the reference date 
(MYEAR = 1985, MONTH = 3, MDAY = 21), and the spacecraft passes Giacobini- 
Zinner 174.22 days (X48) after the reference date and arrives back at Earth 
353 days (X16) after the reference date. 

Ten days after passing Giacobini-Z inner a deep-space burn (of 188.9 
m/sec) makes possible the re-targeting to Earth. This is accomplished by 
inputting TDV = 1. 00010D5 together with X64(2), X65(2), and X66(2) = 1. DO. 
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The deep space burn velocity -increment (X64, X65, X66) initial guess is the 
zero vector, by default. The deep-space burn could just have well occurred 
before arriving at Giacobini-Zinner (e.g. , TDV = 2.000005D5). The iteration 
sequence consists of a 6x6 hunt, with the independent variables being the initial 
heliocentric velocity at Earth (X18, X19, X20) and the deep-space burn velocitj’’ 
increment (X64, X65, X66), and the dependent variables being the position target- 
ing at Giacobini-Zinner (Y41, Y42, Y43) and at Earth (Yl, Y2, Y3)„ The setting 
MAXHAM = 0 is made to avoid unwarranted BAD HAMILTONIAN warning messages 
due to the presence of the deep-space burn. 

The detailed printout of the iteration sequence is omitted by setting 
NPRINT = 3 (compared to the default value NPRINT = 7), and three extra lines 
are added to each print -block by setting MPRINT = 2. 

Finally, the multiple ballistic swingby option is invoked by means of 
MOPT4 = - 3, 42, which directs the spacecraft to swingby the primary target 
(Earth) in such a manner as to re-target to Earth (MOPT4(l) = - 3), and then, 
at the second Earth passage, to again swingby in such a manner as to target to 
the comet Borrelly (MOPT4(2) = 42). Furthermore, the first Earth swingby is 
specified as unpowered (MSWING(lj = - 1) having T2(l) = 480 days as the initial 
guess of the Earth-to-Earth transfer time, and the second Earth swingby is 
specified to be powered (MSWING(2) = - 5) with a specified transfer time from 
Earth to Borrelly (T2(2) = 127 days) and using the input Earth-departure helio- 
centric velocity initial-guess XSWING(1, 2) = .62D0, .91D0, 0.D0. Both swing- 
by iterations converged, and the powered-swingby incremental speed turned out 
to be -13.7 m/sec, the minus sign denoting a braking burn. The default value 
of TGO = - 1.D0 resulted in the printout of several pages describing the tra- 
jectory segments following the time of primary target passage. 

It should be noted that the trajectory simulation would have terminated 
at the primary target (first Earth encounter) if MOPT4(l) were zero, and would 
have terminated at the second Earth encounter had MOPT4(2) been zero. The 
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mission, as simulated, thus consists initially of a comet-fly 1 ^, followed by a 
second comet-flyby making use of a double Earth- swingby to reach the second 
comet. 


The inputs for this computer run are listed below, and the resulting pro- 
gram output is displayed on the following pages. 


&M INPUT X7=l.O0,X2=.lDO,X42=.in0,X8=-l.D3,ALPHAA=0.nO,ALPHAT=0.D0 

X64(2)=1.PO / X65(2)=1.DO/X66(2)=1.DO / X18(2)=1.DO / X19(2)=1.DO / X20(2)=1.DO 

Yl(2)-l.D0/Y2(2)=l.O0/Y3(2)=l.n0,Y41(2)=l.D0,Y42(2)=l.D0/Y43(2)=T.OO 

M YEAR = 19 85,MOMTH-3,MOAY=21,MOPT2=3,MOPT3 = 3,MOPTX=4 1, Xll=l,D-4 / X12 = l.D0 

TDV=1. 0001005 „ XI 5=- 15.D0, X16=3 53. DO, X4 8=1 74 . 22D0 

Xl8=-8. 3276D-2/X19=-l. 0034D0,X20--2.0143D-3 

12=480.00,127.00/ XSW I NG( 1/ 2) =. 6200, .9100, 0.D0 

M0 PT4 =- 3 , 4 2 , MSW I NG =- 1 , - 5 , M PR I NT = 2 ,MAXH AM = 0, NPR I NT=3 A END 
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COMPUTE STEPS. 52 COAST COMPUTE STEPS 
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CASE 1 (CONVERGED) PERFORMANCE SUMMARY 
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